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In  this  work,  we  analyzed  the  anomalous  behavior  of  a  two  dimensional  interacting 

Fermi  system.  The  first  anomaly  that  has  been  investigated  is  related  to  the  low-energy 

property  of  Fermi  liquids  in  two  dimensions.  The  forward  scattering  processes,  in  the 

second  order,  lead  to  a  log  divergence  in  the  imaginary  part  of  the  self-energy  on  the  mass 

shell.  We  found  that  this  divergence  becomes  even  more  singular  at  higher  orders  and 

takes  the  form  of  a  power-law,  potentially  leading  to  the  breakdown  of  Fermi-hquid  theory. 

We  deal  with  the  divergence  by  re-summing  the  maximally  divergent  terms  at  each  order 

of  interaction,  with  the  result  that  the  renormahzed  imaginary  part  of  self-energy  is  well 

defined  on  the  mass  shell.  We  also  investigated  in  detail  the  contribution  of  the  collective 

modes:  zero-sound  mode  (for  a  contact  potential)  and  the  plasmon  mode  (for  a  Coulomb 

potential)  to  the  nonanalytic  temperature  (T^)  dependent  term  in  the  specific  heat. 

We  found  that  the  zero-sound  mode  does  not  make  any  contribution  to  the  nonanalytic 

T^-dependent  term  in  the  specific  heat,  whereas  the  plasmons  do  make  a  nonanalytic 

contribution. 

The  second  anomaly  is  the  lack  of  experimental  observation  of  a  polarization- 
dependent  effective  mass  as  a  function  of  the  applied  parallel  magnetic  field  in  a  Si- 
MOSFET  sample.  We  investigated  this  experimental  observation  in  the  firamework  of  an 
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TV-component  model,  where  A'^  is  an  additional  degree  of  freedom  (which  in  our  case  is  the 
valley  degree  of  freedom).  This  model  gives  the  change  in  the  effective  mass  between  an 
unpolarized  and  a  fully  polarized  state  to  be  not  more  than  1  to  3%,  a  result  consistent 
with  the  experiment. 
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CHAPTER  1 
INTRODUCTION 

1.1  Survey 

The  physics  of  two-dimensional  (2D)  systems  has  been  the  subject  of  intense  research 
both  experimentally  and  theoretically,  over  the  last  few  decades.  Exciting  physics  such 
as  Wigner  crystalUzation  of  electrons  at  low  electron  densities  [1]  and  the  scaling  theory 
of  localization  [2],  were  first  predicted  and  later  confirmed  by  experiments.  On  the  other 
hand,  novel  phenomena  such  as  the  existence  of  a  metal-to-insulator  (MIT)  transition  [3], 
the  quantum  hall  effect  [4],  the  firactional-quantum  hall  effect  [5],  and  radiation- induced 
zero-resistance  states  [6,  7]  were  discovered  first  and  later  given  theoretical  explanation. 
Many  2D  systems  have  been  described  successfully  in  the  framework  of  Fermi-liquid 
theory.  However,  the  validity  of  Fermi-liquid  theory  in  2D  has  been  a  contentious  issue. 
Fermi-liquid  theory  is  known  to  break  down  in  ID.  Anderson  conjectured  that  even  in  2D, 
Fermi  liquid  is  destroyed  for  arbitrarily  weak  interactions  [8]. 

The  basic  postulates  of  the  Landau  Fermi-liquid  theory  state  that,  low  energy 
properties  of  a  interacting  fermionic  system  can  be  described  similar  to  a  non-interacting 
Fermi  system  [9].  This  phenomenological  theory  makes  specific  predictions  for  the  low- 
energy  dependence  of  thermodynamic  and  transport  properties  on  temperature,  frequency 
and  other  small  energy  or  momentum  scales.  In  particular,  the  predictions  for  specific  heat 
and  spin  susceptibifity  are  as  follows.  Specific  heat,  C{T)  scales  linearly  with  temperature 
T,  for  r      0,  whereas  the  spin  susceptibility  approaches  a  constant  as  T  ->  0.  Over  the 
past  five  decades,  this  theory  has  been  enormously  successful  in  describing  the  low-energy 
properties  of  many  fermionic  systems,  particularly  ^He  and  also  a  large  number  of  metals. 
The  common  features  of  a  metallic  Fermi-liquid  system  include,  the  existence  of  a  Fermi 
surface,  a  specific  heat  proportional  to  temperature,  a  constant  spin  susceptibility,  a  finite 
density  of  low-energy  single  particle  excitations  and  a  conductivity  which  increases  with 

decrease  in  temperature.  This  theory,  though,  has  failed  to  account  for  the  behavior  of 
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quite  a  few  heavy  fermion  systems  and  high-Tc  superconductors.  The  deviation  from 
the  simple  Fermi-hquid  picture  for  the  electrical  conductivity  and  other  features  such  as 
the  thermal  conductivity,  the  nuclear  relaxation  rate  etc.,  in  the  normal  state  of  high-Tc 
superconductors  is  quite  perplexing  given  the  fact  that  many  of  these  materials  have  large 
Fermi  surfaces  [10,  11]. 

Issues  related  with  the  validity  of  Fermi-Liquid  theory  in  systems  which  exhibit 
metal-insulator-transition  have  also  come  under  intense  scrutiny.  Two-dimensional  systems 
can  exhibit  insulating  and  metallic  behavior  depending  on  the  concentration  of  impurities 
and  the  density  of  charged  carriers.  Among  the  many  different  theoretical  approaches  to 
understand  2D  systems  in  either  the  metallic  or  the  insulating  regime,  those  based  on  the 
framework  of  Fermi-liquid  theory  have  been  more  successful,  compared  to  other  theories, 
in  describing  the  transport  and  thermodynamic  properties.  Nevertheless,  deviations  from 
the  Fermi-hquid  universahty  have  been  observed,  such  as  the  absence  of  spin-split  mass 
in  a  spin-polarized  system  and  nonanalytic  terms  in  the  thermodynamic  and  transport 
quantities.  Section  1.2  discusses  nonanalj^icities  in  the  thermodynamic  and  transport 
properties,  and  addresses  the  important  question  regarding  the  validity  of  Fermi-hquid 
theory  in  2D.  Section  1.3  discusses  the  metal-insulator-transition  in  2D.  The  issues  related 
to  anomalous  effective  mass  behavior  have  been  taken  up  in  Section  1.4.2. 

1.2    Deviation  from  Fermi-Liquid  Universality 

Recent  studies  on  the  low  energy  behavior  of  Fermi-liquid  systems  away  from  the 
critical  point,  reveal  deviations  from  the  Fermi-liquid  universality  due  to  the  presence 
of  soft-modes  which  are  responsible  for  the  nonanalytic  terms  in  thermodynamic  quan- 
tities hke  specific  heat  and  spin  susceptibility.  These  nonanalytic  terms  survive  near 
the  quantum  critical  point  and  are  responsible  for  the  breakdown  of  the  Hertz-Millis- 
Moriya  (HMM)  theory  [12],  leading  to  either  a  first-order  or  incommensurate  ordering  (or 
both)  [13].  Thus  a  study  of  the  low-energy  physics  of  Fermi  systems,  which  has  its  own 
inherent  merit,  acquires  extra  significance  due  to  its  impUcations  near  the  quantum  critical 
point  [13,  14,  15]. 

Approaches  to  understand  a  system  of  weakly  interacting  fermions  in  3D  through 
perturbative  techniques  on  a  tractable  microscopic  model  of  short-range  interacting 
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fermions  have  added  to  a  regular  T  term  in  the  leading  specific  heat  term  and  a  constant 
term  for  the  spin  susceptibility  as  T  — >  0,  in  agreement  with  the  results  of  a  Fermi 
gas  [16,  17].  These  results,  when  extended  to  arbitrary  dimensionality  D  >  1,  yield  similar 
results  provided  the  interaction  falls  off  rapidly.  The  similarity  in  behavior  of  the  Fermi- 
gas  result  and  the  Fermi-liquid  result  obtained  via  perturbative  calculations,  stop  at  the 
leading  order  corrections.  Perturbative  calculations  beyond  the  leading  term  in  3D  obtain 
a  nonanalytic     InT  term  [18,  19,  20]  and  in  2D  a      term  [21,  22,  23]  for  the  specific 
heat,  different  from  the  regular      term  as  obtained  for  the  Fermi  gas.  The  non-uniform 
susceptibility,  in  terms  of  the  momentum  Q,  gets  a  nonanalytic      In  Q  term  for  3D  [15] 
and  for  2D  max{Q,T)  [13,  15,  22,  24]  instead  of  a  analytic  max{Q'^,T^)  as  the  next-to- 
leading  term.  Some  of  these  predictions,  especially  on  the  specific  heat  behavior,  have 
been  confirmed  by  experiments.  Specific-heat  measurements  on  bulk  ^He  have  revealed 
beyond  the  linear-in-T  term  an  additional     log  T  term  [25],  a  similar     logT  behavior 
was  also  observed  in  heavy-fermion  UPts  [26] .  The  origin  of  the  nonanalytic  terms  were 
attributed  to  nonpertiubative  spin-fluctuation  processes.  Heat  capacity  measurements 
on  monolayers  of  ^He  adsorbed  on  atomically  flat  graphite  yield  a      correction  to  the 
specific  heat  in  agreement  with  theoretical  predictions  [27]. 

For  D  =  1  the  perturbation  theory  is  singular  and  Fermi-hquid  theory  is  destroyed. 
This  break-down  shows  up  already  at  second  order  in  the  interaction  [28].  Apart  from 
the  low  dimensionality  there  are  other  sources  which  can  lead  to  the  breakdown  of  Fermi- 
liquid  theory.  Using  weak  couphng  analysis  one  can  show  that  for  certain  special  shapes 
of  the  Fermi  surface,  for  example,  a  Fermi  surface  satisfying  the  'nesting'  property  {i.e., 

=  ~^p+q)'      Fermi  liquid  is  unstable.  A  'nesting'  Fermi  surface  can  lead  to  the 
breakdown  of  translational  invariance  by  the  formation  of  a  charge/spin  density  wave 
with  a  typical  ordering  vector  Q.  There  have  been  suggestions  that  the  non-Fermi  hquid 
behavior  of  high-Tc  materials  is  due  to  the  shape  of  the  Fermi  surface  [29].  A  dispersion 
relation  in  momentum  space  which  has  saddle  points  on  the  Fermi  surface  can  also 
affect  the  low  energy  properties  of  fermions  (the  velocity  is  zero  at  the  Fermi  surface). 
A  second-order  perturbation  analysis  reveals  that  the  quasi-particle  decay  rate  depends 
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linearly  on  energy  or,  in  other  words,  is  of  the  same  order  as  the  real  part  of  the  self- 
energy  [29],  signahng  a  non-Fermi  liquid  behavior.  The  validity  of  Fermi-liquid  theory  in 
two  dimensions  has  been  a  source  of  intense  research  since  Anderson  suggested  that  in  two 
dimensions  the  Fermi  liquid  is  unstable  even  for  weak  coupUng  [8].  His  argument  was  that 
scattering  of  two  interacting  particles  at  the  same  point  of  the  Fermi  surface  will  cause  a 
finite  phase  shift.  This  means  that  even  in  the  limit  of  the  lowest  energies,  two  particles 
with  opposite  spin  cannot  occupy  the  same  quantum-mechanical  state,  as  is  assumed  in 
the  quasi-particle  picture  of  Landau's  Fermi-liquid  theory.  A  second  argument  was  that 
the  presence  of  an  anti-bound  state  (pole  in  the  particle-particle  channel)  will  lead  to  an 
instability  of  the  Fermi  liquid.  The  effect  of  phase-space  shifts  and  the  anti-bound  state 
has  been  analyzed  in  detail  with  the  conclusion  that  one  cannot  infer  the  breakdown  of 
the  Fermi  hquid  from  these  effects  [30] . 

In  this  work  we  have  considered  a  2D  fermion  with  a  weak  short-range  interaction 
potential.  A  perturbative  analysis  of  such  a  system  with  a  linearized  single-particle 
spectrum  reveals  a  singularity  in  the  imaginary  part  of  the  self-energy  at  second  order 
in  the  interaction  [22,  31,  32,  33].  This  singularity  is  present  on  the  mass  shell  (e  =  ^p) 
and  is  logarithmic  in  nature.  The  divergence  is  even  more  singular  at  higher  orders 
and  starting  from  the  third  order  onward  goes  as  a  power  law.  An  important  point  is 
that  these  divergences  are  present  even  for  an  arbitrarily  weak  interaction.  A  similar 
divergence  on  the  mass  shell  is  observed  in  ID.  This  divergence  is  due  to  what  is  known 
as  the  infra-red  catastrophe,  where  an  on-shell  fermion  can  emit  an  infinite  number  of 
soft  bosons  consisting  of  spin  and  charge  fluctuations,  leading  to  a  power-law  divergence 
on  the  mass  shell  [28] .  At  each  order  in  the  interaction,  the  divergences  can  be  cutoff  by 
invoking  the  finite  curvature  of  the  single-particle  spectrum,  yet  the  series  itself  remain 
non-convergent  at  low  energies  [34].  These  result  seem  to  justify  Anderson's  conjecture 
about  the  invalidity  of  the  Fermi-liquid  picture  in  2D.  This,  however,  is  not  the  complete 
analysis. 

The  divergence  in  the  self-energy  which  at  each  order  in  the  interaction  is  due  to  the 
interaction  of  fermions  with  the  upper  edge  of  the  particle-hole  continuum,  when  summed 
to  all  orders  corresponds  to  the  interaction  of  fermions  with  the  zero-sound  mode.  A 
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careful  selection  of  diagrams  that  at  each  order  are  the  most  divergent,  and  re-summing 
to  all  orders  in  the  interaction  removes  these  divergences,  thus  restoring  the  Fermi-liquid 
picture.  The  resulting  renormahzed  self-energy  with  the  divergences  removed  acquire 
features  which  indicate  a  deviation  from  the  standard  Fermi-liquid  behavior.  Of  these 
the  most  prominent  is  a  non-Lorentzian  shape  of  the  spectral  function.  The  spectral 
function  acquires  a  shoulder  at  the  threshold  for  the  emission  of  zero-sound  bosons  and  a 
nonanalytic  u\u\  term  in  the  real  part  of  the  self-energy.  In  addition  to  the  contribution  to 
the  u)\uj\  term  from  the  nonperturbative  processes,  the  second  order  term  in  perturbation 
also  contributes  to  the  u\ij\  nonanalytic  term  to  the  real  part  of  the  self-energy.  A 
nonanalytic  c<;|c<;|  term  in  the  real  part  of  the  self-energy,  in  principle,  translates  into  a 
nonanalytic      term  in  the  specific  heat. 

The  nonanalyticities  in  the  thermodynamic  and  transport  (for  a  disordered  system) 
properties  are  due  to  what  are  known  in  a  general  sense  as  the  Generic  Scale  Invariance 
(GSI)  of  the  system  [14].  Correlation  functions  away  from  the  critical  point  generally 
decay  exponentially  and  are  short-ranged  with  a  characteristic  length  or  time  scale.  On 
the  other  hand,  near  the  critical  point  the  correlations  become  long  ranged  and  the  decay 
is  only  power-law.  Correlation  functions  with  long-range  order  exhibit  scale  invariance: 
a  change  of  scale,  e.g.,  in  length,  represents  a  change  in  the  correlation  function  by  an 
overall  factor.  For  cases  where  long-range  order  is  obeyed,  even  for  regions  away  from  the 
critical  point,  the  system  is  said  to  exhibit  "generic  scale  invariance"  (GSI).  In  simple 
terms,  correlation  function  that  are  long  ranged  and  belong  to  the  class  of  GSI  (also  called 
soft  modes)  have  their  Fourier  components  oc  1  /k"  (where  a  is  some  positive  exponent  and 
k  is  the  momentum)  for  regions  even  away  from  the  critical  point.  Correlation  functions 
which  are  short-ranged  have  the  form  oc  1/(A;"  +  m),  where  m  is  a  constant  mass  term. 
One  such  example  is  an  0( A'') -symmetric  Heisenberg  ferromagnet  which  exhibits,  in  the 
ordered  phase,  nonanalyticity  in  the  magnetization  as  a  function  of  the  magnetic  field. 
The  ordered  phase  contains  gapless  Goldstone  modes  and  the  corresponding  correlation 
functions  are  long-ranged  and  belong  to  the  class  of  GSI.  These  are  responsible  for  the 
nonanalyticity  in  the  magnetization.  Other  examples  include  the  nonanalytic  behavior  of 
conductivity  as  a  function  of  frequency.  In  a  disordered  system  the  leading  term  in  the 
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conductivity  is  ctq  =  ne^r/me  (Drude  result,  where  n,e,r  and  rUe  are  the  density,  charge, 
scattering  time  of  electrons  with  impurities  and  the  mass  of  electrons  respectively)  and  the 
leading-order  weak-localization  correction  is  given  by  [36] 


where  D,  q,  Q  and  Np  are  the  diffusion  constant,  momentum,  frequency  and  density  of 
state  respectively.  For  D  =  2  this  leads  to  a  logarithmic  correction  and  for  2  <    <  4,  this 
results  in  a  nonanalytic  frequency  dependence 

^^^^  =  1  +  pn^'^-'y'  (1.2) 

cr(0) 

with  0>O. 

The  thermodynamic  properties  acquire  nonanalytic  corrections  even  for  disorder-free, 
point-like  interacting  fermions.  Although  the  bare  interaction  is  short-ranged,  the  effective 
interaction  at  finite  frequencies  is  long-ranged.  At  second-order  in  the  interaction  U  the 
effective  potential  is  U  =  U'^'n.{Q,q),  where  Il{Q,q)  is  the  dynamical  polarization  bubble, 
which  for  the  momentum  range 


where  ao  is  a  constant  and  l/q  is  the  generic  scale  invariant  term  [22].  The  long  range 
interaction  forces,  prevent  a  regular  expansion  of  the  self-energy  in  powers  of  u  leading 
to  a  nonanalytic  expansion  in  frequency.  The  nonanalyticity  that  arises  in  the  self-energy 
is  transferred  to  the  observable  quantities  thus  resulting  in  non-regulax  subleading  terms 
(Section  3.2.1).  In  a  system  with  disorder,  fermions  interact  with  Priedel  oscillations  (these 
oscillations  are  due  to  the  presence  of  a  static  local  perturbation).  This  interaction  falls 
off  as       with  the  distance  from  the  impurity  and  are  the  source  for  nonanalyticities  in 
transport  properties  (these  are  the  interaction  corrections,  unlike  the  earlier  example  of 
weak  localization  correction). 

Chubukov  and  Maslov  [22]  identified  the  nonanalyticities  in  the  self-energy  for 
D  =  2  with  the  three  scattering  processes,  which  are  essentially  one-dimensional  in 
character.  These  are  the  forward-scattering  process  with  the  incoming  momenta  nearly 


(1.1) 


M/vf  <  9  <  fcf  is  given  by  U  =  U^ao—  +  ■■■  , 


Vfq 
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the  same,  undergoing  small  transfer  of  momentum  and  two  back-scattering  processes 
involving  scattering  of  near-opposite-momentum  quasi-particles,  undergoing  either  a  small 
momentum  transfer  or  large  "2A;f"  transfer.  Second-order  perturbative  analysis  reveals 
a  nonanalytic      behavior  of  C{T)  arising  from  the  back-scattering  processes  ,  whereas 
the  forward-scattering  process,  for  a  linearized  spectrum,  contributes  an  unexpected  log 
singularity  at  the  mass  shell  in  the  imaginary  part  of  the  self-energy.  As  discussed  earlier 
this  divergence  is  the  analog  of  a  stronger  power-law-singularity  in  ID  and  is  dealt  with 
by  re-summing  the  perturbative  series  by  including  the  maximally  divergent  diagrams.  An 
offshoot  of  this  re-summation  is  a  collective  mode.  In  the  past,  the  role  of  collective  modes 
in  the  nonanalyticity  in  the  thermodynamic  potential  had  been  treated  by  assuming 
that  the  collective  modes  are  free.  Indeed,  by  invoking  power  counting  arguments  [21], 
it  can  be  shown  that  the  collective  mode  contributes  to  a      term  in  the  nonanalyticity. 
However,  collective  modes  are  not  free  excitations.  The  contribution  of  the  collective  mode 
to  C{T)  has  been  analyzed  in  detail  in  this  work  (Sees.  3.2.1  and  3.2.2  ).  We  show  that 
the  collective  modes  do  not  contribute  to  the      term  in  the  specific  heat  in  2D. 

In  the  next  Section,  we  will  discuss  the  validity  of  Fermi-hquid  theory  in  2D  systems 
which  exhibit  a  metal-to-insulator  transition. 

1.3   Metal-Insulator-Transition  (MIT)  in  2D 

In  the  late  1970s,  an  interesting  theoretical  viewpoint  emerged  on  the  conductance 
of  a  2D  system.  Following  the  concepts  of  Thouless  [35],  Abrahams  et  al.  [2],  developed 
a  single-parameter  scaling  theory  of  localization  (electron-electron  interactions  were  not 
taken  into  account).  The  scaling  equation  being 

^  (1.3) 

where  g  =  G{L) / {e^ /2fi)  is  the  dimensionless  conductance  and  is  the  scaling  parameter, 
and  L  is  the  typical  length  along  each  of  the  directions.  The  asymptotic  values  of  0{g{L)) 
obtained  via  physical  arguments  can  qualitatively  describe  the  flow  equation.  That  is, 
starting  from  a  given  conductance  G{L)  for  a  system  of  size  L  (volume  L^,  where  D  is 
the  system  dimensionahty),  the  flow  equation  can  in  principle  predict  the  conductance 
(for  example,  when  the  linear  dimension  L  is  doubled).  (Of  course  it  is  imphcit  that  the 
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material  property,  the  density  of  disorder  and  the  electron  density,  do  not  change).  In 
the  large-conductance  {G  »  e^/^)  regime,  the  conductance  shows  Ohm's-law  behavior 
with  G  (X  L^~^,  thus  0  =  D  -  2.  In  the  other  extreme  with  the  electrons  localized,  the 
conductance  gets  exponentially  suppressed  with  increasing  length  and  0  =  —L/^,  where 
^  is  the  localization  length.  It  is  assumed  that  /3  varies  monotonically  between  these  two 
asymptotic  limits. 

The  consequence  of  this  theory  is  that  a  2D  system  with  i      oo  is  an  insulator  at 
zero  magnetic  field  and  zero  temperature.  For  a  finite-sized  sample,  the  scaling  relation 
implies  that  an  insulating  behavior  will  be  observed.  If  the  2D  sample  were  dirty  to  begin 
with,  the  system  would  show  exponential  localization  as  a  function  of  system  size.  If  the 
sample  were  relatively  clean  to  begin  with,  then  the  system  on  further  increase  in  size 
would  show  logarithmic  locahzation  (in  3D,  there  is  a  critical  disorder  separating  localized 
and  metallic  behavior  ).  The  form  of  the  2D  conductivity  showing  logarithmic  behavior  is 
given  by 


The  above  result  obtained  for  zero  temperatures  can  describe  the  conductivity  at  nonzero 
temperatures  with  a  slight  modification.  According  to  Gorkov  et  al.  [36]  and  Thouless  [37], 
at  finite  temperatures,  inelastic  scattering  introduces  random  fluctuations  in  the  time 
evolution  of  an  electronic  state  thus  limiting  the  quantum  interference  necessary  for 
locahzation.  Typical  lengths  (cut-off  lengths)  beyond  which  the  phase  is  not  conserved 
are  given  by  Lrh  =  {DTinY^'^,  where  D  is  the  diffusion  constant  and  Tj„  is  the  inelastic 
hfetime.  (  Altshuler  et  al.  [38,  39]  argued  that  for  a  coUision  event  with  small  energy 
change  AE      rg^ ,  Tin  cannot  be  taken  as  the  dephasing  time;  the  dephasing  time  will  be 

~  (A£'Tj„)~2/3,-.^.  )  The  inelastic  hfetime  itself  is  temperature-dependent:        oc  T^, 
with  the  lifetime  decreasing  with  increasing  temperatures.  This  temperature  dependence 
of  the  inelastic  hfetime  is  responsible  for  the  temperature  dependence  of  the  conductivity. 
Thus  Eqn.  (1.5)  the  conductivity  equation  (in  2D) 


(t{L)  =  a(Lo)  - 


ae^  ln(L/Lo) 


(1.4) 


a{L)  =cr(Lo)  + 


ape^  \n{T/To) 


(1.5) 
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acquires  a  nonanalytic  log  T  temperature  dependence.  This  result  was  almost  immediately 
confirmed  by  experiments  performed  on  thin  metallic  films  and  silicon  metal-oxide- 
field-effect- transistor  (Si-MOSFET).  For  relatively  clean  samples,  the  resistivity  showed 
the  expected  logarithmic  increase  with  decreasing  temperature  and  applied  electric 
field  [40,  41,  42]. 

In  an  alternate  theory,  Altshuler  et  al.  [3]  considered  interaction  effects  in  the  frame- 
work of  a  dirty  Fermi  liquid  in  the  regime  when  kpl  »  1  (where  kp  is  the  Fermi  wave 
vector,  I  is  the  electronic  mean-free  path)  and  Tr  <  1  (this  is  the  diffusive  regime,  with 
r  being  the  scattering  time)  and  obtained  similar  nonanalytic  temperature  dependence 
for  the  conductivity.  The  conductivity  obtained  is  expressed  in  terms  of  the  Fermi-hquid 
parameter  F 

e'{l-F)HTT) 
=  2^  '  ^^-^^ 

where  the  factor  F  is  small  for  weak  screening.  However,  Altshuler  et  al.  [3]  did  not 
consider  localization  effects  (phenomena  that  can  occur  without  electron-electron  interac- 
tion). That  is,  the  maximally  crossed  impurity  fines  in  the  current  vertex  that  give  rise  to 
localization  were  disregarded. 

Both  correlation  effects  and  localization  effects  lead  to  logarithmic  temperature 
dependence  of  conductivity.  The  information  about  these  effects  can  be  disentangled 
by  experiments  performed  in  the  presence  of  a  magnetic  field  [43],  since  the  correlation 
effects  and  localization  effects  show  significantly  different  behaviors  [44,  45].  With  good 
agreement  between  theory  and  experiment,  the  nature  of  the  conducting  state  of  a  2D 
system  was  considered  resolved. 

Yet  theoretical  results,  particularly  those  of  Finkelstein  [46,  47]  were  contrary  to  the 
general  viewpoint.  By  performing  field-theoretic  renormahzation-group  calculations,  he 
showed  that  the  resistance  of  the  system  scales  to  weak  coupling  leading  to  a  metafile 
state  at  zero  temperature  (his  theory  explicitly  suppressed  the  maximally  crossed  dia- 
grams that  give  rise  to  logarithms  in  the  locafization  problem).  Altshuler  and  Aronov  [48] 
then  showed  that  by  considering  spin-orbit  or  spin-flip  scattering,  the  Finkelstein  solution 


10 

to  the  resistance  scales  to  strong  coupling,  and  shows  an  insulating  behavior  [48,  49],  thus 
paving  the  way  for  metal-to-insulator  transition  to  be  observed  in  D  >  2  systems. 

With  the  advent  of  better  technology,  it  became  possible  to  fabricate  samples  with 
mobihties  an  order  of  magnitude  higher  than  the  ones  developed  earher.  The  earliest 
experiments  which  reported  the  metallic  behavior  of  resistivity  were  performed  on  high 
mobility  Si-MOSFET's  by  Kravchenko  et  al.   [50,  51].  Both  metallic  and  insulating 
behaviors  were  observed  in  the  same  Si-MOSFET  sample  by  just  varying  the  density  of 
electrons.  Since  then,  numerous  studies  have  observed  the  metallic  state  in  systems  other 
than  Si-MOSFET's:  p-  and  n-Si/SiGe  [52],  p-GaAs  [53],  n-GaAs  [54],  n-AlAs 
[55].  Nearly  all  of  these  studies  have  reported  the  observation  of  some  key  features. 
Metalhc  temperature  dependence  of  the  resistivity  was  observed  for  densities  above  a 
sample-dependent  critical  density  Uc  in  the  ballistic  regime  Tr  >  1,  showing  almost  an 
order  of  magnitude  variation  of  resistivity  with  temperature.  For  densities  below  Uc  an 
insulating  behavior  was  observed.  The  2D  system  shows  an  unusually  sharp  response 
to  the  application  of  a  parallel  magnetic  field.  Conductivity  is  strongly  suppressed  with 
increased  magnetic  field,  up  to  the  saturation  field,  whereupon  conductivity  saturates. 
Interestingly,  the  field  required  for  saturation  of  conductivity  depends  on  the  density 
of  the  2D  sample  and  is  approximately  equal  to  the  field  required  for  complete  spin- 
polarization  of  the  2D  sample  [56] .  Since  the  parallel  magnetic  field  only  couples  to  the 
spin  degree  of  freedom  and  not  with  the  orbital  degree  of  freedom,  the  extreme  sensitivity 
to  parallel  magnetic  fields  implies  that  exchange  and  spin-related  interaxitions  play  a 
crucial  role.  Despite  agreement  as  to  the  general  behavior  of  2D  samples  on  the  insulating 
and  metalhc  side,  there  are  questions  on  which  general  consensus  has  not  been  reached. 
First,  does  the  resistivity  on  the  metalhc  side  continue  to  decrease  {dp/dT  >  0)  down 
to  the  lowest  temperatures  {i.e.,  no  upturn  toward  insulating  behavior  at  the  lowest 
temperatures)?  Prus  et  al.  [57]  have  observed  that  the  resistance  of  some  curves  with 
metaUic  behavior  show  an  upturn  upon  further  decrease  in  temperature,  signifying  a 
greater  role  played  both  by  electron-electron  interaction  and  weak  localization;  and 
signaling  the  possibility  that  the  zero-temperature  state  of  a  2D  sample  is  indeed  an 
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insulator.  In  contrast,  Kravchenko  et  al.  [50,  51],  found  no  resistance  upturn  for  Si- 
MOSFETS  with  similar  mobilities.  The  lowest  temperatures  they  achieved  were  an 
order  of  magnitude  smaller  than  those  achieved  by  Prus  et  al.  [57] ,  indicating  a  zero- 
temperature  MIT.  The  second  question  which  eludes  general  consensus  is  whether  the 
crossover  firom  metal  to  insulator  is  accompanied  by  a  transition  from  a  paramagnetic  to  a 
ferromagnetic  phase  and  if  so,  what  is  the  order  of  this  phase  transition.  This  is  discussed 
in  Section  1.4. 

The  theoretical  approach  to  understand  the  MIT  has  been  dependent  on  underlying 
assumptions  about  the  ground  state  of  2D  systems.  Next  we  briefly  summarize  two 
approaches  based  on  the  Fermi-liquid  picture  and  one  not  based  on  Fermi-liquid  theory, 
which  we  call  the  non-Fermi  liquid  approach. 

The  non-Fermi  liquid  approach.  Chakravarty  et  al.  [58]  assumed  that  the  ground 
state  of  disordered  electrons  is  a  non-Fermi  liquid.  The  single-particle  Green's  function, 
when  analytically  continued  to  the  lower  half-plane,  instead  of  containing  a  quasi-particle 
pole  consists  of  a  branch  point.  The  spectral  function  used  for  the  calculation  was 

M^^k)<x— — ^—  (1.7) 

where  a  is  positive  and  ujc  is  a  high-frequency  scale  proportional  to  the  inverse  of  the 
non-interacting  density  of  states.  For  a  <  1/2  they  found  that  the  non- Fermi  hquid  state 
is  locahzed;  whereas  for  a  >  1/2,  they  found  that  the  disordered  system  is  a  perfect 
conductor,  with  a  stable  non-Fermi  liquid  fixed  point. 

Dobrosavljevic  et  al.  [59]  took  a  phenomenological  approach  and  assumed  that  at 
large  conductance  interactions  can  modify  P  and  make  it  positive,  leading  to  a  metallic 
state.  By  assuming  the  form  of  the  function  P  to  be 


at  large  conductance,  the  conductance  obtained  showed  a  logarithmic  decrease  with 
increase  in  temperature 


P{g)  =  {d-2)  +  A/g^  A>0 


(1.8) 


g  -  (in(ro/r)) 


(1.9) 
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They  argued  that  this  state  cannot  be  a  Fermi  Hquid  because  of  the  unusual  temperature 

dependence  of  conductance  and  because  the  system  will  turn  to  an  Anderson  insulator 

once  the  interactions  are  turned  off. 

Fermi-liquid  approaches.  Among  the  Fermi-hquid  approaches  to  describe  the 

resistivity  of  a  2D  system,  the  most  prominent  are  those  of  Altshuler  et  al.  [3] ,  Finkel- 

stein  [46,  47],  Castellani  et  al.  [49],  and  Zala  et  al.  [60].  We  have  already  considered 

Altshuler  and  colleagues'  solution  which  predicts  localizing  behavior  in  the  diffusive  regime 

with  kpl  >  1.  Next,  we  briefly  review  Finkelstein's  solution  in  the  diffusive  regime  and 

Zala  and  colleagues'  approach  in  the  ballistic  regime,  both  of  which  predict  metallic 

behavior. 

•  The  Diffusive  regime: 
The  range  of  electron-densities  for  which  the  system  is  diffusive  (Tr  <C  1)  and 
also  show  metallic  behavior  is  rather  limited.  For  densities  such  that  p  ~  h/e^  and 
temperatures  T  <^Tf  (where  Tp  is  the  degeneracy  temperature),  the  system  hes  in 
a  diffusive  regime.  The  temperature  dependence  of  the  resistivity  in  this  regime  can 
be  described  by  the  renormalization  group  (RG)  equations  developed  by  Finkelstein. 
The  RG  equation  relating  the  resistivity  and  temperature  is 


i,  +  l-(4n^-l)(l±^ln(l  +  72)-l) 


(1.10) 


where     is  the  valley  degeneracy,  ^  =  -  ln(rr/7i),  the  dimensionless  resistivity  p  is 
in  the  units  of  h/e^-K  and  72  is  the  interaction  parameter  [61].  The  coupling  constant 
also  gets  renormahzed  and  its  equation  is  given  by 

These  RG  equations  [Eqs.  (1.10),  (1.11)]  can  describe  non-monotonic  resistivity 
behavior  as  a  function  of  temperature.  Starting  from  high-temperature  regimes 
where  the  system  shows  insulating  behavior,  a  decrease  in  temperature  causes  the 
resistivity  and  the  interaction  parameter  to  increase  (the  interaction  parameter 
shows  a  monotonic  behavior)  until  a  critical  interaction  strength  7^  is  reached.  At 
this  interaction  strength  the  RHS  of  Eq.  (1.10)  changes  its  sign.  Further  decrease 
in  temperature  causes  the  resistivity  to  decrease  and  hence  show  metallic  behavior. 
This  result  was  not  taken  seriously,  since  for  riy  =  1,  72  diverges  at  low  energies. 
The  situation  however  is  completely  different  for  a  two-valley  system      =  2  (Si- 
MOSFET  along  the  (001)  plane)  where  the  interaction  strengths  required  are  within 
the  perturbative  regime. 
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•  The  Ballistic  regime:  Zala  et  al.  [60]  showed  that  insulating  behavior  in  the  diffusive 
regime  and  metallic  behavior  in  the  ballistic  regime  are  both  governed  by  the  same 
physics:  coherent  scattering  of  electrons  by  the  self-consistent  potential  created 
by  all  the  other  electrons  (Friedel  oscillations).  For  the  baUistic  hmit  they  were 
able  to  show  with  the  usual  quantum-mechanical  arguments  that  the  scattering  of 
electrons  by  impurities  and  by  Friedel  oscillations  created  due  to  impurities  can 
be  coherent,  provided  the  scattering  is  a  back-scattering  event.  As  a  result,  the 
amphtude  of  back-scattering  processes  is  enhanced  compared  to  the  case  with  no 
Friedel  oscillation.  The  correction  to  conductivity  is  hnear  in  the  ballistic  regime  and 


where  the  sign  of  the  correction  to  conductivity  crucially  depends  on  the  Fermi- 
liquid  interaction  parameter  Fq  in  the  triplet  channel  and  near  the  Stoner  instability 
the  corrections  to  conductivity  is  greatly  enhanced.  This  result,  essentially  the 
effect  of  a  single  impurity,  when  extended  to  include  the  case  of  multiple  impurities 
results  in  a  conductivity  correction  that  is  logarithmic  in  temperature  as  obtained 
by  Altshuler  et  al.  [3].  The  nonanalytic  temperature  dependence  of  the  conductivity 
observed  in  both  the  ballistic  and  diffusive  regimes  arise  due  to  the  contributions 
from  long  range  correlations  or  the  soft  modes  [14].  The  soft  modes  which  are 
present  in  the  dirty  system  have  their  analogues  in  a  clean  system  as  well  and  are 
responsible  for  nonanalyticities  in  the  thermodynamic  quantities  hke  specific  heat 
and  spin  susceptibility. 
Although  there  has  been  considerable  progress  in  understanding  the  transport 

properties  of  2D  systems  in  the  ballistic  hmit,  the  question  regarding  the  nature  of  the 

ground  state  of  2D  electrons  has  not  been  resolved.  This  is  because  the  conductivity  can 

show  various  temperature  dependences  depending  on  the  strength  of  disorder  and  also 

on  the  concentration  of  electrons;  this  can  range  from  a  purely  Fermi-liquid  behavior 

with  large  conductivity  to  a  strongly  localized  behavior  (hence  a  non-Fermi  liquid) 

with  exponential  temperature  dependence  of  conductivity  in  the  dirty  hmit  kpl  <S  1. 

For  impurity  and  electron  concentrations  intermediate  between  these  two  limits,  the 

temperature  dependence  of  conductivity  can  show  various  functional  forms  and  the 

ground  state  can  vary  from  bubbles  of  Wigner  crystal  [62],  to  stripe  phase  to  bubbles  of 

Fermi  liquid  [63].  Thus  based  solely  from  the  temperature  dependence  of  resistivity  an 

unambiguous  judgment  on  the  ground  state  of  the  system  cannot  be  made. 


is  given  by 


(1.12) 


By  fitting  the  experimentally  determined  conductivity  with  Zala  et  al.'s  [60]  the- 
ory, which  imphcitly  assumes  Fermi-liquid  ground  state,  Fermi-liquid  parameters  Fq 
and  m*{rs)  were  extracted  [64,  65].  These  were  then  compared  with  the  parameters  ob- 
tained independently  from  other  experimental  techniques,  the  most  prominent  being  the 
Shubnikov-de  Haas  oscillation  in  the  conductivity  [66,  67].  An  excellent  match  between 
the  Landau  parameters  extracted  by  these  two  independent  approaches  points  to  the 
Fermi-liquid  nature  of  the  underlying  state  in  the  ballistic  limit.  At  the  same  time  there 
has  been  a  interesting  observation  which  does  not  fit  well  with  the  predictions  of  the 
Fermi-liquid  theory.  The  theory  predicts  that  the  effective  mass  is  renormahzed  by  the 
density  and  hence  the  Fermi  energy.  A  similar  renormalization  should  also  be  achieved 
by  a  parallel  magnetic  field  since  the  field  alters  the  Fermi  energy  and  the  distribution  of 
electrons.  The  prediction  for  a  field  dependent  effective  mass  is  for  a  clean  system,  still 
one  expects  a  qualitatively  similar  dependence  to  hold  for  a  dirty  system  which  shows 
Fermi-liquid  behavior.  However  no  dependence  of  the  effective  mass  on  the  field  has  been 
observed  [68].  The  issues  related  with  the  extraction  of  the  Fermi-liquid  parameters  from 
the  Shubnikov-de  Haas  data,  the  anomalous  effective  mass  and  the  nature  of  the  ground 
state  of  2D  electrons  are  explored  in  the  next  Section. 

1.4   Quantifying  Fermi-hquid  Parameters 

An  analysis  of  quantum  magneto-oscillation  in  terms  of  the  Lifshitz-Kosevich  (LK) 
formula  yields  information  regarding  the  quasi-particle  effective  mass  m*,  the  Pauli  (spin) 
susceptibility  x*  and  the  Lande  g*  factor.  These  quasi-particle  parameters  are  related  to 
the  harmonics  of  the  Fermi-liquid  interaction  via  the  following  relation: 

=  |:  -  1.  (1-13) 

where      is  the  interaction  term  in  the  triplet  channel,  and 

Ff  =  2(^-1),  (1,4) 

Ff  is  the  first  harmonic  in  the  singlet  channel  and  mt  is  the  band  mass. 

This  analysis  is  based  on  the  assumption  that  the  LK  formula  originally  developed 
to  describe  the  magneto-oscillations  in  3D  Fermi  liquid,  can  be  extended  to  describe  the 
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oscillations  in  2D  systems  by  a  change  in  the  electronic  spectrum  and  in  terms  of  the  zero 
field  Fermi-liquid  parameters  of  the  2D  system.  Recently,  there  have  been  questions  raised 
about  the  validity  of  LK  formula  in  2D  [69,  70].  In  the  following  we  will  discuss  this  issue 
in  more  detail. 

1.4.1    AppUcabihty  of  the  Lifshitz-Kosevich  Formula  in  2D 

In  3D  systems,  it  was  shown  by  Luttinger  [71]  and  also  by  Bychkov  and  Gorkov  [72] 
that  the  presence  of  a  magnetic  field  does  not  affect  the  zero  field  Fermi-hquid  parameters 
to  leading  order  in  1/7V  ,  where  iV  »  1  is  the  number  of  occupied  Landau  levels.  The 
correction  to  the  magneto-oscillations  obtained  by  including  the  oscillations  in  the  self- 
energy  is  of  the  order  of  A''"^  which  is  smaller  by  a  factor  of  y/N  compared  to  the  leading 
oscillatory  part.  Thus  additional  oscillations  of  the  self-energy  can  be  ignored  in  3D. 
The  situation  in  2D  is  different:  the  contribution  coming  from  the  oscillations  in  the  self- 
energy  is  l/N'^  smaller  than  the  monotonic  part  and  is  of  the  same  order  as  the  leading 
oscillatory  term  in  the  magneto-oscillation.  Thus,  for  zero  temperature  the  LK  formula 
in  2D  is  not  apphcable  for  describing  the  oscillations  [69,  70,  73].  On  the  other  hand, 
experiments  performed  at  finite  temperatures  seem  to  agree  well  with  the  LK  formula, 
thus  the  range  of  vahdity  for  the  LK  formula  requires  a  careful  inspection.  The  amplitude 
of  the  k^^  harmonic  of  the  magneto-oscillation  is  given  by 


where  e„  =  n{2n  +  1)T  is  the  Matsubara  energy,  lOc  =  eB/mc  is  the  cyclotron  frequency 
and  S  is  the  self-energy  of  electrons.  The  self-energy  for  a  generic  Fermi  liquid  in  the 
presence  of  short-ranged  impurities  is  given  by  i'Eo{£n,  T)  =  ae„  +  sgn(e„)/2r.  One  obtains 
a  simplified  form  of  the  amplitude  By  expressing  the  effective  mass  as  m*  =  m(l  +  a),  one 
obtains  a  simplified  form  of  the  amplitude 


(1.15) 


(1.16) 
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where  To  =  l/27rT  is  the  Dingle  temperature.  For  temperatures  and  impurity  concentra- 
tions such  that 


2t,\TIuITdIuj,)  >  1, 


(1.17) 


the  amplitudes  of  the  oscillatory  components  will  be  exponentially  suppressed,  the 
suppression  being  stronger  for  higher  amplitude  terms  [70] .  A  similar  argument  holds 
for  the  oscillatory  part  of  the  self-energy,  and  only  the  first  harmonic  needs  to  be  re- 
tained. The  contribution  from  the  first  harmonic  of  the  self-energy  oscillation  to  the 
magneto-oscillations  is  of  the  same  order  as  the  second  harmonic  of  the  quantum  magneto- 
oscillation.  This  imphes  that  the  first  harmonic  of  quantum  magneto-oscillation  does  not 
obtain  any  additional  contributions  from  the  oscillations  in  the  self-energy,  hence,  as  long 
as  only  the  first  harmonic  is  the  major  contributor  to  the  quantum  magneto-oscillation 
one  can  still  use  the  LK  formula  for  analyzing  the  data. 

The  above  analysis  was  performed  for  the  de  Haas- van  Alphen  effect,  however,  it 
is  conjectured  that  many  of  the  arguments,  especially,  the  vaUdity  of  the  first  harmonic 
term  also  holds  for  the  description  of  Shubnikov-de  Haas  oscillations  in  the  conductivity 
for  large  temperatures.  Thus  the  quasi-particle  parameters  which  get  renormalized  due 
to  electron-electron  interactions  (within  the  Landau's  Fermi-liquid  theory)  such  as  the 
spin  susceptibihty  x*,  effective  mass  m*  and  Lande  factor  g*,  can  be  extracted  from  the 
Shubnikov-de  Haas  measurements.  In  this  regard  it  is  worth  mentioning  that  in  a  recent 
experiment,  an  opposite  phase  of  the  second  harmonic  term  as  compared  to  the  one 
predicted  by  LK  formula  was  observed  [74].  This  can  be  attributed  to  the  inapphcability 
of  the  LK  formula  beyond  the  first  harmonic. 
1.4.2    Absence  of  Polarization  Dependent  Effective  Mass 

The  amphtude  of  the  oscillations  encodes  the  effective  mass  and  the  Dingle  tempera- 
ture, both  of  which  are  related  to  the  amphtude  via  the  following  relation 


the  slope  of  ln(fe/p)  vs  T  yields  the  effective  mass,  whereas  the  intercept  gives  the  Dingle 
temperature.  The  amplitude  and  the  oscillatory  part  are  multiplied  by  a  third  term  which 


(1.18) 
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is  also  called  the  Zeeman  term 


Zfc  =  cos 


cos[^],  (1.19) 


eB± 

where  A*  is  the  Zeeman  energy  and  fkOc  is  the  cyclotron  energy.  The  difference  (n|  -nj  oc 
X*Btot,  is  directly  proportional  to  the  total  applied  field  Btot  and  for  fields  with  no  parallel 
component  Zk  is  a  constant.  The  presence  of  non-zero  parallel  component  causes  Zk  to 
be  field  dependent  leading  to  beats.  Prom  the  position  of  the  nodes  in  the  beat  pattern 
X*  is  extracted  and  the  ratio  x*/"^*  gives  g*.  The  nodes  occur  at  regular  half-integer 
values  of  A;A*/cj*  [66].  Perfect  nodes  implies  that  the  oscillatory  terms  due  to  the  up- 
and  the  down-spin  electrons  have  the  same  amplitude  and  hence  the  same  effective  mass, 
a  result  in  contradiction  with  Fermi-liquid  theory.  Indeed,  the  presence  of  a  magnetic 
field  causes  the  Permi  surface  to  spht  into  two,  one  for  the  up-spin  electrons  and  another 
smaller  one  for  the  down-spin  electrons.  A  consequence  of  this  phenomenon  will  be  that 
the  velocities  of  up-spin  electrons  will  be  higher  compared  to  those  of  down-spin  electrons. 
Thus  within  the  frame  work  of  Fermi-hquid  theory  it  is  natural  to  assume  that  there 
are  three  distinct  interaction  terms,  an  interaction  between  two  up-spin  electrons  (fast 
moving  electrons),  one  between  two  down-spin  electrons  (slow  moving  electrons  and  hence 
strong  interaction)  and  the  interaction  between  a  down-spin  and  an  up-spin  electron 
(Section  4.1.1).  These  interaction  terms  in  turn  imply  different  effective  masses  for  the 
up-  and  down-spin  electrons.  This  intuitive  expectation  is  confirmed  by  a  perturbative 
calculation  (Section  4.1.2).  These  results  are  in  contradiction  with  the  experimental 
findings  on  Si-MOSFETs  with  a  (001)  surface,  which  report  absence  of  field  dependence 
of  the  effective  mass  and  also  indicate  that  the  effective  masses  of  up-  and  down-spin 
electrons  are  the  same,      =  [68]. 

To  resolve  the  paradox,  we  have  turned  to  a  large  A''-model.  We  have  found  the  spin- 
valley  degeneracy  A''  plays  an  important  role  in  the  mass  renormalization.  Each  of  the 
components  shows  polaronic  behavior,  with  the  mass  renormalized  primarily  by  the  high 
energy  plasmons.  For  typical     -  2  -  6  (where     =  V^me^/kp)  the  percentage  change 
in  the  effective  mass  between  unpolarized  and  a  fully  polarized  state  is  within  1  -  3%,  a 
result  which  is  consistent  with  the  experiment. 
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1.4.3   Effective  mass  and  Lande-g  Factor  Near  the  Critical  Densities 

A  number  of  groups  have  come  up  with  interesting  conclusions  regarding  the  ultimate 
fate  of  X*  and  m*  as  the  densities  approach  the  critical  one.  Shaskin  et  al.  [75]  and 
Vitkalov  et  al.  [76]  have  both  reported  divergence  of  x*  in  Si-MOSFETS  near  the  critical 
regimes,  with  the  former  arguing  that  the  divergence  in  x*  is  due  to  the  divergence 
in  m*  whereas  the  5*-factor  shows  only  a  nominal  increase.  Vitkalov  et  al.  however 
report  the  divergence  of  g*  as  the  cause  for  the  divergence  in  x*-  These  conflicting 
results  point  to  different  phases  of  the  ground  state  on  approach  to  the  critical  density: 
a  divergence  in  m*  points  toward  Wigner  crystaUization,  whereas  the  divergence  in  g* 
points  toward  ferromagnetism.  The  mutually  opposing  results  of  the  above  experimental 
groups  are  in  further  disagreement  with  the  results  of  a  third  group.  Pudalov  et  al.  [66] 
obtained  an  increase  in  both  g*  and  m*  near  the  critical  density,  however  the  increase 
in  m*  was  by  a  factor  of  3  and  only  a  small  increase  in  the  value  of  g*  was  seen.  Clearly 
the  experimental  evidence  for  the  very  existence  of  phase  transition  and  its  order  is 
contradictory,  nevertheless  it  remains  an  interesting  point  that  if  the  phase  transition  does 
exist  in  2D,  what  is  the  order  of  this  phase  transition. 

The  scenario  of  mass  divergence  near  the  MIT  regime  is  similar  to  the  one  proposed 
by  Brinkman  and  Rice  [77],  who  showed  using  Gutzwiller's  variational  method  [78]  that 
at  half-filling  the  effective  mass  diverges  at  the  critical  value  of  the  interaction.  This 
divergence  is  accompanied  by  the  vanishing  of  the  discontinuity  in  the  single-particle 
occupation  number  thus  indicating  a  deviation  from  the  Fermi-hquid  ground  state. 
Galitski  et  al.  [79]  and  Chubukov  et  al.  [80]  have  proposed  that  at  the  critical  interaction, 
the  quantum  critical  behavior  is  associated  with  a  density-wave  instability  at  a  finite 
wave  vector  Qc.  At  the  critical  point,  the  interaction  diverges  for  finite  wave  vector  Qc 
leading  to  the  divergence  of  the  fermion's  self-energy  and  hence  the  effective  mass.  It  was 
further  argued  by  Chubukov  et  al.  [80]  that  in  the  near  vicinity  of  the  critical  point  the 
self-energy  has  only  frequency  dependence,  and  the  divergence  of  the  effective  mass  right 
at  the  critical  point  and  not  before  it  indicates  that  the  renormalization  factor  Z  remains 
non-zero,  thus  the  Fermi-liquid  state  survives  till  the  critical  point. 


CHAPTER  2 

INTERACTING  FERMIONS  IN  TWO  DIMENSIONS 
2.1  Introduction 

The  nature  of  the  low-energy  electronic  excitation  in  a  2D  system  has  been  of  great 
interest,  particularly  since  the  discovery  of  high  Tc  superconductors.  It  was  suggested 
by  Anderson  that  forward  scattering  of  two  electrons  leads  to  vanishing  renormalization 
factor  Z,  i.e.,  to  the  breakdown  of  the  Fermi- liquid  description  [8].  Subsequently,  the 
effect  of  forward  scattering  has  been  investigated  by  many  authors  within  the  framework 
of  the  f-matrix  approximation,  and  an    log  e  term  as  the  leading  correction  to  the 
imaginary  part  of  the  self-energy,  ImE,  was  obtained.  This  result  indicates  that  the  real 
part  of  the  self  energy  is  still  dominant  compared  to  the  imaginary  part  of  the  self-energy 
and  hence  the  quasi-particles  are  well  defined  [30,  31,  32]. 

In  a  recent  work,  it  was  shown  that  for  fermions  with  a  linearized  dispersion,  ImE 
acquires  a  log  singularity  on  the  mass  shell  even  at  the  second  order  in  the  interaction  [22, 
29,  31,  32].  In  this  thesis  we  show  that  the  divergences  on  the  mass  shell  at  higher  orders 
in  perturbation  theory  become  power-law  and  are  more  singular.  Although  the  divergences 
in  the  individual  terms  are  cut  off  due  to  curvature  effects  (non-linear  dispersion),  the 
series  itself  does  not  converge.  This  result  seems  to  support  the  Anderson's  conjecture 
regarding  the  breakdown  of  Fermi  liquid  in  2D.  However,  this  is  not  the  complete  picture. 

At  each  order  in  the  interaction,  the  singularity  arises  due  to  the  interaction  of 
quasi-particles  with  excitations  near  the  upper-edge  of  the  particle-hole  continuum.  By 
re-summing  the  singular  perturbative  terms  to  all  orders,  we  were  able  to  obtain  a  closed 
form  for  ImE.  This  form  shows  that  the  singularities  are  due  to  the  interaction  of  quasi- 
particles  with  the  collective  excitation  of  the  system,  the  zero-sound  mode  (ZS).  The 
re-summation  regularizes  the  power-law  divergences  whereas  remaining  log  divergences  are 
removed  by  introducing  finite  curvature.  As  a  result  Fermi-liquid  theory  is  restored.  The 

outcome  of  the  re-summation  process  is  nonperturbative  contributions  to  the  self-energy, 
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one  of  them  being  a  zero-sound  contribution  to  the  imaginary  part  of  the  self-energy 
ImEzs  that  is  non-monotonic  as  a  function  of  the  distance  from  the  mass  shell.  The 
non-monotonic  behavior  leads  to  a  non-Lorentzian  shape  of  the  spectral  function. 

The  contribution  from  the  nonperturbative  zero-sound  mode  to  the  real  part  of 
the  self-energy  ReEzs  is  strongly  enhanced  near  the  mass  shell  and  is  responsible  for  a 
nonanalytic  U'^e\s\  contribution.  This  result  suggests  that  the  nonperturbative  processes 
contribute  to  a  nonanalytic      term  in  the  specific  heat.  This  however  is  not  the  case.  We 
discuss  the  zero-sound  contribution  to  the  specific  heat  in  Chapter  3. 

This  Chapter  is  organized  as  follows.  We  begin  in  Section  2.2  by  introducing  the 
relevant  scattering  processes,  these  are  forward-  and  back-scattering  processes.  Section  2.3 
is  essentially  concerned  with  the  self-energy  of  fermions  interacting  with  a  contact 
potential.  We  show  in  Section  2.3.1  and  Section  2.3.2  that  forward  scattering  leads  to  the 
divergence  of  the  imaginary  part  of  the  self-energy  on  the  mass  shell  at  second  and  higher 
orders  in  the  interaction.  The  procedure  used  for  the  resummation  of  diagrams  at  all 
orders  is  discussed  in  Section  2.3.3.  This  is  followed  by  the  calculation  of  the  renormalized 
imaginary  part  of  the  self-energy  in  Section  2.3.4  and  the  real  part  of  the  self-energy 
in  Section  2.3.5.  The  results  from  Sees.  2.3.4  and  2.3.5  are  used  to  obtain  the  spectral 
function  in  Section  2.3.6.  The  spectral  function  obtained  shows  features  that  deviate  from 
the  usual  Lorentzian  function  for  a  quasi-particle.  The  non-Lorentzian  feature  originates 
due  to  the  interaction  between  a  quasi-particle  and  the  zero-sound  mode  which  results 
in  a  non-monotonic  behavior  of  ImE;  the  effects  on  ImS  are  transferred  to  the  spectral 
function.  The  nonperturbative  processes  for  fermions  interacting  via  a  Coulomb  potential 
also  lead  to  non-monotonic  behavior  of  ImE,  as  discussed  in  Section  2.3.7.  The  final  part 
of  this  Section,  Section  2.3.8,  discusses  nonanalytic  corrections  to  the  tunneling  density 
of  states.  We  show  that  the  tunneling  density  of  states  does  not  acquire  a  nonanalytic 
correction  to  order  U^.  We  end  the  Chapter  by  concluding  in  Section  2.4. 

2.2    Scattering  Processes 

Scattering  events  between  two  low-energy  quasi-particles  of  momenta  ki  and  k2 
can  be  divided  into  three  categories.  (1)  Exchange  of  large  momentum  of  the  order 
Sk  =  \ki  -  k2\  oc  kp,  (2)  exchange  of  small  momentum  ^A;  ~  \Q\/vf  <  kp  and,  (3) 
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Figure  2-1:  Scattering  processes  responsible  for  divergent  or  nonanalytic  self-energy  con- 
tributions in  2D.  A)  Parallel- moving  quasi-particles  scatter  at  each  order  by  exchange  of 
a  small  momentum.  B)  The  quasi-particles  are  moving  opposite  to  each  other,  initial  and 
final  momenta  remain  essentially  unchanged.  C)  Each  of  the  oppositely  moving  quasi- 
particles  undergo  momentum  change  close  to  2/cf. 

exchange  of  momentum  near  2/cf  such  that  \\Sk\  —  2/c/r|  ~  IQI/uf  ^  kp.  Scattering 
events  of  the  first  kind  axe  responsible  for  analytic  contributions  to  the  self-energy 
and  thermodynamics(Section  2.3).  Scattering  events  of  the  second  and  third  kinds  are 
responsible  for  nonanalytic  contributions  to  the  self-energy  and  thermodynamics.  The 
effect  of  the  scattering  events  involving  small  momentum  exchange  or  large  2kF  scattering 
becomes  more  pronounced  with  the  lowering  of  the  spatial  dimensionalty,  since  the 
relative  phase-space  for  these  events  increases  in  lower  dimensions  compared  to  the  higher 
dimensions. 

Prom  an  earlier  analysis  of  small-momentum  scattering  events  [22],  it  was  established 
in  2D  that  the  kinematics  responsible  for  nonanalj^icity  in  the  self-energy  are  ID  in 
nature.  These  ID  processes  can  be  subdivided  into  two  kinds:  A)  those  in  which  the  two 
quasi-particles  move  almost  in  the  same  direction  and  undergo  small  momentum  transfer 
so  that  all  four  momenta  are  close  to  each  other  /cj  ~         k[  ~      this  process  is  also 
called  the  ^4  process  and  its  diagrammatic  representation  is  shown  in  Fig.  2-lA,  and 
B)  g2  processes  which  involves  quasi-particles  with  anti-parallel  momentum  undergoing 
small  momentum  transfer  -^1  ~  A;2  ~  -^j  ~  fc^.  The  diagrammatic  representation  of 
a  52  process  is  shown  in  Fig.  2-1 B.  The  third  kind  of  scattering  event  responsible  for  the 
nonanalytic  term  in  the  self-energy  involves  exchange  of  2kF  momentum  and  is  essentially 
ID  in  character  due  to  the  quasi-particle  occupancy  restriction.  The  2fcf  process  is 
depicted  in  Fig.  2-lC. 
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2.3    Self-Energy:  Nonanalyticity  and  Singularity 
In  this  Section,  we  will  first  derive  the  imaginary  part  of  the  self-energy  to  second 
order.  We  obtain  two  interesting  results:  a  nonanalytic  term  in  the  self-energy  and  a 
mass-shell  singularity.  We  will  consider  both  of  these  issues  in  detail.  The  Hamiltonian  of 
our  system  is  the  same  as  for  the  Coulomb  potential  case,  except  that  the  potential  is  now 
short  ranged. 

Before  we  go  into  the  details  of  the  calculation,  we  will  discuss  here  through  some 
simple  arguments  the  origin  of  the  nonanalytic  term  in  the  self-energy.  An  analytic 
self-energy  forms  a  regular  series  in  e,  the  real  part  of  the  self-energy  being  expressed 
as  ReE^(e)  =  aie  +  a2e^  +  a^e^  +  •  •  • ,  whereas  the  imaginary  part  of  the  self-energy 
is  expressed  through  the  series  ImE^(e)  =  bis'^  +  b2e'^  +  635^  +  •  •  • .  An  example  of  a 
nonanalytic  term  in  ImE^,  which  we  will  often  encounter  is  log  |£-|.  Let  us  consider 
ImE^  at  r  =  0, 

ImS^(e)  ^      dioj  d^qlmG^e  -uj,k-  q)lmV^{uj,  q),  (2.1) 

where  ImG'*  and  ImF"  are  the  imaginary  part  of  the  Green's  function  and  the  potential 
respectively.  Prom  this  form,  it  is  clear  that,  if  upon  a  momentum  integration  we  obtain  a 
uj  term,  then  one  ends  up  with  a  regular     term  after  the  remaining  frequency  integration 
is  performed,  whereas  if  the  momentum  integration  yields  a;  log  |a;|  (as  in  2D)  or      (as  in 
3D)  then  the  frequency  integration  yields  a  nonanalytic  e^log  |e|  or       term  respectively. 

Now  to  be  a  bit  more  exphcit,  the  potential  term  V{r,  t)  is  real,  thus  lmV{q,  u)  is  odd 
in  frequency  and  we  can  express  it  in  terms  of  a  characteristic  momentum  Q  of  the  system 
(which  serves  as  an  upper  cut-off)  in  the  following  form 

lmV{q,uj)  =  uF{^,^),  (2.2) 

where  F  is  now  a  even  function  of  frequency.  The  angular  integration  over  ImG^  yields  on 
the  mass  shell  (e  =  ^t;  where  ik  =    jlm  -  k\l2m) 

[ deimG''  =  -  /  den6{u  -  Vf^.q  +  q^ jlm)  =  ^AuC^^^^^Y  (2.3) 

J  J  VfQ  VfQ 
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where  the  subscript  D  stands  for  the  dimension, 

A^{x)   =  9{l-\x\) 

The  function  Ad  primarily  serves  to  impose  a  lower  cutoff  q  >  \uj\/vf  and  we  can 
ignore  the  specific  functional  form.  We  plug  the  result  of  Eq.  (2.2)  and  Eq.  (2.3)  back  in 
Eq.  (2.1)  to  obtain 

ImS«(e)  ~  r  du;u;  [         dqq^-'F{^,  J^).  (2.4) 

Now  if  the  momentum  integral  is  dominated  by  large  momenta  of  the  order  Q,  then  the 
function  F  to  leading  order  can  be  considered  to  be  independent  of  frequency  (since 
\u\/vfQ  ~  £/vfQ  <S  1,  as  Q  the  upper  momentum  cut-off  is  large  and  e  is  small)  and 
one  can  replace  the  lower  limit  by  0  as  well.  The  momentum  and  frequency  integrations 
decouple,  and  one  obtains  an  analytic  £^  term.  Thus  a  process  which  involves  large 
momentum  transfers  yields  analytic  terms.  Let  us  consider  a  case  where  the  decouphng 
between  the  momentum  and  frequency  is  not  possible.  Such  a  scenario  takes  place  at  the 
second  order  in  the  interaction  for  a  contact  potential.  The  retarded  part  of  the  effective 
potential  is  a  constant  times  the  polarization  bubble, 

ImV^iu,  q)  =  -UHiaU'^iu;,  q),  (2.5) 

where  11^  has  the  following  general  structure  in  an  arbitrary  dimensionality 

Imn«(a;,  q)  =  -ud—Bd{~),  (2.6) 
vpq  vpq 

where  Vd  is  the  density  of  state  in  D  dimensions  and  the  dimensionless  function  Bd 
imposes  a  constraint  uj  <  vpq.  For  D  =  3,  we  obtain  from  Eq.  (2.1) 

1  "  " 


2    /     J,  ,    /  j„„2 


ImE(e)  ~       [  duj  [  dqq 

Jo  Jq>\t^\/VF 


qkf 


UJ 


Vpq 


(2.7) 


where  we  have  neglected  the  constants,  and  the  functions  Bz  and  A3.  The  first  square 
bracket  represents  the  result  from  the  angular  integration,  whereas  the  second  one 
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a) 


a2) 


V  bs,V 


b5)  '^^^^^b6)       ^  b7)^T-^5V^ 

Figure  2-2:  Self-energy  diagrams  witii  explicit  and  implicit  polarization  bubbles. 


represents  Imll^.  The  upper  cut-off  is  of  the  order  of  kp-  Performing  the  momentum 
integration  we  obtain 


ImS(£)  f  du)u 

Jo 

~  ae^  -  b\ef, 


FL 


Vf 


Beyond  FL 


(2.8) 


where  the  first  term  originates  from  the  large  momentum  transfer  regime  and  is  the  Fermi- 
Uquid  result,  the  sub-leading  second  term  originates  from  the  small-momentum-transfer 
regime  and  is  nonanalytic.  For  reduced  dimensions,  the  effective  interaction  plays  a  more 
singular  role  with  the  result  that  the  nonanalytic  terms  become  the  leading  terms  in  the 
imaginary  part  of  the  self-energy.  In  2D,  the  nonanalytic  term  behaves  as     log  |e|,  a  more 
dominant  term  than  the  Fermi-liquid     term.  In  the  Sections  to  follow,  we  will  consider 
the  origin  of  nonanalyticity  in  greater  detail. 
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We  will  define  the  Fermion  self-energy  (E)  via  the  Dyson  equation  in  terms  of  the 
exact  (G)  and  bare  (Go)  Green's  function 


1-1 


Go^  +  E,  (2.9) 


note  that  we  have  defined  S  with  an  opposite  sign  compared  to  Refs.  [16,  17]. 
2.3.1    Self-Energy  at  the  Second  Order 

At  second  order  in  the  interaction,  there  are  only  two  non-trivial  diagrams  (Fig.  2- 
2),  the  rest  of  the  irreducible  diagrams  in  effect  renormalize  the  chemical  potential  and 
hence  are  not  important.  For  a  contact  interaction,  V{r)  =  U5{r),  the  two  second  order 
diagrams  in  Fig.  2-2  are  related,  namely,  diagram  (al)  is  —2  times  diagram  (a2).  This  is 
particularly  easy  to  see  in  momentum  space  as  diagram  (al)  =  211"^  Jg  Yl{Q)G{P  +  Q)  and 
diagram  (a2)  =  -U^     J^^  G{P  +  Q)G{P  +  Q  +  Qi)G{P  +  Qi)  =        Jq  n(Q)G(P  +  Q), 
where  Jg  is  the  integration  over  frequency  and  momentum  and  TL{Q)  is  the  polarization 
bubble  without  including  the  spin  summation.  In  the  retarded  formalism,  the  total 
contribution  at  second  order  in  the  interaction  is  given  by 


ImG5'(ei,pi)Imm(n,fc-pi) 
1" 


tt2         p  poo  POO 

E§{e,k)   =   j—^     d'pj  de, 

(27r)%y  n  +  £i-£-id 

X  tanh  (-^^  +  coth  ( | 
\2TJ  \2TJ 


(2.10) 


and  for  T  =  0,  the  imaginary  part  of  the  self-energy  acquires  the  form 

Tj2        /-O  /• 

ImSf  (e,  k)  =  -j^-^  /   dn  /  cPqlmG^ie  +  Q,    +  q)lmU^(n,  q),  (2.11) 
[InyTt  J 

where  ImG^(e  +  Q,k  +  q)  is 

ImG^(£  +  n,k  +  q)^  -7r5{Q  +  e  -  ^fc+g),  (2.12) 

with  =  k'^/2m  -  k%/2m.  We  will  be  interested  in  the  contributions  coming  from 
q  ~  Q,Vf  <C  kp  and  q  close  to  2kF.  The  form  of  the  polarization  bubble  for  small  q  is 
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whereas  for  q  f«  2kF  we  have 


(q-2kF       Uq-2kF\       fQ  +  iO+^^ 
^      y   2kF  2kF   J      V  2A;fVf  ,  ^  ^ 


(2.14) 


The  above  two  results  for  11^  may  be  obtained  by  using  the  result  of  Stern  [81];  alterna- 
tively, one  can  use  the  result  for  the  polarization  bubble  in  the  Matsubara  representation 


n"(ia.,,)  =  ~(^i-^^^^^^^j,  (2.15) 

where  x  =  (qkp)^  -  9^/4  +  m^Q^  and  y  =  q^mVtn,  expand  for  small  q  <^kF  and  q  near  2kF, 
and  perform  an  analytic  continuation  to  real  frequencies. 

The  nonanalyticities  in  Eq.  (2.13)  and  Eq.  (2.14)  represent  Landau  damping  and 
the  Kohn  anomaly,  respectively.  For  Q  <  vpq  damping  of  excitations  with  frequency  Q 
and  momentum  q  can  take  place  only  inside  the  particle-hole  (PH)  continuum,  similar 
to  the  damping  of  a  collective  mode  inside  the  PH  continuum.  The  polarization  bubble 
near  2fcf  shows  different  nonanalyticities  depending  on  whether  Q  <  2kF  or  Q  >  2kF 
(where  Q  represents  the  transferred  momentum).  For  Q  <  2/cf,  the  nonanalytic  part 

<x  iQ/{2kF  -  Qfl'^  and  for  Q  >  2kF  the  nonanalytic  part  11^  oc  (Q  -  2kFY^'^-  In  what 
follows,  we  will  calculate  the  forward  and  back-scattering  contribution  at  2"'^  order. 

Forward-  and  back-scattering  contribution  at  second  order.  In  an  earlier  work 
[15,  22],  the  leading  nonanalytic  contribution  to  the  self-energy  were  identified  to  originate 
from  processes  with  small  and  large  2kF  momentum  transfers.  To  logarithmic  accuracy, 
contributions  from  the  Kohn  (2A;f )  anomaly  to  the  nonanalytic  term  of  the  self-energy 
were  found  to  be  identical  to  the  p2  contributons.  In  the  following,  we  will  show  that  the 
nonanalyticity  which  originates  from  q     2kF  due  to  the  dynamic  nonanalytic  term  in 
the  particle-hole  response  function  is  responsible  for  the  nonanalyticity  in  the  ImE^.  The 
imaginary  part  of  the  polarization  bubble  near  Q  =  2kF+q  <  2kF,  is  given  by  the  dynamic 


term  mQ/4nvF^/\q\kF,  thus  the  contribution  to  the  imaginary  part  of  the  self-energy  from 
the  Q  <  2kF  region  is 

ImSf(e,fc)   =  Tdfi  (2.16) 

\^^)  Jo       J  '-4TrvF\/\q\kF^ 
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where  for  k  =  kp,  6+Q  =  i'^kj-  +  ^fQ  +  g^)/2m  +  {2kF  +  q)kFC0s9/m.  Performing  the 
angular  integration 

2e(H-|a|). 


J  de6{a 


bcos6) 


/ 


QdQ   =   {2kF)  /  dq 


and  keeping  the  most  dominant  terms  (with  typical  q  such  that  \q\vF  3>  Q),  we  obtain 


ImS2«^(£,fc)   =    j^J^dnJ  dq{2kF) 


y/Ak%\q\/rn^^  U7rvF\/i#F 


u2  ^ 


(2.17) 


where  Um/2i:  =  u  and  —W  represents  the  lower  cutoff  on  the  momentum  integration, 
which  is  typically  Ef-  The  other  nonanalj^ic  piece  of  Imll^  cc  [Q  —  2A;f)^/^  —  ^  which 
originates  from  Q  >  2kF  does  not  contribute  to  the  nonanalytic  piece  in  the  self-energy. 
The  contribution  to  the  nonanalytic  term  in  the  second  order  ImEf  [Eq.  (2.11)]  from  the 
q  =  Q  region  (we  need  to  keep  the  nonanalytic  piece  of  the  polarization  bubble  near  q  — 


i.e.  —TnQ/2Try/{vFqy  —  0.'^  )  is  expressed  as 

lmE^{6,k)   =  -^J\nJd\6ie-n-E,+,)[^-^ 


(2.18) 


{vFqy  - 

where  ^k+q  =  k^/2m  +  q^ /2m  +  kqcos6/m  —  A;^/2m.  As  before,  we  perform  the  angular 
integration  on  the  delta  function, 

2tt 


f  ddlmG^  =  — 

J 


^{vfQY  -  (C!  +  A)2 
using  this  result,  we  find  that  Eq.  (2.18)  acquires  the  following  form, 


ImE^(e,A;)  = 


C/2 


(27r)' 


qdq 


^(^)2_(A-n-e,)2 


(2.19) 


I27r ^ {vpqy  -n^^' 

where  =  k'^ /2m  -  k]p/2m  and  A  =  e  -  ^fc  has  the  interpretation  of  the  distance  from  the 
mass  shell. 
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For  A  =  0  and  neglecting  the  curvature  (e,  =  0)  we  obtain 


{^T^Y  Jo  JmivF 


mivF 


mQ 


1 


(2.20) 


i.e.,  an  unexpected  log  singularity  on  the  mass  shell.  The  origin  of  this  divergence  is  due 
to  the  form  of  the  polarization  bubble  at  small  momenta  which  is  square  root  divergent 
at  \Q\  =  VpQ  and  due  to  the  divergent  form  of  the  Green's  function  obtained  after 
performing  the  angular  integration.  Combining  the  two  contributions  and  performing  the 
d^q  integration  one  obtains  a  logarithmic  singularity  due  to  the  soft  modes  on  the  mass 
shell  (A  =  0).  Higher  orders  in  the  interaction  have  larger  number  of  particle-hole  bubbles 
and  hence  stronger  singularities.  These  case  will  be  considered  in  the  next  Section. 
By  keeping  A  finite,  we  obtain  the  following  form 


ImEf  (e,  k) 


dQQ  log 


El 


{A-kYvI 


A^  log 


L|A(2fi- A)| 
A  -2e 


A 


El 


mU 


+  2e(e  +  A)  +  4.^1og|^ll^_2^^|j 


167r3  ^ 
and  near  the  mass  shell 

ImEf(e,A;)  = 


"  r/  2    ^  ^  ,  /  2     ^  ^  1 


(2.21) 


u 


(2.22) 


It  has  been  estabhshed  that  the  two  log  contributions  in  Eq.  (2.22)  have  different  ori- 
gins [22] .  The  first  term  is  due  to  a  process  with  two  incoming  particles  of  almost  identical 
momentum  undergoing  small  angle  {q  ^  0)  scattering  (in  ID  terminology  it  is  the  g4 
process),  whereas  the  second  term  is  due  to  a  process  in  which  two  incoming  particles  of 
almost  opposite  momentum  undergo  small  momentum  scattering  (a  g2  process).  The  fol- 
lowing analysis  will  convince  us  that  there  are  only  two,  q     0  processes  which  contribute 
to  the     log  e  term.  We  will  reexamine  the  contribution  to  ImE^ 

Q 

in/vf  a/ {kpqy  -  ' 


ImE^   oc    f  dn  f     qdqde6{Q  +  A  -  kq/mcose) 

Jo  Jii/VF 
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where  k  is  the  momentum  of  the  incoming  particle.  The  imaginary  part  of  the  Green's 
function  yields  &  l/q  term  and  for  kpq  2>    the  imaginary  part  of  the  polarization  bubble 
yields  al/q  term,  the        term  obtained  is  responsible  for  the  nonanalytic  term.  The 
above  expression  reduces  to 


r dnn  f    ^  f  d66(^{n  +  A)-cose), 

Jo         Jn/VF  Q  J        ^kq  > 


for  small  r?  and  A  typical  angles  near  ±7r/2  will  yield  the  nonanalytic  e^loge  term,  this  is 
the  angle  between  the  incoming  momentum  k  and  the  transferred  momentum  q.  A  similar 
relationship  exists  between  the  momentum  transfer  and  the  incoming  momentum  p  of  the 
"other"  electron.  The  information  about  the  "other"  electron  is  in  the  imaginary  part  of 
the  polarization  bubble  and  is  given  by 

Imn^(fJ, 9)  (X  j    pdpd(f)^5{Q  —  £q— pq cos 4))  —  5{Q.  + £q— pq cos 4))^ 
=   -  f  dpddlsi— — —  -  cos(f>)  -  ^9  _  cose!)) I 

V¥W^'  ^^'^^^ 

thus  for  pq  >  Q,  we  find  ImlT^  ocQ/q  (the  l/q  form  of  the  polarization  bubble  is  crucial 
for  the  nonanalytic  term),  clearly  the  typical  angles  involved  in  this  contribution  are 
(f)     ±n/2.  Thus  both  the  fermionic  momenta  p  and  k  are  perpendicular  to  q,  implying 
two  processes,     and  g2,  are  responsible  for  the  nonanalytic  contribution  to  the  self-energy 
from  q  ^  0  regimes.  In  this  work,  we  will  classify  52  as  a  back-scattering  contribution  since 
it  involves  scattering  of  two  incoming  momenta  with  total  momentum  close  to  zero.  For 
a  contact  interaction,  the  total  back-scattering  contribution  which  is  a  sum  of  52  and  2kF 
processes  is  given  by 

Since  there  is  no  mass-shell  singularity  from  the  back-scattering  contribution,  we  do  not 
need  to  consider  higher  order  terms  and  hence  we  can  stop  the  series  at  the  second  order 
in  the  interaction  U.  We  will,  however,  pursue  higher  order  contribution  from  forward 
scattering  process  because  of  the  mass-shell  singularity  in  this  channel. 
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2.3.2    Higher-Order  Forward-Scattering  Contributions 

In  an  earlier  analysis,  we  have  established  that  the  singularities  from  the  soft  mode 
{q  ~  0)  in  the  particle-hole  bubble  causes  a  log  singularity  of  the  ImE^  at  the  second 
order.  At  higher  order  in  the  interactions,  these  singularities  proliferate  leading  to  a 
power-law  divergence  of  ImE^.  The  third  order  diagram  of  the  self-energy  consists  of 
a  maximum  of  two  polarization  bubbles,  having  the  same  momentum  and  frequency 
(Fig.  2-2), 

Q  l2 


(n«)2  =  i^? 


1  +  i 
1  +  i 


2Q 


Thus  at  the  third  order,  the  most  singular  term  in  Im(n^)^  is  given  by 

^2 


Im[n«]2  =  -£^[7rQ|fi|<5(fi2_452)]_ 


The  singular  term  at  the  third  order  is  given  by 


\mL^p{e,k)   =  U' 


dn  r  qdq 
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_Uhrf_  r 
327r3t;2 
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qvf) 
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vW^^A)A 


C/3 


(2.25) 


where  in  the  last  step  of  the  above  equation  the  singular  l/\fK  term  has  been  extracted. 
The  additional  0  function  signifies  that  the  expression  is  zero  if  A  and  e  are  of  the 
opposite  sign.  At  higher  orders  we  obtain  •  ■  .  •  ; 


rrn^n/2+l       .  ^  , 


(2.26) 
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for  n  >  2.  Collecting  forward-scattering  contributions  to  all  orders  in  u,  we  obtain 

'     '  n=l 

where  C„  are  the  numerical  coefficients.  This  series  converges  only  for  (u^e/A)  < 
1,  outside  this  regime  the  series  will  not  converge.  Thus  there  is  a  clear  need  for  re- 
summation  of  the  perturbative  series.  This  is  the  topic  of  our  next  Section. 
2.3.3   Re-summation  of  Forward-Scattering  Contributions 

A  proper  re-summation  of  the  self-energy  series  involves  keeping  at  each  order  in  the 
interaction  the  most  singular  terms.  We  find  that  such  terms  correspond  to  diagrams  with 
the  maximum  number  of  particle- hole  bubbles  (with  small  frequency  and  momentum). 
A  convenient  way  to  re-sum  the  contributions  will  be  to  first  identify  the  singular  four- 
fermion  vertex  T  and  then  relate  the  fermion  vertex  to  the  self-energy  E  via  the  Dyson 
equation. 

Four-fermion  vertex.  The  diagrams  that  form  paxticle-hole  bubbles  are  of  two  kinds 
(for  a  short  range  interaction):  the  ones  which  have  explicit  bubbles  and  the  others  with 
implicit  bubbles  (those  which  are  not  easily  identifiable).  These  explicit  and  impUcit 
particle-hole  bubbles  also  appear  in  the  self-energy.  A  typical  example  is  the  2"'^  order 
self-energy  where  we  consider  two  bubbles,  one  of  which  has  an  explicit  bubble  Fig.  2- 
2(al)  and  the  other  an  implicit  bubble  Fig.  2-2(a2).  An  explicit  bubble  comes  with  an 
additional  factor  of  -2  from  the  spin  summation.  Due  to  these  differences,  the  procedure 
of  finding  an  overall  pre-factor  at  order  u  is  rather  involved,  as  it  requires  counting  the 
number  of  diagrams  at  the  same  order  in  a  conventional  diagrammatic  technique  operating 
with  a  non-symmetrized  vertex,  f .  Once  we  obtain  the  non-symmetrized  vertex,  we  anti- 
symmetrize  it,  and  use  the  result  in  the  Dyson  equation  to  obtain  the  self-energy.  In  the 
following,  we  derive  an  expression  for  the  forward-scattering  part  of  the  non-symmetrized 
vertex,  F,  summing  up  diagrams  with  the  maximum  number  of  polarization  bubbles  to  all 
orders  in  the  contact  interaction,  U.  The  diagrams  for  a  non-symmetrized  vertex  up  to 
third  order  are  presented  in  Fig.  2-3. 


32 
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Figure  2-3:  Second  and  third  order  Vertex  diagrams  with  maximum  number  of  particle- 
hole  bubbles.  Additional  diagrams,  obtained  from  those  in  the  second  column  by  a  permu- 
tation a  — >  /3, 7     e,  are  not  shown. 
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In  the  Matsubara  technique,  we  associate  a  factor  of  —  f/  with  each  of  the  interaction 
Hnes,  and  a  factor  —2  with  each  of  the  polarization  bubble.  There  is  also  an  extra  factor 
of  —1  for  exchange  processes  in  which  the  two  outgoing  legs  are  permuted  (the  last 
diagrams  of  second  and  third  order  in  Fig.  2-3).  We  present  here  a  general  recipe  for 
calculating  the  u*''^  {v  >  1)  order  vertex  diagram.  The  vertex  consists  of  two  parts  (Fig.  2- 
3).  The  first  part  comes  from  the  direct  interaction  and  contains  a  spin  factor  Sa-ySp^-  The 
second  part  is  due  to  the  exchange  interaction,  and  comes  with  a  spin  factor  daeSp-f.  At 
each  order,  there  is  only  one  exchange  diagram  whose  contribution  is  —{—U)''Il''~^6a£Si3^ 
(At  second  and  third  orders,  these  are  the  first  and  third  diagrams  in  the  third  column  of 
Fig.  2-3,  respectively). 

The  rest  of  the  diagrams  are  due  to  the  direct  interaction  and  contain  various  number 
of  bubbles.  At  order  u  (number  of  interaction  lines  =  u),  the  number  of  bubbles  in  the 
diagram  (both  explicit  and  implicit)  will  he  u  —  1,  thus  the  number  of  explicit  bubbles  or 
the  rings  (R)  can  vary  from  0  to    -  1.  We  consider  a  diagram  with  R  rings  and  i/  -  1 
interaction  lines,  for  this  diagram  R  +  1  interaction  lines  will  be  used  up  in  connecting 
the  rings  (we  are  considering  connected  diagrams)  and  also  in  connecting  the  rings  to  the 
two  external  solid  lines.  The  remaining  N  =  v  —  R  -  I  interaction  lines  can  be  arranged 
an5rwhere  either  at  the  two  ends  of  the  chain  of  bubbles  (or  in  other  words,  dressing  up 
the  two  interaction  lines  which  connect  the  solid  line)  or  inside  the  R  ring  bubbles.  Thus 
there  will  be  5  =     +  2  sites  (two  ends  and  inside  the  R  rings)  where  N  interaction  lines 
can  be  placed.  The  number  of  diagrams  with  R  ring  diagrams  is  equal  to  the  number  of 
ways  to  arrange  A'^  lines  among  S  sites: 


(5  +  A^-  1)! 


(2.28) 


(5-l)!A!       (i?  +  i)!(j._i?_i)r 


Consequently,  the  contribution  to  T  from  diagrams  with  R  bubbles  is 


{R  +  l)\{u  -  R-l)\ 


(2.29) 


The  total  contribution  from  all  bubble  diagrams  at  the  order  u  is  then 


(2.30) 


=  (i-2r  =  (-ir. 


where  we  have  used  an  identity 

Adding  up  direct  and  exchange  terms,  we  obtain  the  following  form  for  the  non- 
symmetrized  vertex  at  order  y  >  I: 

The  vertex  function  can  now  be  readily  summed  to  all  orders,  with  the  result 

oo 
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(2.31) 


(2.32) 


■ye 


1/=1 

u 


(2.33) 


'''i-([/n)2  ""'-'^  1  +  uu{qy 

where  by  including  the  additional  factor  -USaeSg^  we  are  able  to  extend  the  summation  of 
the  vertex  function  to    =  1.  Using  an  SU(2)  identity 


SaeSp,  =  (1/2)  {a°^^a%  +  6,^60,) , 
and  introducing  the  dimensionless  spin  and  charge  vertices 


(2.34) 


1  1 


21 +  UU{qy 
G  =1-_L_ 

"  2i-t/n(9)' 

we  obtain  for  the  non-symmetrized  vertex  the  following  form 


(2.35) 


^a0,ie{P\,P2,Pi- q,P2  +  q)   =  T{q) 

=  -U 


.  (2.36) 


Finally  we  anti-symmetrize  the  vertex  by  the  following  procedure 


^a0;-re{q)  =  ^aP;-yeiq)  "  ^a^Mq)- 


(2.37) 
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For  a  repulsive  potential  the  retarded  charge  vertex,  Qp,  has  pole  at  the  zeros  of  1  - 
UU^iq,  Q)  =  0.  The  pole  exists  even  for  an  infinitesimally  small  U,  since  for  \Q\  >  \vfQ\, 
is  real  and  can  be  arbitrarily  large  as  |Q|  — >•  {vp'Q)'^. 

The  collective  mode  which  corresponds  to  the  short-range  interaction  is  called  a  zero- 
sound  mode  whose  dispersion  is  |Q|  —  cQ  and  the  zero  sound  velocity  is 

where  the  expansion  is  valid  for  m  <  1.  The  charge  vertex  exhibits  singular  behavior  near 
the  zero-sound  mode  and  has  the  following  form 

1 


+ 


2(Q  -  cg)|j[l  -  UU{Q,  Q)]  n=cQ    2(Q  +  cQ)§^[l  -  UU{Q,  Q)] 


n=-cQ 


For  small  values  of  u  the  zero-sound  mode  cQ  is  just  above  the  upper  boundary  of  the 
particle-hole  continuum,  given  by  VpQ.  The  residue  of  the  charge  vertex  is  not  a  constant 
and  is  proportional  to  Q'^,  thus  the  quanta  of  the  zero-sound  mode  are  not  free.  The  spin 
vertex      on  the  contrary  has  a  pole  in  the  imaginary  axis  and  remains  over-damped. 

Dyson  equation.  The  self-energy  due  to  forward  scattering  is  related  to  the  vertex 
function  via  the  Dyson  equation  [17] 

^F,a0{p)  =  Sa0  f  UG{p-q)-  [     Ur,,,,0{p,p'  +  p"  -p;p'y)Gip')Gip")G{p'  +  p"  -p), 

(2.38) 

where 

/■••^^E  /  d'k/{2i:f...  (2.39) 

Although  the  Green's  functions  in  the  Dyson  equation  are  exact  ones,  we  can  safely  ignore 
the  self-energy  insertions  into  the  diagrams  which  diverge  near  the  mass  shell  as  these 
insertions  do  not  lead  to  any  additional  singularities  near  the  mass  shell.  As  we  keep 
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u={Um/2'rc)  small,  regular  corrections  are  thus  irrelevant,  and  we  can  safely  use  bare  G's 
instead  of  the  exact  ones  in  Eq.  (2.38).  Substituting  Eq.  (2.37)  into  the  Dyson  equation 
(2.38),  we  obtain  for  the  self-energy 

S„^(p)   =  6^0  fuG{p~q)-  [  UT^a;-y0G{p'Mp-p')+  [  UT,c;g,G{p")U{p-p") 

Jq  Jp'  Jp" 

Rearranging  the  result,  we  obtain  the  self-energy  in  the  following  form  =  ^a/jSp 

where 


J  a 


^  1    WW{q)        3  U^n^iq) 


U  +  U'li{q)  + 


G{p-q).  (2.41) 


2{l-UIi{q))      2{l  +  UYl{q))_ 
We  remind  the  reader  that  in  the  perturbation  theory  ImE^  diverges  upon  approaching 
the  mass  shell-  logarithmically  to  second  order  in  f/,  and  as  l/^e  -     to  third  order.  It  is 
convenient  to  rearrange  the  terms  of  Eq.  (2.41)  and  decompose  Sp  into  three  parts  making 
use  of  the  charge  and  spin  vertices,  introduced  in  Eq.  (2.35),  as 

Sf(p)  =  (p)  +  (p)  +  Sex;  (2.42a) 
Sp(p)   =   U  !G,{q)G{p-q)-  (2.42b) 

Jq 

lla  (p)   =   SU  [g„  (q)  Gip-q);  (2.42c) 

Jq 

=    f[2U-  U^Uiq)]  G{p-q).  (2.42d) 

Jq 

Terms  Hp  and  E^^  correspond  to  the  interaction  in  the  charge  and  spin  channels,  respec- 
tively, and  are  summed  to  all  orders  in  U.  The  remainder,  Eex,  contains  extra  contribu- 
tions of  the  first  and  second  orders  in  U,  which  are  not  included  in  the  spin  and  charge 
terms  and  are  required  to  reproduce  the  correct  expansion  of  the  first  and  second  order 
term  of  the  self-energy. 
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A  similar  expression  was  obtained  for  the  self-energy  in  the  paramagnon  model  (the 
spin-fluctuation  model)  [20],  although  the  contribution  from  the  remainder  term,  Eex, 
was  not  taken  into  account.  The  paramagnon  model  includes  only  two  kinds  of  diagrams, 
RPA  and  ladder  diagrams  corresponding  to  the  charge  channel  and  the  spin  channel, 
respectively.  In  contrast,  our  choice  of  diagrams  is  based  on  a  broader  consideration  of 
including  the  maximum  number  of  particle  hole  bubbles  at  each  order  which  for  forward 
scattering  leadto  mass-shell  singularity,  thus  we  have  included  non-RPA  and  non-ladder- 
tj^e  diagrams  as  well.  Including  the  maximally  divergent  diagrams  also  means  that 
the  perturbative  expansion  near  the  mass  shell  is  in  terms  of  the  diverging  parameter 
(([/m/27r)^e/A)  and  any  expansion  which  involves  the  regular  f/m/27r  is  neglected.  After 
the  complete  resummation,  it  turns  out  that  the  expression  for  the  self-energy  is  very 
similar  to  the  one  obtained  in  Ref.  [20]. 
2.3.4    Renormalized  Imaginary  Part  of  the  Self-Energy 

The  imaginary  part  of  the  self-energy  acquires  non-zero  contributions  from  the 
regions  in  momentum  and  frequency  space  where  the  effective  interaction  has  a  non-zero 
imaginary  part.  There  are  only  two  such  regions,  one  of  them  at  the  collective  mode 
\Q\  =  cQ  where  the  charge  channel      has  a  pole,  and  the  other  region  where  Imll^  ^  0, 
i.e.,  where  particle-hole  excitations  are  permissible.  All  three  terms-the  charge-vertex  part, 
the  spin- vertex  part  and  the  extra  term—  acquire  contributions  from  the  particle-hole 
region,  whereas  only  the  charge- vertex  part  has  a  contribution  from  the  collective  mode 
excitation.  We  can  represent  the  imaginary  part  of  the  self-energy  as  follows 

ImE^   =  ImE^  +  ImE^  +  ImSf, 

=  (ImE«  +  ImE^  +  ImS^)pH  +  (ImEj)zs,  (2.43) 

where  the  subscripts  PH  and  ZS  stand  for  the  particle-hole  and  zero-sound  contribution 
respectively.  In  the  following  analysis  we  will  investigate  the  forward  scattering  contri- 
bution from  the  particle-hole  excitation  and  the  zero-sound  mode,  and  show  that  ImE  is 
completely  renormalized. 
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Particle-hole  contribution.  The  first  term  of      is  real  [Eq.  (2.42d)]  and  does  not 
contribute  to  the  imaginary  part,  whereas  the  second  term  differs  from  the  second  order 
self-energy  term,      ,  by  only  a  sign.  Thus  we  can  immediately  infer  ImE^  from  the 
forward  scattering  part  of  ImE^  [Eq.  (2.22)], 

ImEj  =  -taE«=-^-l„g^.  (2.44) 

The  imaginary  part  of  Ep  and      in  the  retarded  formahsm  is  given  by 

—  J  j^lmG^ie  +  n,k  +  q)lmgp  (2.45) 

ImE^   =   -3C/£^y^ImGo^(e  +  Q,k-t-q)Ima.  (2.46) 

where  for  q  <  pf,  Rell^  -m/2i:,  thus  for  weak  interaction  Um/2-K  <  1,  the  terms  Im^?p 
and  ImGc  will  be  equal  in  the  leading  order  of  interaction 

.  ^         1  UlmU^  1  C/Imn« 

Im^?p   =  - 


2  (1  -  i7Ren«)2  +  {uimmy  2  1  +  {uimmy 
1  t/imn^  1  f/Imn^ 


2(1  +  t/Ren«)2  +  {uimmy^    2 1  +  (f/imn«)2  ■ 

Substituting  the  value  for  ImG^  and  Imll^  from  Eqs.  (2.12)  and  (2.13)  in  Eqs.  (2.45) 
and  (2.46)  we  obtain 


ImGpH   =   ImEp  +  ImE^ 

2  2 

4nE^ 


1      Ef  /  2A 

log  ~m  +  Gi 


v?\e\  vw^kl 


(2.47) 


where 


Gi{x)  =  2  log  2  -  1/2  +  log        -  2Re        /^^     log  ^  +  ^}-^/' 

Jo    y/l-x/z        1-  y/1  -  x/z 

In  the  limit  |x|  >  1,  the  scaling  function  behaves  as  G/(x)  «  log        and  hence 


(2.48) 
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In  the  opposite  limit  |a;|  <  1,  function  G]{x)  vanishes  as  xln  |x|.  As  a  result,  net  ImSpji 
remains  finite  at  A  =  0,  and  for  \x\  <  1  (i.e.,  for  |A|  <C         it  behaves  as 


y2  ^2 

ImSSo  =  —  — 

47r  Ef 


'  Ef  2A,  |A| 
In  ^rfr  +  —  In  '  ' 


(2.50) 


Comparing  the  hmiting  forms  of  Eq.  (2.49)  and  Eq.  (2.50),  we  see  that  higher  order  terms 
in  u  simply  cut  the  logarithmic  divergence  in  ImEp^  for  |A|  <  tx^|e|.  To  logarithmic 
accuracy,  one  can  then  approximate  ImEpjj  by 

ImE^H  =  ^^|?-ln^,  (2.51) 
™     AttEf     \w\  ^  ' 

where  w  =  max  (|A|,  u^|e|). 

Zero-sound  contribution.  Unlike  the  particle-hole  region,  which  contributes  to  all  the 
terms  in  Eq.  (2.43),  the  zero-sound  mode  makes  a  contribution  only  to  ImE^ : 

ImS«_^5  =  U      ~  j  ^linG^{e  +  Q,  k  +  q)Imap-zs,  (2.52) 

where  the  subscript  p  —  ZS,  stands  for  the  contribution  to  the  charge  channel  from  the 
zero-sound  mode  region.  The  collective  mode  Im^^_^s  has  a  pole  given  by 

We  will  plug  the  form  of  interaction  given  in  Eq.  (2.53)  into  Eq.  (2.52)  and  perform  the 
angular  integration  on      to  obtain 


QdQ  27r 


)'  VMY      +  A)2 

—-^[6{n~cQ)-6{Q  +  cQ)] 
on 

167r3  7o     ^{vfQY  -  {-QcTW 
near  the  mass-shell  QcA  >  A^  and  (?  -      =  {Um/2-Kf.  After  some  manipulations,  on 
Eq.  (2.54),  we  obtain 
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Figure  2-4:  Scaling  functions  Fj  and  G/  as  a  function  of  x.  Note  the  strong  asymmetry  of 
Fj  about  X  =  0. 


The  above  result  was  for  e  >  0.  For  e  <  0  we  have 


QdQ 


2A   _  I  Um\2 
V  27r  ) 


(2.56) 


Qvf      ^  2-K 

Thus  for  any  e  near  the  mass  shell,  the  contribution  to  the  imaginary  part  of  the  self- 
energy  is  given  by 

QdQ 


^2jj2^2  r  /■min{l,vS} 


2n 


f 

Jo 


y 


/o  ^/x  -  y2J 

The  integral  of  Eq.  (2.57)  can  be  performed  using  standard  Mathematica  software,  the 
result  is 


(2.57) 


r 

rix)  =  2tv  / 
Jo 


min{l,v^} 


y^dy 
y/x-y"^ 


SttItVS,  for  0  <  X  <  1; 


^sin-i  1 


,  for  a:  >  1. 


(2.58) 


Here  x  =  and  for  z  >  1,  Fi{x)  k,  27r/5i/i.  From  the  plot  of  F/(x),  Fig.  2-4,  it  is 

clear  that  ImEf  ^  remains  finite  in  the  entire  range  of  x.  Combining  the  asymptotic  value 
of  Fi{x)  and  Gi{x)  for  x  >  1,  we  recover  the  3'"'^  and  higher  order  results  obtained  via  the 
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perturbative  process, 


20Ef 


20Ef 


VSvh^  IS 
20Ef 


It  should  be  noted  that  the  above  form  of  the  third  order  self-energy  is  vahd  only  for 
regions  far  away  from  the  mass  shell. 

Unhke  the  perturbative  case,  where  the  regions  near  the  mass  shell  are  inaccessible 
due  to  divergence,  no  such  divergence  takes  place  for  Fi{x).  At  |A|  =  u'^\e\/2,  Fj{x)  has 
a  maximum  and  on  the  mass  shell  (A  =  0),  ImEf^  vanishes.  The  vanishing  of  ImSf^ 
on  the  mass  shell  is  due  to  the  Cherenkov  type  restriction  where  the  quasi-particle  with 
velocity  Vp  <  c  cannot  emit  a  zero-sound  boson,  hence  the  zero-sound  mode  cannot  decay 
the  quasi-particle.  The  asymmetry  of  Fi{x),  Fi{x)  =  0  for  x  <  0,  can  be  understood  if  we 
refer  to  Eq.  (2.54),  for  e  >  0  the  integral  will  have  non-zero  contribution  as  long  as  A  >  0. 
For  a  detailed  analysis  let  us  consider  Eq.  2.52.  The  delta  functions  originating  from  ImE^ 
and  IrnGp  impose  severe  restrictions  on  the  phase  space  for  the  emission  of  bosons.  We 
find  the  following  condition  should  hold  for  non-zero  contribution  to  ImS^,^^  or  for  the 
emission  of  Bosons  (in  the  following  we  will  use       vf{1  +  u^/2)) 

VJF. 

c 

1 


e-ik   =   f^H  r^cos^ 

c 

A 

-Kcos^   -  ^ 


0< 


0< 


2  <  1 

^  2 


<   A  <  2Q  <  2£. 


<  ^<A<1 


4e   -    2  -  v?e 

This  condition  clearly  suggests  that  e  and  A  should  have  the  same  sign.  To  study  the 
decay  of  a  quasi-particle,  we  consider  a  fermion  with  energy  e  such  that  u^|e|/2  <  A  or 
a;  >  1,  in  this  case  all  fi's  satisfy  Vt/Ae  <  x/2.  As  x  is  increased  further  the  condition 
x/2  <  n/u^e  imposes  constraints  on  the  low-energy  bosons,  thus  fewer  and  fewer  bosons 
are  able  to  contribute  to  the  decay  processes,  and  hence  F[{x)  reduces.  For  A  such  that 
uh/2  >  A,  i.e,  x  <  1  decreasing  A  (or  x     0)  implies  the  frequency  range  available  to  the 
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bosons  that  contribute  in  the  decay  process  shrinks,  thus  F/(x)  decreases.  This  explains 
the  non-monotonic  variation  of  ImEf^  as  a  function  of  the  parameter  x  =  2A/u^e. 

An  important  thing  to  notice  is  that  the  contribution  from  the  zero-sound  region 
ImEf  5  is  small  compared  to  the  contribution  from  the  particle-hole  region, 

ImEfg  ^  Fijx) 

since  Fi{x)  ~  1,  ImSp^  is  log-dominant  with  respect  to  ImEfg. 

Finite-curvature  effects.  From  the  analysis  thus  far,  we  conclude  that  the  divergences 
present  from  the  3'"''  order  onwards  in  the  perturbative  expansion  of  ImS^  [Eq.  (2.27)] 
have  been  regularized  by  resumming  the  series.  Nevertheless,  there  still  remains  a  log 
divergence  at  second  order  due  to  the  extra  term  ImE^,  which  must  be  cut  off  due  by  the 
finite  curvature  of  the  fermionic  dispersion.  The  log-singularity  is  cut  at    /2m  ~  Ac  = 
/Ep,  and  we  obtain  from  Eq.  (2.44) 

Thus  all  the  terms  of  Eq.  (2.43)  are  renormalized  and  ImS  is  finite  at  the  mass  shell.  A 
natural  question  that  arises  is:  is  it  possible  to  remove  the  mass-shell  divergence  in  the 
perturbative  expansion  of  ImE^  [Eq.  (2.27)]  without  resorting  to  resummation  and  instead 
removing  the  divergence  by  including  the  cut-off  due  to  finite  cm-vature?  Let  us  insert  the 
cut-off  Ac  =  e^/Ep  in  Eq.  (2.27).  We  obtain 

This  series  does  not  converge  for  \e\  <  Wc,  where  Uc  =  u^Ep,  thus  the  finite  curvature  does 
not  lead  to  the  convergence  of  the  series  and  the  need  for  resummation  remains. 

For  A  »  Ac  curvature  effects  can  be  totally  neglected  (this  is  the  regime  that  has 
been  considered  in  our  work)  and  the  results  for  Imf  ^  and  Imp^:^  will  remain  the  same  as 
has  been  derived  by  the  resummation  procedure.  For  A  <  Ac  the  results  will  change, 
since  Ac  can  no  longer  be  ignored  in  the  integrations.  To  obtain  the  functional  form  all 
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one  has  to  do  is  replace  A  by  Ac  in  the  functions 

/2Ae 


Fi 
Gi 


2A, 


Gi 


2|g| 
2\e\ 


We  will  use  this  information  to  obtain  ImSp^  and  ImSf^  on  and  near  the  mass  shell.  Let 
us  consider  the  limit  |A|  <^  Ac.  Following  our  earlier  analysis  we  will  make  the  necessary 
changes  ImE^^(|A|)  — >•  ImE^jj(|Ac|),  obtaining 


ImE^^ 


log 

Ef 

+  G/ 

f2Ac 

9 1  1 

U''\£\ 

log 

Ef 

+  Gi 

^2|e|\ 

AttEf 

9 1  1 

U''\£\ 

^Uc  ) 

As  a  reminder  the  scaling  function:  Gi{2\£\/u!c)  =  Gj{x)  behaves  as 

Gi{x)  = 
Prom  this  result  we  obtain 


xln|x|,    for   X  <C  1,    |e|  "C 
In  A,    for    X  >  1,    lei  >  Uc- 


(2.61) 


^InJ^,    for  |£|«a;c 


"      In  — 


(2.62) 


e  >  cjc- 


There  is  no  such  scaling  relation  for  ImEf^(A)  ImEf^(Ac)  and  the  functional  form 
remains  the  same  for  either  limit,  £  >  Wc  or  e  <  Wc, 


ImSf 


1 

4rr£p^°^M- 


The  zero-sound  contribution  to  ImEf^  obeys  a  scaling  relation  governed  by  the  scaling 
function  Fi{x).  For  x  <  1,  F/(x)  scales  as  STr^xVS  and  for  x  »  1  scales  as  2'K/by/x,  this 
then  implies 


ImEfs  = 


A-kEf 


for  |e|<Cu;c; 

— -f ,      |£|  3>  UJc- 


(2.63) 
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Adding  up  contributions  from  all  forward  scattering  terms,  the  final  form  reads 

logu^l,    for    l^l  <^  Uci 


ImS«  =  ImS^^  +  ImEf,  +  ImE«  =     '^"^l  /    "  '  (2.64) 

Notice  since         ~  1  for  all  ranges  of  x,  the  contribution  from  ImEp;^  and  ImE^  will  be 
dominant  by  a  large  log  term  compared  to  ImEf  and  has  therefore  not  been  included  in 
Eq.  (2.64). 

ImS^  on  the  mass  shell:  final  result.  The  leading  contribution  to  the  ImE^  comes 
from  both  forward-  and  back-scattering  contributions.  We  have  successfully  renormalized 
the  forward-scattering  contribution  and  as  discussed  before  no  such  renormalization 
scheme  is  required  for  the  back-scattering  contribution,  so  the  total  contribution  is 


ImS^  =  ImS«  +  lmS«  =  ^ln^ 

znEf  \e\ 


5,    for  |c|<a;c; 
1,    |e|  >  ujc. 


(2.65) 


Thus  near  the  mass  shell  ImE  behaves  as  the  usual  log 
2.3.5    Renormalized  Real  Part  of  the  Self-Energy 

In  the  following  analysis,  we  will  derive  the  real  part  of  the  self-energy  near  the  mass 
shell  from  the  back-scattering  process  by  making  use  of  the  Kramers-Kronig  relationship 
between  the  real  and  the  imaginary  part  of  the  retarded  self-energy.  A  detailed  calculation 
of  the  zero-sound  contribution  to  the  real  part  of  the  self-energy  is  also  presented,  where 
instead  of  the  Kramers-Kronig  relationship  we  have  considered  the  explicit  form  of  the 
ReE  in  the  retarded  formalism. 

Back  scattering.  In  an  earlier  analysis  we  obtained  ImE^^.^  and  ImS^.  A  sum  of 
these  two  back-scattering  contributions  is 

,,2  ^2 

^"^^(^^  =  4;^^1«SR-  (2.66) 
where  14^  is  an  upper  cut-off  typically  of  the  order  of  Ep.  We  will  make  use  of  the 
Kramers-Kronig  relation  to  obtain  the  real  part  of  the  self-energy  from  the  imaginary 
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paxt. 

n    J  E-e  87r4u2     w^ooj^  -  ) 

Focussing  only  on  the  cut-ofF  independent  and  nonanalytic  term,  we  obtain 

^'        Stt^^    w^ooJq  ^E\  E^-e^J 

lim  / 


mU'^e\e\^  /-^/I'l  ,  logx 
— — -5— lim   /  dx-r — - 

elel 


8 

Forward  scattering.  The  real  part  of  the  self-energy  consists  of  three  contributions: 
from  the  remainder  term  (S^),  from  the  particle-hole  continuum  (Epn),  and  from  the 
collective  mode  (Sfg)-  The  contribution  from  the  remainder  term  and  the  particle-hole 
term  is  extracted  similar  to  the  previous  approach,  that  is  through  the  Kramers-Kronig 
relation,  we  obtain  for  the  remainder  term 


(2.67) 


ReSf,  =  ^^+0(A2logA)  (2.68) 


and  for  the  particle-hole  term 


adding  up  the  contributions  from  Eqs.  (2.68)  and  (2.69)  we  obtain 

ReE« +ReS^H  =  -^-  (2.70) 

Forward  scattering:  zero-sound  contribution.  In  this  Section,  we  will  evaluate  the 
nonanalytic  e\e\  contribution  from  the  zero-sound  mode.  The  explicit  form  of  ReS  in  the 
retarded  formahsm  is 

^   '  J  (27r)2         TT         27r  ,  uj  +  ei  -  e 

tanh-+coth- 
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where  Im^^  is  given  by  Eq.  (2.52).  Following  the  usual  angular  integration,  we  perform 
the  dpi  integration  utilizing  the  delta  function  term  from  ImQp.  This  yields 

yoo  roo  ^2  /  ^  ^  \ 

ReE»(.,rt  =  a  I   ^  _     ^         ^      (tanh  ^  +  coth  -j 


,0  r 

a  j     djjj  I 

J—oo       J —oo 


OEx  — ,    tanh  — — — h  coth 


7(^)2 -(£i+A)2(a;  +  ei)  V         2r  2T 

where  a  =  u^UlAi:.  For  T  =  0,  the  trigonometric  functions  impose  constraints  on  e\  and 
u;.  Taking  the  zero  temperature  limit, 

ReEofe,  p)   =  a       duj       dei — .   

+a  /     da;  /  dei- 

J —  oo        J —c 


c 

=   a       dijj  i  dsi 


-a  I    duj  I  de\ 


Since  we  are  working  in  the  near  vicinity  of  the  mass  shell  the  condition  e  !»  A  holds, 
thus 


dw  j  de 


'   2u^  


c 

oo  ne 


2a;2 


=  a       duj  dei 
Jo  Jo 

+a  /    du      dei  ,  

Jo       Jo       ^{^r-{-ei  +  Anuj-€i) 

=  h  +  h. 

Evaluation  of  both  7i  and  I2  integrals  is  similar,  we  will  choose  to  integrate  I2  and  simply 
state  the  result  for  Ii.  We  impose  an  upper-cutoff  n  for  the  fi-equency  integration  (this 
constraint  represents  the  damping  of  zero-sound  modes  at  very  high  energies),  and  a  lower 
cut-off  at  ui  to  keep  the  integral  real,  contributions  from  the  high  frequency  regions  will  be 
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Figure  2-5:  Scaling  function  Fr.  Note  the  strong  asynmietry  of  Fr  about  x  =  0. 


irrelevant  and  will  not  be  considered. 


Q  /  (fel  / 
^0  Jq 


djjj 


>/(^)^-(-£i  +  A)2(a;-£0 


2c  r  e^log4 

=   -a—  1  — 

p  [    4  4 


+2a 


/ 


0  V(^)'-(-^i  +  A)^ 


loe 


-2 


(ei  -  A)2  -  (p/c)2c^ei  1 


It  can  be  shown  that  the  term  inside  the  log  and  the  square  root  in  the  denominator  are 
purely  imaginary,  this  further  simplifies  the  above  expression 


h  =  -a 


2c 
V 


Similarly  we  obtain  for  /j 


I\  =  —a 


2c 
P 


£:^log4 
"T^  4 


log  4 
4 


+  TTCiRe 


v'ti%?-2eiA' 


Adding  up  /i  and  /2  we  obtain 


/i  +  /2  =  TraRe  /  —^=jfiL=  +  naRe  f 
Jo  ^/u^ej-2eiA  Jo 


lA  Jo  Vu%f  +  2eiA 
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The  first  term  is  more  dominant  and  we  will  retain  it,  so  the  final  expression  for  ReS^^  is 


ReEfs 


47r 
8E, 


Re 

qdq 

Jo  sju^ 

-  2A/vFq 

-Fr 

\  u^e  J  ' 

(2.71) 


where  the  scaling  function  Fr(x)  is, 
Fr{x)  =  Re 


(1  +  -x)  +  -x^  log  j=— 


(2.72) 


As  shown  in  Fig.  2-5,  the  scaling  function  is  asymmetric  with  respect  to  x  — >  — x,  for 
X  >  1  the  function  vanishes  and  for  large  and  negative  x, 

Fr{x)  4Re- 


5v— X 

The  real  part  of  the  self-energy  away  from  the  mass  shell  is 


ReS5 


x/2-3 


zs 


Re 


20  Ep' 

the  ZS  contribution  away  from  the  mass  shell  is  smaller  by  a  factor  of  u  compared  to  the 
\P'e\e\  contribution  to  the  self-energy  from  the  back-scattering  process.  However,  right  on 
the  mass  shell  the  interaction  with  the  zero-sound  mode  is  strongly  enhanced,  the  function 
Fr{.^)      1  and  ReEf^  =  \P'e\e\l%EY,  this  result  is  of  the  same  order  as  the  contribution 
from  the  PH  region.  The  small  x  expansion  of  Fr{x)  is  given  by 


3x^  1 
Fr(x)  =  1  +  X  +  -—  log  -  +  ■  ■ 
4  X 


(2.73) 


Using  this  expansion  in  Eq.  (2.71)  we  obtain  to  linear  order  in  A,  the  following  correction 
to  the  self-energy 


8  Ep 


+ 


Adding  up  contributions  from  Eqs.  (2.70)  and  (2.74)  we  obtain 

ReE^  =  ReE^„  +  ReEf.  +  ReEf^  = 


(2.74) 


(2.75) 


49 

Final  result  for  the  real  part  of  the  self-energy.  Collecting  up  all  nonanalytic 
contributions  to  the  ReSf  from  the  forward  scattering  processes, 


ReEf  =  ReEf  +  ReEf^  +  ReEf  5 

^    ^  2  11 

W  e\e\ 


(2.76) 


<-    -  8  AEf 

we  find  the  above  result  from  the  forward-scattering  process  differs  from  the  back- 
scattering  result  by  a  sign.  Thus  near  the  mass  shell 

ReE^  =  ReEi  +  ReEf  =  0{u^A^  log  A).  (2.77) 

The  conclusion  from  this  result  is  that,  for  a  short  range  potential  the  nonanalytic  u^e\e\ 

term  is  absent. 

2.3.6    Spectral  Function 

Having  determined  ImE^  and  ReE^,  we  are  now  in  a  position  to  evaluate  the  spectral 
function.  The  spectral  function  in  terms  of  the  Green's  function  is  expressed  as 

A{e,k)  =  -hmG'^{e,k) 

1  ImE^(5,fc)  

TT  [A  +  ReE«(£,  fc)]2  +  [ImE«(e,  k)]'^ '  ^  ^  ^ 

We  will  plot  the  spectral  function  as  a  function  of  x  =  2A/u^e,  keeping  e  fixed  and 
varying  A  —  e  —     through  its  momentum  dependence.  We  will  also  consider  a  situation 
where  A  ~         >  Ac  (as  a  reminder  Ac  is  the  energy  scale  at  which  the  curvature 
effects  become  important),  so  that  we  can  neglect  curvature  effects  on  the  spectral 
function.  Near  the  mass  shell,  ReE^  in  the  denominator  can  be  neglected  compared  to  A 
as  ReE^  ^  u^A^  log  A  <  u^uP'  log  |a;|.  Thus  the  spectral  function  reduces  to 

^^''^^   -    7rA2  +  [ImE«(e,/t)]2-  ^^"^^^ 

The  variation  of  ImE'*  as  a  function  of  e  or  ^fe  cannot  be  neglected,  and  is  responsible  for 
non-Lorentzian  type  features  in  the  spectral  function.  The  imaginary  part  of  the  retarded 
self-energy  is  a  sum  of  four  different  contributions, 

ImE^  =  ImEf  +  ImE^^  +  ImEf^  +  ImEfs, 
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Figure  2-6:  The  log-plot  of  spectral  function  A{e,  k)  in  units  of  I/tt'^u^Ef.  A{£,  k)  is  plot- 
ted as  a  function  of  x  =  2(e  -  ik)lu^e  for  \og{EFlu^\e\)  =  2  and  el2nEF  =  0.025.  A  kink 
at  X  =  1  is  due  to  the  interaction  of  fermions  with  the  zero-sound  mode.  Inset:  part  of  the 
spectral  function  Ai{e,k)  for  £/2'kEf  =  0.25.  A  maximum  in  Ai  sX  x  =  1  gives  rise  to  a 
kink  in  total  A  (main  panel). 

where  the  last  three  terms  are  due  to  the  forward  scattering  contribution.  Near  the 
mass  shell,  both  ImEg  and  ImSf^  are  independent  of  A  and  to  leading  order  the  sum 
IinSfx  +         vanishes,  thus  ImS^  shows  dependence  on  A  through  the  remaining  terms, 


ImS^ 


zs 


ImE^ 


PH 


Fj{x) 


AttEf 

The  spectral  function  acquires  becomes 

1  logJ^  +  G/(x)  +  F,(x) 


(2.80) 


A{e,k)  = 


TT^U^El 


+  (a)  h  ( Jr)  +  Gi{x)  +  Fi{x) 


(2.81) 


The  scaling  function  Gi{x)  exhibits  weak  dependence  on  x  for  x  ~  1.  Whereas,  Fi{x) 
is  strongly  dependent  on  x  with  a  sharp  peak  at  x  =  1.  The  spectral  function  shows  a 
narrow  peak  at  x  =  0  and  a  well-pronounced  kink  at  z  =  1,  the  kink  corresponding  to  the 
sharp  peak  of  F{x).  A  plot  of  the  spectral  function  A{x)  as  a  function  of  x  is  shown  in 
Fig.  2-6. 
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For  typical  A  ^  w^lej  ov  x  ^  I  ,  logiEp/u^lel)  >  [Gi{x)  +  Fi{x)].  Therefore,  we  can 
split  the  spectral-function  into  two  parts:  a  part  which  describes  the  regular  quasi-particle 
peak  (Aq)  and  a  part  with  scaling  behavior  (Ai) 

A{e,k)   =  Ao{6,k)  +  Ai{e,k), 

1  \og{EF/u^\e\) 


Ao{e,  k) 


Me,k)   =  -^[Grix)  +  F,{x)]-  " (2.82) 


''''''^^^'+(2^fclog(EWn2|£|))" 

-  {^logjEpm)) 
(x^+{^log{EFm)) 

Near  x  =  1 ,     (x)  is  a  smooth  function  whereas  Ai  (x)  shows  a  rapid  variation.  This 
variation  associated  with  Ai{x)  is  due  to  the  fact  that,  for  |A|  >  u^lel  (x  >  1),  a  fermion 
with  frequency  e  >  0  can  emit  a  zero-sound  boson  of  any  frequency  from  0  <  Q  <  e 
whereas  for  |A|  <  u^\e\,  a  Cherenkov  type  restriction  makes  it  impossible  to  emit  and 
absorb  bosons  with  frequencies  above  |A|/(1  —  vp/c)  —  A/u^. 
2.3.7    Coulomb  Potential  '  ■ 

Next,  we  will  consider  the  interaction  between  the  fermions  and  the  plasmon  modes 
and  show  that  the  kink-like  behavior  in  the  spectral  function,  which  was  obtained  for 
fermions  interacting  with  the  zero-sound  mode  in  the  case  of  a  short-range  potential,  is 
also  present  for  the  case  of  electrons  interacting  with  plasmon  modes.  Herein,  we  will 
consider  the  low  energy  limit:  with  |e|,^fc  <  ujpu  where  lOj,i  =  y^rgEp  is  the  maximum 
(undamped)  plasmon  frequency.  Near  the  plasmon  mode  the  potential  is  given  by 

V{q) 


Viuj,q) 


l-Viq)U{q,io) 
^-Ke^u?  jq 


(a;  +  zO+)2-fio(g)2' 
^  27re^mvp 

{oj  +  iO+)2  -  Qo(g)2  ^^-^^^ 

where  Qo(q)  =  (e^mvlqy/'^  is  the  plasmon  dispersion  and  we  have  used  the  following  form 
for  the  polarization  operator: 

tR         ITT'  I .  w 


n«  =  --  1 
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Using  Eq.  (2.83)  as  the  potential,  we  can  calculate  the  imaginary  part  of  the  self  energy 
(we  will  keep  e  >  0),  i.e.,  the  damping  due  to  the  interaction  of  the  quasi-particle  with  the 
plasmons,  r 

-       dfl  j  (fqS{£-(k-^ -vfQ  cos  e)6{Q:^ -Qliq)),  (2.84) 

where  the  first  (5-function  is  due  to  the  Green's  function  and  the  second  is  due  to  ImV^. 
The  second  delta  function  forces  the  transferred  boson  momentum  to  be  small,  i.e., 
q  <C  1^^0(9)  =  ^\/vf,  thus  we  can  neglect  the  momentum  dependence  in  the  first  delta 
function.  The  rest  of  the  analysis  is  trivial  and  we  obtain 

ImE«(£,6)   =   7r^^'~^^'    for    0  <     <  e 
hp 

=   0    otherwise  (2.85) 

For  a  fixed  value  of  e,  ImE^  has  a  kink  at     =  0,  also  the  kink  is  independent  of  the 
coupling  constant.  A  kink  in  the  spectral  function  could,  in  principle  be  detected  in 
photoemission  experiments  on  layered  materials  or  in  a  momentum  conserving  tunneling 
between  two  parallel  2D  gases  [82]. 
2.3.8    Corrections  to  the  Tunneling  Density  of  States 

The  presence  of  |£|A  term  in  the  real  part  of  the  self-energy  of  the  particle-hole 
[Eq.  (2.70)]  and  the  ZS  [Eq.  (2.74)]  contributions  can  in  principle  give  rise  to  a  nonan- 
alytic  correction  to  the  tunnehng  density  of  states.  The  tunnehng  density  of  states  is 
represented  as 

Ni')  =  -l  J  -0^^^G^i^^k).  (2.86) 

The  term 

a\£\—  =  a\£\—^ 

H/p  Hip 

in  the  real  part  of  the  self-energy  is  absorbed  by  the  Z-factor  in  the  Green  function, 
yielding  a  simple  form  for  the  interacting  Green's  function  in  terms  of  the  free  Green's 
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function 


{k,w)  = 


1  1  Z{e) 


(2.87) 


The  tunneling  density  of  state  for  a  free  system  is  m/n.  Thus  we  can  deduce  the  tunnehng 
density  of  states  for  an  interacting  system  to  be 


This  would  be  the  case  if  a  nonanalytic  |e|A  were  present  in  the  real  part  of  the  self- 
energy.  However,  we  see  from  Eqs.  (2.70)  and  (2.74)  that  the  |£|A  terms  cancel  out, 
moreover,  the  back-scattering  term  of  the  ReS  does  not  contain  a  |£:|A  term.  Thus  the 
tunneling  density  of  state  does  not  get  any  nonanalytic  correction.  The  above  result, 
vaUd  for  a  contact  potential,  will  also  hold  for  a  finite  range  potential,  since  as  before 
the  nonanalytic  |e|A  terms  cancel  between  the  different  forward  scattering  contributions, 
whereas  the  term  is  absent  in  the  back-scattering  contribution. 


We  will  conclude  this  Chapter  by  restating  our  main  results.  Second-order  pertur- 
bation theory  yields  two  interesting  results,  namely,  a  nonanalytic  term  in  the  imaginary 
part  of  the  self-energy  and  a  logarithmic  singularity  on  the  mass  shell  for  fermions  with 
a  Unearized  dispersion.  The  nonanalytic  contributions  originate  from  processes  involving 
two  incoming  particles  of  nearly  opposite  momenta  undergoing  either  a.  q  ^  0  or  q  ^  2kF 
scattering,  and  also  from  the  scattering  process  between  two  incoming  particles  with 
nearly  identical  momentum  that  undergo  q  k.  Q  momentum  transfer.  The  latter  process 
(the  54  process  in  ID  terminology)  is  also  responsible  for  a  logarithmic  singularity  on  the 
mass  shell,  indicating  that  the  quasi-particles  are  ill-defined.  We  find  that  at  higher  orders 
in  the  interaction  the  singularity  is  enhanced  and  becomes  a  power  law,  causing  a  formal 
breakdown  of  the  perturbation  theory.  These  results,  which  are  obtained  for  fermions  with 
a  linear  dispersion,  are  modified  by  restoring  finite  curvature  of  the  dispersion.  The  sin- 
gularities from  the  individual  terms  in  the  perturbation  theory  are  removed,  yet  the  series 
itself  remain  divergent.  A  re-summation,  performed  by  selecting  the  maximally  divergent 
diagrams,  removes  these  singularities  even  for  a  linearized  dispersion,  with  the  result  that 


(2.88) 


2.4  Conclusion 


Fermi-liquid  theory  is  restored.  However,  we  demonstrate  that  the  interaction  of  fermions 
gives  rise  to  a  kink  in  the  spectral  function  at  A  =  The  non-Lorentzian  shape  of 

the  spectral  function  should  be  amenable  to  a  direct  check  in  photoemission  measurements 
or  in  momentum-conserved-tunnehng  between  two  parallel  2D  electron  gases. 


CHAPTER  3 
SPECIFIC  HEAT  OF  A  2D  FERMI  LIQUID 

3.1  Introduction 

Thermodynamic  quantities  in  a  Fermi  gas  form  a  regular  and  analytic  series  as  a 
function  of  temperature  T  and  momentum  q.  For  example, 

C{T)/T  =  7  +  aoT^  +  aiT^  +  .. 
Xs{T,q  =  0)   -   x1{0)  +  boT'  +  .. 

Xs{T^O,q)   =   x°(0)  +  Cog2  +  ..  (3.1) 

where  the  even  powers  of  temperature  T  are  due  to  the  particle-hole  symmetry  of  the 
Fermi  function  about  the  Fermi  level  and  the  even  powers  of  q  are  due  to  the  quadratic 
dispersion  of  the  free  particles.  This  result,  valid  for  a  Fermi  gas,  does  not  necessarily  hold 
for  interacting  fermions.  Fermi-liquid  theory  predicts  only  that 

C(T)/r  =  7*, 
Xs{T,q)  =  xm, 

where  the  constants  7*  and  x*,  differ  from  the  free  Fermi  gas  result  and  in  general 
depend  on  the  interaction  parameters.  To  obtain  higher  orders  in  the  temperature-  and 
momentum-expansion  of  Cy  and  Xs  for  an  interacting  system,  the  usual  approach  is 
to  consider  microscopic  models  (for  example,  fermions  with  short-range  repulsion)  and 
employ  perturbation  theory.  Corrections  to  the  linear-in  temperature  specific-heat  term 
are  conveniently  obtained  by  calculating  the  self-energy  on  the  mass  shell,  and  relating 
the  self-energy  to  the  specific  heat.  A  nonanalytic  correction  to  the  self-energy  leads  to 
nonanalytic  terms  in  the  thermodynamic  quantities,  for  example,  in  3D  an  e^log  |£|  term 
in  the  real  part  of  the  self-energy  leads  to  a     log  T  term  and  in  2D  an  £\e\  term  in  the 
self-energy  gives  a  nonanalytic      term  [19,  20,  21,  22,  23,  24]. 
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In  the  last  Chapter  we  investigated  the  nonanalytic  contributions  to  the  self-energy 
from  perturbative  and  nonperturbative  processes.  It  was  shown  that  upon  resummation, 
the  nonperturbative  contributions  to  ReE^  become  of  the  order      near  the  mass  shell 
and  exactly  cancel  the  second  order  U'^ele]  term.  This  would  then  mean  an  absence 
of  nonanalytic      terms  in  the  specific  heat  at  second  order  in  the  interaction.  We 
investigate  the  role  of  contributions  of  nonperturbative  processes  to  the  specific  heat  and 
find  that  the  imaginary  part  of  the  self-energy  contributes  to  the  nonanalytic  term  in 
the  specific  heat,  in  such  a  way  that  it  cancels  the  contribution  from  the  nonperturbative 
real  part  of  the  self-energy,  thus  keeping  the  second-order  result  intact.  We  confirm 
our  result  in  three  different  ways.  In  Section  3.2  a  direct  calculation  of  the  specific 
heat  from  the  self-energy  is  performed,  and  we  find  that  the  contributions  from  the 
nonperturbative  processes  cancel  out.  In  Section  3.3  we  use  the  Luttinger-Ward  formula 
for  the  thermodynamic  potential  to  show  using  real  frequencies  (Section  3.3.1)  and 
imaginary  frequencies  (Section  3.3.2)  that  the  nonperturbative  processes  do  not  contribute 
to  the      term  in  the  specific  heat  through  order  U"^.  We  will  consider  the  long-range 
Coulomb  potential  in  Section  3.4  and  show  that  the  plasmons  (nonperturbative  processes) 
make  an  interaction  independent  contribution  to  the  specific  heat.  We  end  the  Chapter 
with  a  conclusion  in  Section  3.5. 


We  will  use  the  result  of  Ref.  [16]  which  relates  the  thermodynamic  potential  to  the 
Green's  function  and  estabhsh  a  relationship  between  the  self-energy  and  the  specific  heat. 
The  thermodynamic  potential  reads, 


where  the  Green's  function  is  expressed  in  terms  of  Matsubara  frequencies  ujm  and  has  no 
additional  temperature  dependence.  We  can  convert  the  summation  over  frequencies  into 
a  contour  integration  obtaining 


3.2    Specific  Heat  Calculation  From  the  Self-Energy 


(3.2) 


The  thermodynamic  potential  and  the  specific  heat  are  related  via 


C  =  -T 

using  Eq.  (3.2)  and  Eq.  (3.3)  one  obtains, 

2T2   a  ri    /■  (Pk 
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(3.3) 


k) 


(3.4) 


For  the  case  of  a  weak  interaction  such  that  |E|  -Co;,  one  can  expand  the  Green's  function 
in  Eq.  (3.4)  in  powers  of  E.  To  the  lowest  order  in  S  we  obtain 


C{T)  =  C^g{T)  +  5C{T), 


(3.5) 


where  Cfg(T')  =  mnT/3  is  the  specific  heat  for  free  fermions  in  2D  and  SC{T)  is  given  by 


5C{T)/T- 


2  d_ 


(3.6) 


There  are  two  contributions  to  6C{T)-one  from  Re  and  another  from  ImS^-which  are 
labelled  Ci{T)  and  C^iT),  respectively: 


SC{T) 
Ci{T)/T 


Ci  [T)  +  C2  (T) ; 


C2{T)/T  =  - 
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TT  dT 
2  d 


(Lju^lmG^iu,  A;)ReE«(a;,  A;) 
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n  dT 


f  J  J2^  J_ 

i  rd..§^P^ImE>,.)  . 

T  J  (27r)^  5a;     w-Cfc         \  \ 


(3.7a) 

(3.7b) 
(3.7c) 


Eq.  (3.7b)  implies  that  to  calculate  the  specific  heat  contribution  from  the  real  part  of 
the  self-energy,  one  needs  only  the  mass-shell  form  of  the  real  part  of  the  self-energy.  On 
the  other  hand,  Eq.  (3.7c)  implies  that  the  contribution  from  the  imaginary  part  of  the 
self-energy  is  zero  if  the  specific  heat  is  independent  of  momentum  k,  since  the  integral 

1 


/ 


(fk- 


vanishes  if  one  assumes  the  density  of  states  to  be  constant. 

In  the  discussion  to  follow  we  will  show  that  the  nonperturbative  contribution  to  the 
f/^T^  term  in  the  specific  heat  vanishes.  Thus  contrary  to  earUer  results  [21],  the  t/^T^ 
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term  is  not  modified  by  the  additional  contribution  from  the  nonperturbative  process  and 
remains  exactly  the  same  as  obtained  from  the  second  order  perturbative  calculation  [22]. 
3.2.1     Zero-Sound  Mode  Contribution  to  the  Real  Part  of  the  Self-Energy 

The  zero-sound  mode  contributes  both  to  the  real  and  imaginary  parts  of  the  self- 
energy.  In  this  Section,  we  will  calculate  nonanalytic  contributions  from  the  real  part 
of  the  self-energy  ReEf^  =  u^Lo\u!\/8Ef  to  the  specific  heat.  Notice,  for  calculating  the 
contributions  to  the  specific  heat  from  the  real  part  of  the  self-energy  we  need  only  the 
mass-shell  form  of  the  self-energy.  The  contribution  is, 


<5Ci(r)/r  = 
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iEfdT 
ANu'^T 

9  ) 
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where 
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Jo  cosl 


dx 


9C(3) 


(3.8) 


(3.9) 


cosh  {x)  8 

3.2.2     Zero  Sound  Mode  Contribution  to  the  Imaginary  Part  of  the  Self-Energy 

For  local  self-energy  (i.e.,  independent  of  momentum)  the  imaginary  part  of  the 
self-energy  does  not  contribute  to  the  specific  heat.  In  general,  the  imaginary  part  of  the 
self-energy  need  not  be  local,  this  is  the  case  for  ImEf^.  Near  the  mass  shell  |a;| 
thus 


2y3  M/c 
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QdQ 


2A 

Qvf 


'0  ^Jsgn{uj) 

—Re  ^"^^ 
m  Jo 


2A 
QvF 


(3.10) 
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The  specific  heat  contribution  from  the  above  term  can  be  written  as, 
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where  we  have  made  the  substitution  z  =  2A/vf.  We  will  first  calculate  the  following 
integral, 


(3.11) 


f^ReT 
Jo     z  Jo 
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By  splitting  into  regions  in  the  (2,  Q)  plane  where  the  integrand  is  real,  we  obtain 
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Using  this  value  back  in  Eq.  (3.11)  yields 
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2uvp ' 


ANu'^T 


TTVt 


(3.12) 


Adding  up  the  two  contributions  from  the  zero-sound  region,  Eqs.  (3.8)and  (3.12),  we 
obtain 


SC2  SCi_ 


(3.13) 


thus  we  find  the  zero  sound  makes  no  contribution  to  the  leading  nonanalytic  term  in  the 
specific  heat. 
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3.2.3    Finite  Range  Potential 

The  presence  of  a  finite-range  potential  U  {q)  does  not  modify  our  earlier  result 
regarding  the  cancellation  of  the  nonperturbative  contributions.  This  is  because  the 
nonperturbative  terms  arise  due  to  forward  scattering  processes  with  small  momentum 
transfer  hence  U{q)  can  be  replaced  by  U{0).  The  analysis  for  the  cancellation  of  non- 
perturbative contribution  to  the  f/^T^  term  in  the  specific  remains  the  same.  As  a  result, 
the  nonanalytic  contribution  to  the  C/^T^  term  in  the  specific  heat  is  the  same  as  that 
obtained  by  the  perturbative  approach  and  is  given  by  [22] 


6C{T)  9C(3) 


2 

.2^. 

U''{0)  +  U'\2kF)  -  UiO)U{2kF)  (3.14) 


3.3    Specific  Heat  Calculation  from  the  Luttinger-Ward  Formalism 
The  approach  used  in  the  previous  Section  to  show  the  absence  of  a  nonperturbative 
contribution  to  the  7^  term  is  quite  involved,  so  we  will  use  an  alternative  method,  which 
is  easier  to  execute,  and  in  the  process  double  check  our  earher  result.  We  will  use  the 
thermodynamic  potential  which  is  expressed  in  terms  of  the  Luttinger-Ward  formalism  to 
calculate  the  specific  heat.  The  most  commonly  cited  expression  in  the  literature  for  the 
thermodynamic  potential  of  a  system  of  interacting  fermions  with  a  potential  gV{x)  is 

^(9)  =  ^0  +  ^  r^Yl  exp{iUr,m{9',k,uJn)G{g',k,Un),  (3.15) 
P  Jo    9  - 

where  g  is  the  coupling  constant,  E  and  G  are  the  full  self-energy  and  Green's  function 
respectively  [83],  and  Qq  is  the  thermodynamic  constant  in  the  absence  of  mteractions 

^o  =  -jY,  log{l  +  exp[-/J(e,^  -  //)]}. 
ic 

As  shown  by  Luttinger  and  Ward  [84],  Eq.  (3.15)  can  be  written  in  a  compact  form 

2V 


2V  r 


+  y  E  E  (2^)  ^rn{iuJn,  k)G{ioJn,  k),  (3.16) 


where  in  the  last  expression  Sm(fe,w„)  stands  for  the  m^^  order  self-energy  diagram 
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Figure  3-2:  Diagrams  through  the  third  order  that  contribute  to  the  thermodynamic 
potential. 

determined  from  the  skeleton  diagram  with  fully  interacting  Green's  function.  An  example 
of  a  second  order  skeleton  diagram  is  shown  in  Fig.  3-1. 

To  calculate  the  specific  heat  from  the  thermodynamic  potential  we  will  be  interested 
in  those  diagrams,  that,  at  each  order  of  interaction,  have  the  maximum  number  of 
particle-hole  bubbles  (these  are  the  most  divergent  diagrams  of  the  thermodynamic  series). 
Typical  examples  of  third-order  diagrams  are  shown  in  Fig.  3-2(3).  It  might  seem  that 
only  one  of  the  diagrams  in  Fig.  3-2(3)  belong  to  the  category  of  "maximum  number  of 
particle-hole  bubbles" ,  however  a  more  careful  analysis  shows  that  the  rest  of  the  diagrams 
also  have  equal  number  of  bubbles,  though  they  are  'implicit'  (  they  are  obtained  by 
integrating  over  the  fermionic  degrees  of  freedom).  Furthermore,  it  is  easy  to  show  that 
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the  diagrams  of  the  kind  shown  in  Fig.  3-2  (which  are  "part"  of  the  skeleton  diagrams) 
are  obtained  from  the  last  term  of  Eq.  (3.16).  It  turns  out  that,  even  for  higher  orders, 
the  third  expression  is  the  only  one  which  is  relevant.  Following  the  prescription  described 
above,  we  obtain  the  following  series  for  the  thermodynamic  potential  of  a  short  range 
potential, 


no  +  vT 


(27r)2 
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5  3 


(3.17) 


3.3.1    Real-Frequency  Approach 

In  an  earlier  work  on  nonanalytic  contributions  to  the  specific  heat  [15,  22],  it  was 
shown  that  the  nonanalytic  U'^T^  contribution  to  the  specific  heat  arises  firom  two  regions 
in  momentum  space,  from  k  <^  kp  and  k  ^  2kF-  In  the  following  analysis,  we  re-sum 
higher-order  terms  so  as  to  obtain  the  nonperturbative  contributions  from  the  A;  <  jt^ 
regime.  We  can  split  the  thermodynamic  potential  into  two  parts,  a  part  containing 
contributions  coming  from  the  small  momentum  region  (summed  to  all  orders)  and  a  part 
which  has  contributions  from  the  large- momentum  region  (only  the  second-order  term). 
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(3.18) 


The  above  expression  for  the  thermodynamic  potential  can  be  rearranged  to  obtain  a 
compact  form 

spin-channel  charge-channel 
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{(f/n)2  +  -  iog(i  +  f/n)  +  ^  iog(i  -  t/n)}] 


nonperturbative  Contribution 


(3.19) 


where  the  terms  inside  the  curly  brackets  are  the  nonperturbative  contributions.  The 
summations  on  discrete  Matsubara  firequencies  can  be  converted  to  integrals  through  the 
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Figure  3-3:  Contour  for  summing  up  the  Log.  The  thick  hne  represents  the  branch  cut 
and  the  small  circles  represent  the  discrete  Matsubara  poles. 

usual  trick  employed  in  [16].  We  will  show  the  procedure  of  summation  for  the  log  terms 
of  Eq.  (3.19).  By  making  use  of  the  relation  between  the  Matsubara  dielectric  function 
D^{k,iWm)  and  the  retarded  and  advanced  functions 

D'^{k,iu  ^LO  +  iQ-)  =  D^{k,uj), 

the  discrete  sum  can  be  expressed  as  an  integral  with  contours,  Ci  and  C2,  as  shown  in 
Fig.  3-3.  The  discrete  summation  is 
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■^TJ2  hog{D''{k,iUn))   =  ■^T\ogD''{k,Q) 
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Jc  ^i'^Biuj)\og  D^ik,Lo),  (3.20) 


where  715(0;)  is  the  Bose  distribution  function.  The  retarded  and  the  advanced  function 
have  neither  poles  nor  zeroes  in  the  upper  and  lower  half  plane  respectively.  Also  the 
function,  dnBiuj)/dT,  decays  rapidly  as  w  ^  ±00.  These  properties  of  the  integrand 
allows  us  to  "open"  the  contour  and  express  the  integrals  in  terms  of  contours  d  and 
C'2.  Notice,  the  small  bumps  in  the  contours  C'l  and  C'2  cancel  the  first  term  of  the 
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above  equation.  So  the  above  equation  with  contours  C'l  and  C'2  minus  the  bumpy  part 
is 


d_ 
dT 


Next,  we  will  make  use  of  the  property  of  the  log  and  the  property  of  the  function  on  the 
real  axis       =  D^* ,  to  obtain  a  simplified  form.  The  contributions  from  the  spin  and 
charge  channels  wiU  be  evaluated  using, 
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The  contributions  that  will  be  considered  are  those  from  the  spin  channel,  the  charge 
channel  and  the  Y\?  term.  Out  of  these  only  the  charge  channel  will  get  contributions  from 
the  zero-sound  mode  region.  Let  us  first  consider  the  spin  channel: 
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where  the  polarization  operator  in  the  particle-hole  region  is, 


Ll  +  f/Ren«J' 


(3.22) 


n«  =  -^(i  +  i-^  


w 


The  above  term  does  not  get  any  contribution  from  the  zero-sound  mode  region.  We 
will  obtain  the  contributions  from  the  paxticle-hole  region  under  the  assumption  of  weak 
coupling  mU  <C  1 
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We  are  interested  in  the  term,  hence  we  will  only  be  interested  in  the  term  from  the 
dk"^  integration.  The  XiP'  term  from  the  dk"^  integration  is 
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dfc^tan  ^ 

vp 

Therefore,  the      contribution  is 
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contributions  from  the  charge  channel  (particle  hole  region)  can  be  obtained  by 
utilizing  the  same  approach  as  for  the  spin  channel.  We  obtain 

iw„,k 

Unlike  the  spin  channel,  the  charge  channel  does  get  a  contribution  from  the  zero-sound 
mode  region 

where 
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Go 


the  integral  of  Eq.  (3.26)  is  easy  to  perform  as  the  argument  of  Qp  in  the  momentum 
region  u/c  <  k  <  u/vp  is  constant.  In  this  region  Re(^p)  <  0  and  Im(^p)  — >  0-,  giving  an 
argument  of  -tt.  Thus, 
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where  c  =  vp/y/l  -  {Um/2'Ky  is  the  zero  sound  velocity.  This  gives 
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The  contributions  coming  from  the  small  momentum  region  to  the  second  order  term  is 
d  1  ^  f  cPk    /■+~  dw  2U^w 


d  1  ^  ^2.r2  f  /■+~  dw  2U^W         ,  ,  , 
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The  imaginary  part  of  (H^)^  is 


The  first  term  of  the  above  expression  gives  contributions  to  the  T  and  the  term,  the  5 
function  term  is  responsible  for  the  term 
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We  can  now  add  up  Eqs.  (3.24),  (3.25),  (3.28),  and  (3.31)  to  obtain  the  nonperturbative 
contribution  of  Eq.  (3.19).  The  result  is  that  the  nonperturbative  contributions  add 
up  to  zero,  thus  there  is  no  nonperturbative  contribution  to  the      term  in  the  specific 
heat,  and,  only  the  perturbative  second  order  in  the  interaction  strength,      terms  of 
Eq.  (3.19),  contribute  to  the  f/^T^  term  in  the  specific  heat.  We  will  show  below  that  the 
k  ^0  and  k  «  2kF  regions  give  equal  contributions  to  the  U^T'^  term  in  the  specific  heat. 
3.3.2    Matsubara  Formalism 

In  this  final  approach,  we  will  keep  Eq.  (3.17)  as  it  is  and  perform  the  summation 
directly  on  the  Matsubara  frequencies.  A  great  advantage  with  this  approach  is  that  it  is  a 
priori  apparent,  for  weak  interaction  u  <C  1,  that  we  do  not  need  to  consider  terms  higher 
than  the  second  order  term  for  the  [/^T^  contribution,  as  none  of  the  terms  in  the  series 
are  divergent  and  the  series  converges.  On  the  other  hand,  a  resummation  will  indeed  be 
required  if  the  interaction  is  strong,  u  ~  1,  or  long  ranged.  In  the  following  discussion  we 
will  consider  the  two  cases:  weak  potential  u  <  1  and  the  strong  potential  u  ~  1.  We  will 
then  extend  our  result  to  include  a  potential  with  an  explicit  momentum  dependence. 

Contact  interaction:  u<^l.  For  the  weak  potential  it  suffices  to  consider  terms 
through  second  order  in  the  interaction.  The  first  order  term  does  not  contribute  to  the 
nonanalytic  temperature  dependence  since  the  nonanalyticity  arises  only  from  the  2"'* 
order  onwards. 

We  will  calculate  the  contribution  fi-om  the  11^  term  which  originates  from  A;  <  A;^ 
regions,  followed  by  contributions  ft-om  the  A;     2kF  region.  The  Matsubara  polarization 
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operator  in  the  small  frequency  and  momentum  region  is 
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Using  this  form  of  the  polarization  operator  and  performing  the  momentum  integration, 
we  obtain 
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where  iV  >  1  is  the  upper  cutoff  and  X  =  {kFVFl2-KTf  >  1.  Notice,  the  term 
has  been  extracted  in  front.  A      term  in  the  thermodynamic  potential  corresponds 
to  a      term  in  the  specific-heat.  Thus,  we  will  be  interested  only  in  extracting  the 
term  independent  of  temperature  from  the  square  bracket.  To  this  end  we  will  use  the 
Euler-Maclaurin  formula,  which  reads 
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where  are  the  Bernoulli's  coefficient  and  F*^''  is  the  2k  —  l"*  derivative  of  the 
function  F,  and  the  function  F  in  our  case  is 

,n2  +  X, 


F{7i)        +  4n{n  -  VnF+X)  +     log  | 


Prom  the  above  expression,  it  is  clear  that  the  only  term  which  can  contribute  to  a  cut  off 
and  temperature  independent  term  will  be  the     log        term.  We  express  the  series  as 

N 


2^{nf  login)   =   -2      x'^logxdx  -  N^logN 

n=l  •'I 


N 


/c=l 


(x^logx)t^^-^)(7V)  -  (a:Mog3r)(^^-i)(l) 
i2k)\ 
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thus  the  constant  term  is 

N 


AT 


and  we  obtain 


constant 


-2  ^     log  n 


-2  1  

^9     12  360 


C/m,2  7rr3C(3) 


03) 
27r2' 


(3.34) 
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The  remaining  f/^r^  contribution  to  the  specific  heat  comes  from  the  Ikp  region,  to 
evaluate  this  contribution  we  need  the  asymptotic  form  of  the  polarization  operator  given 
by 


where  x  —  {k  -  Ikp^jlkp.  The  square  of  the  second  term  in  the  bubble  yields  another 
r?  log  n  term,  which  as  before  is  responsible  for  the      term  in  the  specific  heat.  The 
square  of  the  second  term  in  11  yields, 


-o  y  n't/' 


iWn,k 


'^'^^       W  27r 


1  Um^ 


1  U'm 


2/3  (2 


dx 


X  +  \x\  + 


8\x\klvl\  ' 


(3.35) 


where  the  precise  upper  Umits  on  x  are  not  important,  the  log-divergence  of  a;  (at  x  =  0) 
is  cut  off  at  x  =  uJn/2kFVF-  We  obtain  a  similar  log  term. 


20 


2n 


vl  27r2' 


(3.36) 


Thus  both  k  <^  kjp  and  k  «  2fc/r  yield  identical  nonanalytic  contributions  to  the  specific 
heat.  The  total  contribution  is  twice  the  contribution  from  either  one  of  the  regions. 
3.3.3    Contact  Interaction:  Beyond  Second  Order 

In  this  Section  we  will  sum  up  the      contributions  from  all  the  terms  (t/";  n  >  2) 
of  the  thermodynamic  potential  Eq.  (3.17).  We  will  widen  the  scope  of  our  result  by 
including  self-energy  insertions  in  to  the  bare  Green's  function.  Although  the  general 
structure  of  such  a  Green's  function  is  very  complicated,  near  the  Fermi  surface,  one  can 
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Figure  3-4:  Plot  of  ul^  (u*)  vs  u* 

approximate  the  exact  Green's  function  by  its  expression  near  the  quasi-particle  pole, 

Z 


G{iuJn,k) 


where  Z  is  the  renormahzation  factor,  el  =  v*p{k  —  kp),  v*p  —  kf/m*,  and  m*  is  the 
renormahzed  mass.  Parameters  Z  and  m*  are  some  functions  of  the  bare  interaction  u, 
whose  forms,  in  general,  are  not  known.  This  amounts  to  replacing  the  prefactor  and  the 
Fermi  velocity  in  Eq.  (3.32) 

m*Z^ 


ll{ujm,k)     U*{um,k)  = 


2n 


1 


(3.37) 


As  noted  in  the  previous  Section  a      term  in  the  thermodynamic  potential  originates 
from  a  u;^  log  u  term,  which  is  in  turn  obtained  after  the  momentum  integration  of 
(-W^si=)2.     a  j^j  jjg^         i7"n*",  the  term  (  ,  j"'"'!    y  is  obtained  from  the  3'*'' 
term  of  the  binomial  expansion: 


^     27r  ' 


1-n- 


+ 


n{n  —  1) 


Using  this  result  in  conjugation  with  Eq.  (3.17),  we  obtain 

6C  (T)  /T 


+ 


(3.38) 


■^^2  C*  /E* 
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(3.39) 
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Here       =  7rm*r/3,  E*^  = 


fcf  and 


iu*f  [1  +  2u*  +  6{u*f  +  A{u*f  +  ...] 


-  1 


(3.40) 


where 


u*  =  Z'^m*U/2TT. 


(3.41) 


This  result  is  valid  for  0  <  «*  <  1.  The  divergence  of  ulf^  at  u*  =  1,  resulting  from  the 
spin-channel,  signals  an  instability  toward  a  magnetically-ordered  state.  A  plot  of  u^fj  (u*) 
is  presented  in  Fig.  3-4. 
3.3.4    Generic  Interaction 

In  Section  3.3.3,  we  obtained  the  specific  heat  of  a  contact  interaction  for  arbitrary 
interaction  strengths,  this  result  has  a  rather  limited  use  as  most  realistic  interactions 
have  a  finite  range.  To  obtain  a  form  that  is  valid  for  finite-range  potential  we  begin  by 
investigating  the  structure  of  vertices. 

Let  us  consider  higher-order  terms,  which  can  lead  to  two  kinds  of  corrections;  self- 
energy  insertions  and  vertex  corrections.  The  effect  of  self-energy  corrections  is  reflected  in 
a  renormahzed  effective  mass  and  a  Z-factor  for  the  quasi-particle.  The  second  correction 
is  attributed  to  vertex  terms.  A  typical  diagram  of  n"*  order  consists  of  n  bubbles.  To 
obtain  a     contribution  to  the  specific  heat  one  needs  a  product  of  two  \0.m\/Q  terms 
which  are  obtained  ft-om  any  two  out  of  the  n  bubbles.  Any  extra  factor  of  \Qm\/Q 
will  not  lead  to  a  nonanalytic  term,  thus  we  can  set       =  0  and  Q  — »  0  in  the  rest 
of  the  diagram.  Thus,  even  at  arbitrary  order  we  are  dealing  with  only  two  bubbles. 
The  effect  of  higher  order  terms  can  be  hidden  in  the  static  vertices  (these  are  static 
because  we  need  to  first  set       =  0  and  then  Q  ->  0)  joining  the  two  bubbles.  These 
vertices  are  of  two  types:  r*(fc,  -A;;  k,  -k)  and  r*^(A;,  -k;  -k,  k).  What  we  are  left  with 
is  a  "2"'*"  order  diagram,  where  the  wavy  lines  are  replaced  by  two  blocks  represented 
either  by  a  r'=(A;,  -k;  k,  -k)  or  a  r*(A;,  -k;  -k,  k).  In  the  conventional  notation  the  vertices 
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T''{k,  -k;  k,  -k)  and  r*'(A;,  -A;;  -A;,  A;)  are  related  to  r'^(^)  via, 

t''{e)=    lim   r{kFhuO;kFn2,0\Q,n),  (3.42) 

|12|/Q  ^0 

where  r{ki,u!i;k2,uj2\Q,Q.)  is  the  vertex  for  a  process  (ki,ui;k2,i^2  — +  fci  —  Q,^^l  - 
fi;  k2  +  Q,uj2  +     and  6  is  the  angle  the  two  incoming  quasi-particles.  Since  the  incoming 
momenta  are  nearly  anti-parallel  to  each  other,  we  can  set  the  angle  to  be  equal  to  ^  =  tt. 
The  vertex  with  the  full  spin  representation  is  then, 

^U,-rsi^)  =        -k;  k,  -k)5aj50s  -  T^{k,  -k;  -k,  k)SasS0y.  (3.43) 

We  can  express  the  scattering  amplitude  in  terms  of  the  vertex  function  via  the  relation, 

A.0,,s{Tr)  -  Z'Ti^^^sin),  (3.44) 

where  the  scattering  amphtude  itself  is  represented  in  terms  of  spin  and  charge  compo- 
nents, 


7r 


rrr 


n 


=   :;^[iMT^)-MT^))Sa^S0s  +  2As{Tr)6^5S0y].  (3.45) 

Prom  Eqs.  (3.43),  (3.44)  and  (3.45)  we  obtain  the  relationship  between  the  spin  and 
charge  components  and  the  vertex  function. 


Z'r'{k,-~k;k,-k)   =        [A,{Tr)  -  AsiTT 


vrr 


ZV{k,-k-k,-k)   =  -2—Asi-n). 


(3.46) 


Finally  we  will  substitute  Z'^T>'{k,  -k;  k,  -k)  for  C/(0)  and  Z^r>'{k,  -k;  -k,  k)  for  U{2kF) 
in  Eq.  (3.14),  to  obtain  the  nonanalytic  form  of  the  specific  heat  in  a  generic  Fermi 
liquid  [85], 

{Z^r\k,  -k;  k,  -k)f  +  (Z^T^ik,  -k;  -k,  k)f 


5C{T) 

9C(3) 

"  m  ■ 
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-{Z^V\k,  -k;  k,  -k)){Z^r'{k,  -k;  -k,  k)) 
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(3.47) 
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Notice  the  universal  sub-leading  term  in  the  specific  heat  is  not  an  angular  average 
over  all  the  angles  and  is  different  from  the  leading,  Fermi-liquid  specific  heat  term 
{C{T)/T  oc  constant)  which  is  expressed  as  an  "angular-average- ((Fc(^)  cos  ^))"  over  the 
Fermi-liquid  interaction  function. 

The  result  for  Cv  has  been  expressed  in  terms  of  the  scattering  amphtude.  Alterna- 
tively, the  specific  heat  can  be  expressed  in  terms  of  the  harmonics  of  the  quasi-particle 
interaction  function  F{0).  The  harmonics  of  the  scattering  amplitude  and  the  interaction 
function  are  related  via 


In  contrast  to  the  T  term  in  the  specific  heat  which  requires  only  one  harmonic,  F^,  the 

term  requires  the  full  set  of  harmonics.  The  parameters  of  E)q.  (3.47)  can  be  extracted 
from  the  specific  heat  measurement  on  the  fluid  monolayer  of  ^He  [27].  To  a  reasonable 
accuracy,  the  data  can  be  fitted  by  a  form  C/{NTE*p)  =  -f{T/E*p),  where  N  is  the  density 
per  unit  area  in  a  fluid  monolayer,  Ep  =  EF{m/m*)  and  7(x)  ^  a-bx.  A  fit  to  the  data 
then  yields  b  =  0.9{Al{n)  +  3Al[n]]  «  9  -  14,  thus  [A^^{tt)  +  SAUtt]]  «  11  -  15.5.  In  the 
case  of  "almost  locaUzed  fermions"  [86],  the  charge  channel  is  the  dominant  interaction 
and  Fi"^  >  1.  Thus  the  scattering  amplitude  term,  Ac{n),  is  expected  to  be  small  due  to 
cancellation  between  the  subsequent  terms  in  the  expansion.  Neglecting  Ac{ti),  the  spin 
part  of  the  amplitude  is  obtained  to  be  \As{7t)\      1.9  -  2.3.  In  addition,  if  the  n  =  0 
harmonic  of  Fg  is  the  leading  contributor  to  ^^(Tr),  then  Fi°^  «  -(0.66  -  0.7).  This  is 
consistent  with  the  3D  value  F^^  ^  -0.75  [25]. 


We  have  shown  earlier  that  the  collective  mode  for  a  finite  range  potential  leads  to  a 
nonanalytic  behavior  of  the  self-energy  near  the  mass  shell.  However,  this  nonanalyticity 
is  not  manifested  in  the  contribution  to  the      term  in  the  specific  heat.  For  a  long 


c/s 


thus 


3.4    Specific  Heat  for  the  Coulomb  Potential 
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range  Coulomb  potential,  we  obtain  a     correction  to  the  specific  heat  both  from  the 
particle-hole  and  from  the  plasmon  regime.  In  the  temperature  range  T  <^  nvp  <SC  Ep, 
where  k  =  me^  is  the  screening  wavenumber,  the  correction  is  independent  of  the 
interaction  strength.  For  the  Coulomb  potential,  the  singular  nature  of  the  interaction 
for  small  momentum  transfer  allows  one  to  consider  an  effective  RPA  potential  for  the 
self-energy,  and  for  the  thermodynamic  potential  we  need  to  consider  the  usual  series  of 
ring  diagrams.  A  knowledge  of  the  thermodynamic  potential  allows  us  to  evaluate  the 
entropy  and  the  specific  heat.  The  finite  temperature  thermodynamic  potential  in  the  ring 
diagram  series  is  given  by 

T   =  7^1og(l-^Wnm(fc.i^n)),  (3.48) 

iWn  ^  ' 

where  the  summation  is  over  the  bosonic  frequencies,  Wm  is  the  polarization  operator  (we 
include  contributions  from  both  up  and  down-spin  electrons),  and  V'(fc)  =  the  2D 

Coulomb  potential. 

Using  Eq.  (3.21)  to  convert  the  frequency  summation  into  frequency  integrals,  we 
obtain 

_i  cdu^     ^         ,amD^{kM]  ..... 

2  y  (27r)2  7   TT  sinh2(a;/2T)         lReD'i{k,co)\'  ^^'^^^ 

where  Z)^  =  1  -  V{k)Um{k,iWn)-  The  regions  contributing  to  the  entropy  are  the  plasmon 
region  and  the  particle-hole  region  (Fig.  3-5).  Near  the  plasmon  region, 

g-'^i-vMnV^j^'"^'"?-"'.  (3.50) 

thus  ImJ?^  ^  0+  and  (1  -  V{k)ReIl^)  <  0  giving  a  value  of  n  for  the  argument  in  the 
integrand.  The  contribution  from  the  plasmon  region  is  thus 

CI  _       r^dw  ^  kdk 

Jo  smh^{w/2T)  JjuL  2tt 
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co<k) 


A  B 

Figure  3-5:  Contour  for  summing  up  the  discrete  sum  of  the  Matsubara  thermodynamical 
potential.  A)  Matsubara  poles  at  the  Matsubara  frequencies  and  at  the  plasmon  positions 
—LOpi  and  LOpi  (for  the  integration  variable  gg  );  the  branch  cut  from  — u-'fc  to  uJk  represents 
the  particle-hole  region.  B)  Particle-hole  spectrum  for  2D. 

The  other  contribution  to  the  entropy  is  from  the  particle-hole  region.  The  argument 
[Eqn.  3.49]  decays  as  one  goes  deep  into  the  particle-hole  region.  In  the  particle-hole 
region  the  polarization  bubble  is  given  by 


n>,g)  =  --(l+7  ^ 


(3.51) 


thus, 


(3.52) 


\particle-hole  ~  I  7;; — TIT  /  .  ,  9  ,  7r;rr  tan 

J  (27r)^  io  sinli^(u;/2r) 


2me^u'  /  y/k-  —  w- 


k  +  '2me^ 


,  (3.53) 


We  will  make  the  substitutions  w/2T  ->  z  and  VFk/2T  k',  and  also  a  change  in  the 
order  of  integration  to  obtain 


S\ 


particle— hole 


J_2r  dzz  /^,^,^Ytan- 


),  (3.54) 


l2Tk'/iirs  +  l- 

where  r«  =  n/kp  and  //  =  A;fUf/2.  It  will  be  convenient  to  spht  the  integrals  into  two 
parts,  the  first  part  gives  the  hnear  in  T  correction  to  the  entropy  and  a      term,  whereas 
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the  second  integral  gives      correction  to  the  entropy. 

^Iparticie-noie   -      2n^v,Uo    7^ sinh\z)  J ,        Tk' /{fir,)  +  1 
2n  Vf    Jq  7rsinh^(2) 

the  expression  in  the  big  bracket  is  convergent  and  the  main  contribution  comes  from 
A;'  ~  1,  hence  for  T  <§:  kvf  we  can  make  a  further  simphfication, 

1  ,2T.o  dzz 


'  vf    Jo  ttsu 


Jo  7rsinh^(z) 

X 


Finally  we  have 

_   r,log(r,)r      1   2r  23C(3)      1  .2r        3C(3)  3C(3)x 

=  +  (3.56) 

3m  Vf  iir 

So  the  linear  in  T  contribution  to  the  specific  heat  from  the  particle-hole  region  is 
C  =  rs\og{rs)T/{2,m).  The  7^  term  due  to  the  particle-hole  term  is 

S\partide-hole     =  ~47r^'  (3.57) 

Adding  the  contributions  from  the  plasmon  and  the  particle-hole  region,  we  get  for  the 
nonanalytic  correction  to  the  entropy 

S\parUde-hole+pla3mon    ~    ~  (^^^  ~4tt    '  (3.58) 

Hence  the      contribution  to  the  specific  heat  C  =  TdS/dT\^,A  is 


We  notice  that  in  this  temperature  range  the  specific  heat  is  independent  of  the  inter- 
action strength.  A  further  discussion  is  provided  in  Appendix  A  where  we  reconsider 
the  plasmon-pole  and  particle-hole  contribution  to  the  specific  heat  from  an  alternate 
approach.  We  obtain  the  same  results  as  before 

3.5  Conclusion 

In  this  Chapter,  the  important  question  regarding  the  role  of  nonperturbative 
interaction  effects  in  the  specific  heat  has  been  considered.  Often  while  evaluating  the 
contributions  to  the  specific  heat  firom  the  self-energy,  the  contribution  from  the  imaginary 
part  of  the  self-energy  is  neglected.  We  have  found  that  for  a  short  range  potential,  the 
imaginary  part  of  the  self-energy  emanating  from  the  zero-sound  mode  region,  ImE^s* 
does  contribute  to  the  nonanalytic  term  in  the  self-energy.  In  fact,  this  contribution 
exactly  cancels  the  contribution  coming  from  the  real  part  of  the  self-energy  (whose 
contribution  also  originates  from  the  zero-sound  region).  This  result  implies  that  a 
nonanalytic  correction  to  the  self-energy  arises  only  due  to  perturbative  processes.  We 
have  verified  this  result  directly  by  calculating  the  specific  heat  from  the  thermodynamic 
potential.  To  this  end  use  has  been  made  of  both  the  real  and  Matsubara  frequency 
formalism.  An  advantage  of  using  the  Matsubara  formahsm  is  the  ease  in  identifying 
the  terms  which  gives  rise  to  the  nonanalytic  contributions  to  the  specific  heat.  This  has 
helped  us  in  extending  our  second  order  result  to  all  orders.  This  extension  to  all  orders  is 
required  for  cases  when  the  interaction  strengths  are  not  weak.  In  a  further  extension  of 
this  work,  a  finite-range  potential  was  considered.  We  have  been  able  to  obtain  a  closed 
form  solution  for  the  contributions  to  the  specific  heat  to  this  problem  as  well.  Finally,  we 
have  considered  the  case  of  long  range  Coulomb  potential  and  shown  that  perturbative 
process  do  contribute  to  the  specific  heat  and,  in  the  low-T  limit,  this  contribution  is 
interaction-independent. 


CHAPTER  4 
ANOMALOUS  EFFECTIVE  MASS 

The  observation  of  a  metallic-like  resistivity  and  an  apparent  metal-insulator  tran- 
sition in  high-mobility  Si-MOSFETs  in  1994  [50,  51]  challenged  the  scahng  theory  of 
locahzation  and  has  since  led  a  number  of  groups  to  study  both  the  transport  and  ther- 
modynamic properties  of  Si-MOSFET's  and  other  semi-conductor  heterostructures.  The 
origin  of  the  anomalous  metallic  behavior  still  remains  a  subject  of  active  research.  As 
discussed  in  Section  1.3,  although  models  based  on  the  conventional  dirty  Fermi-liquid  can 
account  for  many  observed  effects  both  qualitatively  and  quantitatively,  there  are  quite 
a  few  proposals  for  the  non-Fermi  liquid  origin  of  the  anomalous  metalhc  state.  On  the 
experimental  side,  the  main  qualitative  argument  for  the  Fermi-liquid  nature  of  the  metal- 
lic state  is  the  observation  of  conventional  Shubnikov-de  Haas  oscillations,  which  implies 
an  existence  of  well  defined  quasi-particles.  The  parameters  of  the  2DES,  the  eflFective 
mass  m*  and  the  Lande  p-factor  are  obtained  by  fitting  the  Lifshitz-Kosevich  formula 
to  the  Shubnikov-de  Haas  (SdH)  oscillation  data.  Excellent  fit  to  the  SdH  data  using  -  . 
the  renormahzed  parameters  and  the  matching  of  these  parameters  with  the  parameters 
(X*  oc  m*g*)  obtained  via  other  experimental  techniques  Uke  saturation  of  magnetoresis- 
tance  approach  [64,  66,  67],  give  strong  evidence  for  the  Fermi-liquid  behavior  of  the  2D 
system.  In  particular,  studies  indicate  the  spin  susceptibility  x*  shows  a  strong  increase 
as  the  density  approaches  the  critical  one  (insulating  regime)  with  a  large  fraction  of  this 
increase  attributed  to  the  increase  in  the  effective  mass  term  than  g',  an  attribute  very 
similar  to  the  He^  case. 

The  magnetic  field  dependence  of  the  effective  mass  serves  as  an  additional  check  for 

Fermi-liquid  behavior.  Experimental  observation  of  such  behavior  in  2D  Si-MOSFET's 

in  the  presence  of  a  parallel  magnetic  field  shows  that,  contrary  to  the  expectations, 

the  effective  mass  shows  no  dependence  on  the  degree  of  polarization  [68]  (Fig.  4-1). 

Although  the  effective  mass  exhibits  strong  dependence  on  the  density  and  hence  the 
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Figure  4-1:  Effective  mass  vs  the  degree  of  spin-polarization.  The  electron  densities  in 
units  of  10"cm~^  are:  1.32  (dots),  1.47  (squares),  2.07(diamonds),  and  2.67  (triangles). 

Fermi  energy  Ep,  the  lack  of  dependence  on  the  magnetic-field,  which  alters  the  Fermi 
energy,  is  quite  surprising  as  it  does  not  follow  the  prediction  of  Fermi-liquid  theory.  That 
the  masses  should  depend  on  polarization  can  be  seen  from  considering  two  limiting  cases: 
of  zero-  and  100%  polarization.  At  fixed  total  density,  the  Fermi  energy  is  doubled  by  fully 
polarizing  the  2D  system.  Thus  the  ratio  of  Coulomb  to  Fermi  energy  g  e^^/im/EF, 
which  is  a  measure  of  interaction  strength  is  halved  compared  to  the  unpolarized  case. 
Yet  the  effective  mass  remains  polarization  independent.  In  this  Chapter  we  analyze  this 
paradox  in  detail. 

The  present  Chapter  is  divided  into  two  parts.  Section  4.1  deals  with  the  polariza- 
tion dependence  of  the  effective  mass.  We  will  consider  two  independent  approaches, 
Landau's  phenomenological  approach  in  Section  4.1.1  and  the  weak  coupling  approach  in 
Section  4.1.2,  to  show  that  in  a  Fermi  liquid  the  effective  mass  not  only  acquires  a  field  de- 
pendence but  also  exhibits  spin-splitting.  In  the  second  part  of  this  Chapter,  Section  4.2.1 
and  Section  4.2.2,  we  examine  the  influence  of  valley  degeneracy  on  the  effective  mass 
and  on  the  spin  susceptibility,  respectively  The  importance  of  valley  degeneracy  on  the 
transport  properties  in  the  (001)  plane  of  Si-MOSFET  has  been  emphasized  in  a  recent 
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study  by  Punnoose  and  Finkelstein  [61].  The  diffuson  propagator  in  a  spin- 1/2  system 
has  4  channels,  1  singlet  and  3  triplet,  whereas  for  a  system  with  spin-1/2  and  2-vaUey  de- 
generacy (total  degeneracy  is  4),  as  in  the  (001)  plane  of  Si-MOSFET,  there  are  1  singlet 
and  15  triplet  channels.  These  additional  channels  have  been  held  responsible  for  a  large 
(almost  an  order  of  magnitude)  drop  in  resistivity  near  the  critical  region.  Motivated  by 
these  findings,  we  consider  a  Coulomb  gas  in  the  large- A''  hmit  (for  Si-MOSFETs,  A''  =  4). 
In  a  system  with  no  valley  degeneracy,  the  effective  mass  is  renormalized  by  emission 
of  both  virtual  electron-hole  pairs  and  plasmons,  in  the  large  N-limit  of  the  Coulomb 
gas,  the  effective  mass  is  renormalized  primarily  due  to  a  polaronic  effect:  via  emission 
and  absorpbtion  of  (virtual)  high  energy  plasmons,  while  electron-hole  pairs  play  only  a 
nominal  role.  As  plasmons  are  classical  objects,  the  quantum  degeneracy  and,  hence  the 
polarization,  does  not  affect  the  effective  mass  to  the  leading  order  in         The  large- 
expansion  obtained  is  rapidly  convergent  even  for  non-valley  degenerate  system  (A''  =  2) 
and,  as  such,  it  provides  a  non-trivial  way  of  going  beyond  the  weak  coupling  limit. 

In  Section  4.2.2  we  investigate  the  behavior  of  the  spin  susceptibility  in  the  large- A'^ 
limit.  We  find  that  the  spin  susceptibility  is  renormalized;  however  this  renormalization 
in  leading  and  subleading  order  is  due  to  the  renormalized  effective  mass,  whereas  the 
Lande-g  factor  remains  unrenormalized.  This  indicates  that  the  enhancement  of  the  spin 
susceptibility  does  not  reflect  the  tendency  towards  ferromagnetic  ordering. 

4.1    Spin-Polarized  Effective  Mass 

In  the  following  two  Sections  we  show  how  a  standard  treatment  results  in  a  polar- 
ization dependent  effective  mass  in  the  presence  of  a  magnetic  field  applied  parallel  to  the 
2D  surface.  The  parallel  magnetic  field  is  assumed  to  couple  only  to  the  spin  degree  of 
freedom  and  not  to  the  orbital  degree  of  freedom,  and  acts  as  a  knob  which  can  be  used 
to  change  the  population  of  either  kind  of  electrons,  while  keeping  the  total  density  of 
electrons  fixed. 

4.1.1    Landau's  Phenomenological  Approach 

First,  we  derive  an  effective  mass  expression  for  a  partially  spin-polarized  Fermi  liquid 
system,  generahzing  the  approach  used  for  the  single  component  case  [87].  The  form  of  the 
Landau  interaction  function  for  a  spin-polarized  case  is,  however,  much  more  comphcated 
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due  to  the  breakdown  of  SU(2)  symmetry.  We  will  show  that  this  interaction  function 
effectively  reduces  to  three  independent  terms  which  describe  the  interaction  between  two 
up-spin  electrons,  between  two  down-spin  electrons,  and  between  an  up-spin  electron  and 
a  down-spin  electron.  The  strategy  to  obtain  the  mass  is  to  compare  the  change  in  the 
quasi-particle  energy  between  a  moving  frame  and  a  fixed  frame  via  two  different  ways. 
In  the  first  way,  Galilean  invariance  is  utilized  to  find  the  difference  in  the  energy  of  a 
quasi-particle  between  the  rest  frame  and  a  moving  frame.  In  the  second  way,  one  notices 
that  the  distribution  functions  are  different  in  different  frames,  and  thus  we  can  calculate 
the  change  in  the  quasi-particle  energy  which  arises  due  to  a  variation  in  the  distribution 
function. 

The  general  form  of  the  interaction  function  in  a  spin-polarized  system  is  given 
by  [88] 

/a,;i;/3,i/(P,  P ')    =    F{p,  p  ')6c,g6^^u  +  G{p,  p  ')aa,0.a^^^  '  +  Z.  [X{p,  p  ')aa,(3Sf,,^ 

+X{p',p)d^^JaA  +  yip,P')[doc,0-i][d^,u.i\,  (4.1) 

where  the  field  is  in  the  z  direction,  and  a  are  the  Pauli  matrices. 

The  change  in  quasi-particle  energy  due  to  a  change  in  the  quasi-particle  distribution 
function  is  given  by 

5e{p)  ^Tr'l f{p,a;p',a')Sn{p')0^^,  (4.2) 
and  the  distribution  matrix  5n  is 

y   0     dni  j  ' 

where  the  change  in  distribution  function  for  the  majority  (spin-up)  quasi-particle  is  given 
by  Sn-^  and  for  the  minority  (spin-down)  quasi-particle  by  Sni-  We  explicitly  calculate  the 
product  of  the  density  matrix  Jn^,^  with  the  interaction  function,  to  obtain  the  following 
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terms: 


7.X{p,p')da,05^,Mn,u  =  X{p,p')(ylA5n^  +  5ni\ 


and 


zX{p',p)d^^Jn^,^Ja,0   =   X{p',p)alJnf,Jc,0  =  X{p',p)[6n^-Sni%,0 
Yip,p')al0alJn,,.   =   Y{p,p ')al0[6n^  -  5n^].  (4.3) 

It  is  easy  to  deduce  that  the  quasi-particle  energy  matrix  is  in  a  diagonal  form  since  6a,i3, 
0  and  Sna^g  are  all  in  a  diagonal  form.  Thus  for  a  spin-up  quasi-particle  we  obtain  the 
following  relation 

Se^ip)   =  1 0^  Sn^{p'){F{p,p')  +  G{p,p')+X{p,p')  +  X{p',p)  +  Y{p,p')] 

+5n^ip')[F{p,p')-G{p,p')+X{p,p')-X{p',p)-Yip,p')] 
Similarly,  for  a  spin-down  quasi-particle 

Ssiip)   ^  J  0^  6n^{p'){F{p,p')-G{p,p')-Xip,p')+X{p\p}-Y{p,p')] 

+5n^{p')[Fip,p')  +  G{p,p')  -  X{p,p')  -  X{p',p)  +  y(p,p ')} 

To  simplify  the  algebra  we  will  define  a  new  interaction  function  F*^{p,p'),  where 
hj  =T)  i  such  that 

cPp' 


Se^ip)   =  J 


(27r)2 


F^'^  {p,  p  ')6n^  {p  ')  +  F^'^  {p,  p  ')Sni  {p  ') 
F^'T  {p,  p  ')Sn^  {p')  +  F^'^  {p,  p  ')5n^ {p  ') 
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where 


FT'T(p,p')  =  F{p,p')  +  G{p,p')  +  X{p,p')+X{p',p)  +  Y{p,p') 

F^'^{p,p')  =  F{p,p')-G{p,p')  +  X{p,p')-Xip',p)-Y{p,p') 

F^'^{p,p')  =  Fip,p')-G{p,p')-Xip,p')  +  X{p',p)-Y{p,p') 

F^^^{p,p')  =  F{p,p')  +  G{p,p')-X{p,p')-X{p',p)  +  Y{p,p'). 

Thus,  we  can  invoke  GaUlean  invariance  to  obtain  the  effective  mass  [87].  Using  Galilean 
invariance  we  can  relate  the  quasi-particle  energy  in  a  moving  frame  e'  with  the  quasi- 
particle  energy  in  a  laboratory  frame  e.  Let  us  assume  the  moving-frame  has  a  velocity  u 
with  respect  to  the  lab  frame.  Introducing  a  quasi-particle  with  momentum  p  in  the  lab 
frame,  we  increase  the  total  mass  of  the  system  by  the  bare  mass  m,  since  a  new  particle 
has  been  added.  The  quasi-particle  energy  in  the  moving  frame  is 

2 

e'ip-mu)   =  e{p)-pu-{ — — 


p'  =  p  -  mu, 

where  p  (p  ')  is  the  quasi-particle  momentum  in  the  lab  frame  (moving  frame)  and  m  is 
the  bare  mass  of  the  particle.  We  thus  obtain 

e'ip)  =  e{p  +  mu)  -  {p  +  mu)  ■  u  H — — . 

The  right  hand  side  of  the  equation  can  be  expanded  to  linear  order  in  u  to  obtain  the 
relation 


t,-^         /-A  (m-m*)p-u 


m* 

Hence  for  the  particle  of  the  first  kind  (spin-up)  we  have 


m*^ 


(4.4) 


where  the  bare  mass  for  spin-up  and  spin-down  quasi-particles  is  the  same  and  equals 
m.  A  second  relation  is  obtained  by  observing  that  the  quasi-paricle  energy  e'ip)  in  the 
moving  frame  differs  from  the  quasi-paricle  energy  eijp)  in  the  lab  frame  because  of  a 
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n',(F) 


Figure  4-2:  Distribution  function  in  the  primed  frame  and  lab  frame  are  represented  by 
dashed  and  full  circle  respectively. 

different  distribution  function.  The  shift  in  the  energy  spectrum  5e{p)  =  e'(p)  -  e{p)  due 
to  a  different  distribution  [the  new  distribution  is  shifted  by  —mu  in  momentum  space  for 
both  up-  and  down-spin  particles,  (Fig.  4-2)]  is 

(fp' 


Hp) 


=  tr'y"  fip,p')6hp. 


(4.5) 


where  tr'  is  the  trace.  For  spin-up 

(Pp' 


and  for  spin  down 


F^'^  {p,  p  ')Sn^  ip')  +  F^'^  {p,  p  ')5ni  {p  ') 


(4.6) 


where  we  have  used 


»  ^  / 


cPp' 


F^'\p,p  ')5n0 ')  +  F^'^ip, p  ')5n^{p  ') 


(4.7) 


5n^{p')  =  n\{p')-n]{p'). 


(4.8) 


There  is  a  very  simple  relationship  between  the  distribution  function  in  the  primed  frame 
and  that  in  the  lab  frame:  a  particle  moving  with  momentum  p  '  in  the  primed  frame  will 
have  momentum  p  '  +  mu  in  the  rest  frame,  thus  the  distribution  functions  in  the  different 
frames  will  obey  the  relation  n\{p')  =  n^(p  '  +  mu),  where  n\  (n^)  is  the  distribution  in 
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the  primed  (rest)  frame.  The  two  distribution  functions  axe  shown  in  Fig.  4-2.  Continuing 
with  Eq.  (4.8)  we  obtain 

n'T(p')-n?(p')   =  n'^ip'  +  mu)-n^p') 

p'dn'^ip') 


=  mu  ■ 


m*  de^ip') 
mu  ■  p  ' 


m^ 


«5(et(p?t)-ef(p')),  (4.9) 


where  pp^^  is  the  Fermi  wave  vector  of  the  up-spin  electron.  Similarly  for  the  spin-down 
particle,  the  change  in  distribution  is 

.  ^  p'dn\{p')        mu-p'..    _  . 

6n,ip  )  =  mu.  =  --—5{e,{pn)  -  e,{p  ')), 

where  p>i  is  the  Fermi  wave  vector  of  the  down-spin  electron.  For  a  spin-polarized  system 
the  Fermi  wave  vector  of  the  up-spin  electron,  pp],  will  be  greater  than  that  of  the  down- 
spin  electron,  pf[.  Substituting  the  simplified  form  of  5n-^{p')  and  5n^{p')  in  Eq.  (4.6) 
and  Eq.  (4.7),  and  using  the  relations 

<fp  m*  dedO 

(27r)2  " 

/-aTr 


/  de'F{p-p')u-p'  =  /  de'F{\p'\\p\cose')\u\\p'\cos{e-e') 

Jo  Jo 

=   |w||p'|cos(^)  /    de'F{\p'\\p\  COS  9')  COS  9',  (4.10) 
Jo 


we  obtain 


Se-t  =  - 


m  _ 
27r 


m  _ 


/dO' 
-^cos9'F^'^PF^,pj.^) 

-u-PFi  J  ^  cose' F^^^{Pfj,Pf^  ') 

^  =   -^"-my  ^cosm^(pfi,p>') 

m  ^    ^    f  d9'       ,   ,  , 
-^^■PFl  J  ^cos9'F^'^{pFi,pFi'),  (4.11) 


where  p>|  and  p>|  point  in  the  same  direction,  and  in  general  -F^'^(pft,Pft  ') 
F^'^(Pfi,Pfi  ')  7^  F^'\pFl,pF^  ')•  Comparing  Eqs.  (4.4)  and  (4.11)  we  obtain  for  the 
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up-spin  fermions 


It  is  clear  from  the  above  equations  that  the  effective  mass  for  the  up  and  down-spin 
electrons  is,  generally,  diflferent.  An  explicit  evaluation  of  the  effective  mass  requires  a 
knowledge  of  the  F  function  which,  however,  is  not  provided  within  Landau's  phenomeno- 
logicaJ  theory.  In  the  following  Section,  we  will  perturbatively  evaluate  the  effective  mass 
of  a  system  of  electrons  interacting  with  the  long  range  Coulomb  interaction  and  confirm 
our  prediction  regarding  the  spin  sphtting  of  the  effective  mass. 
4.1.2   Weak-Coupling  Approach 

A  knowledge  of  the  exact  eigen  states  of  a  Hamiltonian  allows  us  to  calculate  relevant 
physical  quantities.  In  many  cases,  though,  it  is  simply  too  difficult  to  calculate  these 
eigen  states.  In  these  situations  one  resorts  to  perturbative  techniques  and  obtains 
the  physical  quantities  as  a  perturbative  expansion  in  terms  of  a  small  parameter  of 
the  system.  We  resort  to  this  approach  to  calculate  the  effective  mass  in  terms  of  the 
dimensionless  constant      =  ^/2me^/kp  (a  measure  of  the  interaction  strength).  The 
stronger  the  interaction,  the  more  the  deviation  of  the  electron  effective  mass  from  the 
bare  electron  mass.  For  a  system  with  different  numbers  of  up-  and  down-spin  electrons, 
there  will  be  unequal  remormalization  of  the  masses  and  hence  spin-split  masses. 

The  Hamiltonian  of  our  system  is  given  by 


where  the  summation  on  repeated  spin  indices  is  implied  and  U{r  -  r')  is  the  usual 
Coulomb  interaction  potential.  The  Fermi  wave  vectors  for  spin-up  and  spin-down 
electrons  are  different:  k]}  =  kFy/T±J,  where    =  {n^-n^)/n  is  the  degree  of  polarization. 


similarly  for  the  down-spin  fermions 


H  = 


(4.12) 
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Figure  4-3:  Self-Energy  of  spin-up  and  spin-down  electrons. 


We  wiU  consider  all  interactions  present  in  this  system:  between  two  up-spin  electrons, 
between  two  down-spin  electrons  and  between  an  up-spin  and  a  down-spin  electron. 

The  electrons  interact  via  a  Coulomb  potential.  This  potential  is  screened  due  to  the 
creation  of  electron-hole  pairs.  To  calculate  the  self-energy  (required  to  obtain  the  effective 
mass),  we  choose  those  diagrams  which  at  each  order  in  the  interaction  have  maximum 
divergence.  Self-energy  diagrams  given  in  Fig.  4-3  contain  all  such  contributions.  The 
series  when  summed  to  all  orders  [the  choice  of  this  particular  series  is  known  as  the 
Random  Phase  Approximation  (RPA)]  yields  an  effective  potential  which  is  dynamically 
screened. 

The  effective  mass  in  terms  of  the  self-energy  is  defined  as 

1 


m 


Oil' 


m 


*n 


1  - 


diw 


(4.13) 


where  E^^  is  the  self-energy  of  up  and  down-spin  electrons  with  momentum  p,  frequency 
iW  and  quasi-particle  energy  ^p.  For  p  ~  A;Ji^,  ^Ti  =  /jU^^  -  v]}),  where      and  Vp  are 
the  Fermi  momenta  and  velocities  respectively.  Let  us  consider  the  self-energy  diagram  for 
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spin-up  electrons  given  in  Fig.  4-3.  Because  the  Fermi  momenta  of  spin  up  and  spin  down 
electrons  are  different,  a  probe  electron  causes  different  variations  in  the  densities  of  the 
two  spin  orientations.  The  resulting  potential  contains  the  sum  of  two  contributions  from 
both  spin  orientations.  Summing  up  the  RPA  series,  we  obtain  the  following  expression  for 
the  self-energy: 

_  r+°°dw        kdk  V^{k)Y[{k,iw) 
J-oo  STrio    {27ryi -V{k)Uik,iw) 

I    d9G\Up  +  k,  iw  +  iw),  (4.14) 


X 


where  V{k)  =  2-Ke^/k  is  the  bare  potential,  G  is  the  free  single-particle  Green's  function 
and  n  is  the  polarization  operator.  The  momentum-frequency  integrals  in  Eq.  (4.14)  are 
dominated  by  the  regions  k  ^  k      kp,  where  k  ~  rgkp  is  the  screening  length  and 
uj  ~  kvp.  Therefore  the  long- wavelength  limit  of  the  polarization  operator  n(A;,u;)  can  be 
used.  This  form  is 

n  =  n+  +  n-,  n^Uw^k)  =  -^(i  -   ). 

^      '    '  27rV  ^(fct,±)2+u,2'' 

At  this  point  we  will  summarize  the  approach  that  has  been  taken  to  calculate  the 
effective  mass.  In  effect,  we  need  to  evaluate  two  quantities,  a  self-energy  differentiated 
with  respect  to  the  external  frequency  and  a  self-energy  differentiated  with  respect  to 
We  find  that  when  differentiating  Eq.  (4.14)  with  respect  to  Cfc,  the  two  individual  terms 
diverge,  however  the  divergences  cancel  each  other.  The  first  term  of  the  self-energy  yields 
a  log  singularity  which  is  exactly  cancelled  by  a  similar  log  singularity  arising  from  the 
second  term.  Similar  care  needs  to  be  taken  when  differentiating  Eq.  (4.14)  with  respect 
to  the  fi-equency.  Details  of  the  calculation  are  provided  in  Appendix  B.  The  result 
essentially  states  that  the  effective  mass  of  spin-up  and  spin-down  electrons  depends  on 
the  interaction  parameter     and  the  degree  of  polarization, 

/I     .     1  X    .  .Tl 
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Figure  4-4:  EflFective  mass  for  a  spin-up  and  a  spin-down  electron  at  polarization  ^  =  0.3. 
hence  expanding  to  the  lowest  order  in  polarization  we  obtain 

The  general  behavior  of  the  effective  mass  for  ^  =  0,  even  the  decrease  of  the  effective 
mass  compared  to  the  bare  mass  at  high  density,  is  in  agreement  with  a  recent  experiment 
on  the  2D  electron  system  in  an  GaAs/AlGaAs  heterostructure  [89].  On  the  other  hand 
the  above  theoretical  predictions  regarding  the  presence  of  a  spin-split  effective  mass 
(Fig.  4-4)  is  not  consistent  with  the  experimental  findings.  The  LK  formula  for  the  case  of 
magnetic  field  tilted  with  respect  to  the  2D  electron  gas  plane  reads 

Pxi/po   =   ^AkCos[2n^khn„/eBj_- n],  (4.17) 

where  a  =],  [  is  the  spin  index,  and  Ak  is  the  ampHtude  of  the  k^^  harmonic  which  is  a 
function  of  the  electron  mass  and  is  assumed  to  be  the  same  for  spin-up  and  spin-down 
electrons.  If  however,  the  masses  were  different  then  the  corresponding  amphtudes  would 
differ  and  the  LK  formula  would  change  to 


—  =  y"^(m^)cos[27r^s/in<,/eBx  -  tt] 

+  ^{Mm^)  -  Asim^))  cos[^^^-^  -  tt].  (4.18) 


The  frequency  of  the  oscillation  is  dependent  on  the  density  of  the  components  and  for  a 
spin-polarized  system  a  mismatch  in  frequency  results  in  a  beat  pattern  to  be  observed  in 
the  experimental  data.  Due  to  the  difference  in  amplitude  term  {AAs),  the  beats  will  not 
be  complete.  Nevertheless  the  experimental  data  show  full  beats,  implying  that  AAg  =  0 
and  that  the  effective  masses  for  the  up-spin  and  down-spin  electrons  are  equal  [66].  On 
a  similar  note,  the  experimental  group  of  Shashkin  et  al.  [68]  found  no  dependence  of  the 
effective  mass  on  the  degree  of  polarization  (Fig.  4-1). 


The  absence  of  polarization  dependence  of  the  effective  mass  suggests  that  m* 
is  renormalized  via  an  interaction  with  some  classical  degree  of  freedom  which  is  not 
affected  by  the  quantum  degeneracy  of  the  electron  state.  As  will  be  shown  below,  such 
a  mechanism  is  provided  by  the  renormalization  of  electrons  by  the  virtual  emission  and 
absorption  of  high  energy  plasmons.  We  reanalyze  the  issue  of  spin  independent  effective 
mass  by  considering  the  spin-valley  degeneracy  in  the  2D  system.  A  valley  degeneracy 
provides  an  additional  parameter  (N)  to  the  problem.  This  parameter  is  4  for  the  case 
of  a  2D  Si-MOSFET  (2  spins  (g)  2  valleys)  grown  along  a  (001)  plane.  The  additional 
parameter     =  4  turns  out  to  be  critical  in  explaining  transport  data  in  both  the  diffusive 
and  ballistic  hmits.  We  will  consider  the  effect  of  valley  degeneracy  in  the  context  of  the 
effective  mass  problem  and  will  also  calculate  the  spin  susceptibility.  To  this  end  our 
Hamiltonian  is  defined  as  follows 


where  a,^  =  l....^',  are  the  spin- valley  indices. 
4.2.1    Effective  Mass 

In  the  long-wavelength  limit  the  Matsubara  polarization  operator  at  T  =  0  for  a 
system  with  degeneracy  N  {N  components)  is 


4.2    Multi- Valley  System 


(4.19) 


mN  I  


where  the  Fermi  wave  vector  in  terms  of  the  density  of  electrons  n  is  pp  =  2^yim/N,  i.e., 
the  Fermi  momentum  decreases  by  a  factor  of  \/N.  This  decrease  in  the  Fermi  momentum 
is  due  to  the  fact  that  the  electrons  are  distributed  in  A'^-valleys.  Within  the  random 
phase  approximation  (RPA),  the  typical  momentum  transfer  is  of  the  order  of  k  ~  e^mN. 
The  ratio  a  =  k/pf  =  is  a  parameter  describing  the  crossover  between  the 

regimes  of  weak  (a  <  1)  and  strong  (a  »  1)  screening.  Both  weak-  and  strong-screening 
regimes  are  within  the  reach  of  perturbation  theory  {Ts  <  1),  provided  the  degeneracy 
iV»  1. 

In  the  weak-screening  regime,  the  scattering  processes  leading  to  mass  renormaliza- 
tion  are  essentially  elastic  and  the  typical  momentum  transferred  g  ~  k  is  small.  The 
main  contribution  to  the  mass  renormalization  is  via  the  interaction  of  electrons  with  the 
particle-holes  [shaded  region  in  Fig.  4-5A],  whereas  the  emission/absorption  of  plasmons 
provides  subleading  contributions.  The  effective  mass  is 

^  =  l  +  ^ln(r.iV^/^)  +  0(r.). 

Thus  in  the  weak  screening  regime,  the  ratio  of  the  change  in  effective  mass  in  a  fully 
spin-polarized  case  (A^  to  that  in  the  unpolarized  case  is  rather  large.  This  is 

similar  to  the  earlier  considered  case  of  A''  =  2. 

In  the  strong-screening  regime,  r^A^/^  »  1,  the  t}T>ical  momentum  transfers  are 
much  greater  than  kp.  In  this  regime,  as  the  momentum  transferred  is  large  the  scattering 
is  almost  isotropic  (s-wave).  The  particle-hole  continuum  contribution  to  the  mass 
renormalization  is  greatly  reduced  for  s-wave  scattering.  However,  the  interaction  of  the 
quasi-particle  with  high  energy  plasmons  [shaded  region  in  Fig.  4-5B]  plays  the  dominant 
role  and  leads  to  the  enhancement  of  m*.  As  a  plasmon  is  a  classical  collective  mode,  it  is 
not  affected  by  a  change  in  N.  Consequently,  the  leading  term  in  the  A^^^  expansion  for 
m*  does  not  depend  on  N,  whereas  the  next-to-leading  term  is  numerically  small. 

The  Matsubara  polarization  operator  for  q  »  pp  is 

de  drp  1 


^2  {ie  +  iuj-  ^p+g)iie  -  ^p) 
=   -2n6g/{el  +  uj^),  (4.20) 
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Weak  Screening 


1/21/2 

Figure  4-5:  Particle-hole  spectrum  for  a  multi-component  and  a  2-component  system. 


where     =  p'^/2m  ~  Ep,  eg  =    /2m  and  Ep  is  the  Fermi  energy. 

The  dressed  interaction  lines  consisting  of  the  RPA  diagrams  still  remain  the  most 
dominant,  due  to  the  large  A^-factor  associated  with  each  polarization  bubble.  The 
modified  RPA  potential  with  the  above  polarization  operator 


V{q,  iuj) 


2w£ 

1 


1  + 


27re2  2ne,. 


(4.21) 


has  dominant  contributions  coming  from      ^  rl'^y/Npp,  with  the  screening  radius 
r  ~  (ri^^v^)-!  >  y/n~^  larger  than  the  typical  interelectron  distance  yet  much  smaller 
compared  to  the  screening  radius  r  ~  (r^v/ri)-^  for  the  case  with  only  spin-degeneracy. 
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The  zeroes  of  1  —  {2iTe^/q)U^{q,uj)  determine  the  plasmon  spectrum.  As  shown  in 
Fig.  4-5,  the  long- wavelength  hmit  of  plasmon  spectrum  (proportional  to  y^)  depends 
only  on  the  density  of  electrons  and  is  independent  of  the  spin-valley  degeneracy.  For 
momenta  above  ,  the  plasmon  dispersion  for  the  2-component  system  is  completely 

damped  by  particle-hole  excitations.  In  contrast,  for  multi-component  systems  due  to  the 
considerably  smaller  velocities  of  the  components,  the  plasmons  remain  undamped  until 


large  momenta  of  the  order  s/r^y/Wn.  The  plasmon  dispersion  in  the  large  momentum 
region  is  given  by  ujp  =  ^e'^  +  2'Ke'^nq/m.  The  slow-moving  electrons  in  the  multi- valley 
system  are  better  screened  by  the  high-energy  plasmon  excitations,  with  the  particle-hole 
excitation  playing  only  a  nominal  role.  The  area  of  the  particle-hole  continuum  in  the 
(w,  q)  plane  is  reduced  and  its  contribution  to  the  real  part  of  the  self-energy  (via  ImVfi) 
and  hence  to  the  mass  renormalization  is  correspondingly  small  [16]. 

We  obtain  the  quasi-particle  effective  mass  from  the  self-energy  via  the  standard 
relation 


*  /  1  /'^ 
m  /m  ^  1  -  i- 


+ 


d^p         '  die 

where  AE^(e)  =  E^(e)  -  S^(0)  is  the  difference  of  the  Matsubara  self-energies. 
The  difference  of  the  static  and  dynamic  self-energy  is, 


(4.22) 


AEl,,,(ep)  +  AE^„(^p„z£) 
did  cPq 


I 
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2n  (27r)2 

duj  (fq 
2^  (27r)2 


V{q,0) 
Viq,0) 


G{^p+q,  ie  +  iu)  -  G{^pp+q,  iu) 
{e  +  u;  +  i^p+q){uj  +  i^p^+g) 


m 


2^2^pJ^  d0coseV{2kFsm9/2,0) 


(4.23) 


ASt(ep,^•e)   =  lt-0f2iyi^M-V{q,O)] 


G(^p+q,  ie  +  iu)  -  G{^pj,+q,  iu) 


I 


du  (fq 
2^(2^ 


{V{q,  iu)  -  V{q,  Q)){ip+g  -  ip-,+q  -  ie) 
{e  +  u  +  i^p+g){u  +  i^p^+g) 


(4.24) 
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Figure  4-6:  Diagrammatic  corrections.  A)  The  RPA  potential.  B)  Corrections  to  the 
polarization  bubble  due  to  the  interaction  terms.  C)  Vertex  corrections  to  the  self-energy. 


The  pole  at  a;  =  ireg,  where  r  =  y^l  +  ine'^n/qeg,  provides  the  dominant  contribution  to 
the  frequency  integration  of  the  dynamic  term  in  Eq.  (4.24).  Taking  the  derivative  of  the 
dynamic  term  we  obtain 


dynamic 

die 


dynamic 


—  —  2n7re^ 


/ 


qdq 


r{r  +  lyel 


,/3r(i/3)r(7/6) 


(4.25) 


The  result  of  Eq.  (4.25)  is  larger  by  a  factor  of  {rsN^''^fl^  compared  to  that  from  the 
static  part  and  will  contribute  to  the  leading-order  correction  to  the  effective  mass.  Thus 
the  leading  term  in  the  1/iV-expansion  is  (this  form  was  obtained  earlier  by  lordanskii  and 
Kashuba  [90]) 


m 
m 


o/3r(i/3)r(7/6) 


(4.26) 


The  form  of  the  polarization  operator  as  given  by  Eq.  (4.20)  is  insensitive  to  the  degree 
of  polarization,  thus  the  mass  renormalization  in  the  leading  order  remains  unchanged 
even  for  full  polarization.  Corrections  to  this  leading  order  term  is  obtained  by  including, 
A)  corrections  to  the  polarization  bubble  due  to  the  contributions  coming  from  small 
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momentum  regions,  B)  vertex  corrections  to  the  bubble  [Fig.  4-6B],  and  C)  vertex 
corrections  to  the  self-energy  [Fig.  4-6C]. 

From  a  scaling  analysis  we  can  infer  that  the  vertex  corrections  from  B)  and  C) 
contribute  to  a  rt^^  term  to  the  leading  result.  The  correction  from  B)  is  Sm^/m  oc 
J du(PqG'^V^6U.,  where  V  ~  £qo/n,  q      qo,  lo  ~      and  SU  ~  rV^n/eqa,  this 

yields  a  ri^^  correction  term.  The  corrections  from  C)  is  6ml/m  oc  /  dui(fqidu)2(fq2G'^V'^, 
where  Ui  ~  W2  ~        ~  ^90;  9  ~  9o  and  V  ~  Sgo/n,  this  also  results  in  a  rt^^  term. 
We  have  checked  exphcitly  that  the  scaling  analysis  holds.  These  corrections  remain 
independent  of  the  degree  of  polarization.  We  will  next  consider  contributions  from  A) 
which  gives  polarization-dependent  corrections  and  form  the  next-to-leading  term  in  the 
1/A^-expaiision. 

Contributions  from  the  small-momentum  region  to  the  polarization  bubble  are  given 
by  Eq.  4.27, 


6U^{q,zu)   =   27V I  ^-^-^"^ 


(27r)2 


(3a;%-g3)(ffl^)2 

((^2+^2)3 


(4.27) 


mN   (a;2  +  £2)3  ' 

The  correction  in  the  polarization  bubble  modifies  the  position  of  the  plasmon-pole  from 

uj  =  irsg  to  LU  =  ir*eg.  Where  r*'^£l  =  r^s^  +  A^,  with  rsg  >  A  and 

qmN  (u;2  +  e2)2~  (^^2-62)2-  (4-28) 

Thus  the  additional  corrections  to  the  mass  term  due  to  the  new  position  of  the  pole  will 
be 


3r  -  2r2 


die  ^  d^p 


_2_  rdq  I +  31 
N  J  "7  r3(l 


+  r)3  '  (4.29) 


where  the  partial  derivatives  of  6E  are  taken  at  (e,^p)      0.  The  UV-divergence  in  the 
above  expression  is  cut  off  at  typical  momentum  transfers  q  ~  \/rJP/^pF,  as  at  these 
momentum  the  plasmon  spectrum  is  damped  due  to  the  particle-hole  excitations.  The 
leading-order  corrections  to  the  effective  mass  from  the  corrections  to  the  polarization 
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Figure  4-7:  Percentage  change  in  effective  mass  between  a  fully  polarized  and  an  unpolar- 
ized  system  as  a  function  of  r^. 


bubble  will  thus  be 


d 


d 


-die  ^  d^p 


=  — L  log  iVr2/3. 


This  contribution  is  larger  by  a  factor  of  log  rl^^N  compared  to  that  from  the  static  part 
[Eq.  (4.23)]  and  is  the  dominant  subleading  correction  to  the  effective  mass.  Thus  we 
obtain  the  mass  to  next-to-leading  order  in  l/N  as 


^  =  1  +  0.14r,^/3  +  ^  logr.Ar3/^  +  o(_i_). 


m 


(4.30) 


Because  of  the  small  coefficient  in  front  of  the  leading  term  in  Eq.  (4.30),  the  actual 
condition  for  small     is  rather  weak:  a  two-fold  enhancement  of  the  mass  occurs  only  for 
rs  ^  20,  thus  smaller  values  of     allows  a  reasonable  description  within  the  mean-field 
theory.  Equation  (4.30)  practically  coincides  with  the  observed  m*{rs)  dependence  in  the 
interval     =  2  -  6.  Beyond     =  6  the  experimental  value  (Pudalov  et  al.  [66])  exceeds 
the  theoretical  value  (Fig.  4-8).  For  full  polarization  N  -»  N/2,  unlike  the  earlier  model 
for  a  2-component  system  (spin-degeneracy  only),  the  percentage  change  in  the  ratio  of 
effective  mass  for  a  fully  polarized  to  the  unpolarized  case  is  small.  This  ratio  is  given  by 
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Figure  4-8:  Plot  of  the  effective  mass  formula  obtained  from  the  weak-screening  case  and 
from  the  strong  screening  case. 


Eq.  (4.31), 

Am        m*{^  =  1)  —  m* 


■m%.  (4.31) 


rriavg     {m*{^  =  l)  +  m*)/2 
For  an  Si-MOSFET  sample  N  changes  from  4      2  at  full  polarization,  the  calculated 
effective  mass  remains  almost  the  same  and  for  the  range  of     from  2-5,  Am  is  not  more 
than  1-3%  (Fig.  4-7),  a  similar  result  is  obtained  for  iV  =  2  ^  TV  =  1.  The  difference  in 
effective  mass  is  particularly  small  for  the     =  4  ->     =  2  degenerate  system  due  to  the 
peculiar  cancellation  of  the  log  N/N  term. 

The  exphcit  polarization  dependence  of  the  effective  mass,  given  by 
m* 

 =  i  +  U.14r:'"  H  —  log!  - 

m  '         8N  ^Vl+^^ 

is  plotted  in  Fig.  4-9.  The  effective  mass  remains  essentially  constant  throughout  and 


1  _L  ^2         ,  „2/3  ^ 

-.l  +  0.14.V3  +  l^,„g(ri^),  (4.32) 
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Figure  4-9:  EfTective  mass  as  a  function  of  degree  of  polarization  for  different  values  of  r^, 
starting  from  the  lower  one     =  2  to     =  3, 4  and  5. 

there  is  no  spin-sphtting  of  the  masses  at  this  order  in  1/A''  Eq.  (4.32)  is  valid  as  long  as 
the  separation  between  two  spin-down  electrons  is  less  than  the  screening  length.  The 
condition  for  the  applicability  of  mean-field  theory  is  1  —  ^  »  r^^^  ] .  This  result  is  in 
qualitative  agreement  with  the  experimental  finding  on  Si-MOSFETs  [68],  where  the 
effective  mass  showed  no  dependence  on  the  degree  of  polarization. 
4.2.2    spin  susceptibility 

The  second  Fermi-liquid  parameter  which  is  of  great  interest  in  the  context  of  a 
possible  ferromagnetic  transition  and  is  experimentally  determined,  is  the  Lande-^*  factor. 
The  g*-factor  is  extracted  from  a  knowledge  of  the  static  spin  susceptibility  x*  via  the 
relation  x*/x  =  'ni*9* /mg.  In  the  previous  Section,  we  obtained  the  effective  mass  as  a 
function  of  interax;tion  strength     and  the  valley  degeneracy  parameter  N.  As  before,  we 
will  use  perturbative  techniques  to  obtain  the  spin  susceptibihty  as  a  function  of  the  above 
two  parameters. 

The  spin  susceptibility  is  obtained  by  studying  the  response  of  the  system  to  the 
external  magnetic  field  and  in  the  Matsubara  formalism  is  given  by: 


(4.33) 


^..x..  98 
where  m  =  9Ms(V'qC^q^^^/3)/2  is  the  magnetization,  a  is  the  Pauh  matrix,  \ib  is  the  Bohr 
magneton,  a,  h  are  the  vector  indices  and  g    2  is  the  Lande-^  factor.  Plugging  the  form  of 
the  magnetization  in  Eq.  (4.33),  one  obtains 

X^^'^H^^n)  =  f^^  e-"^(r,^,(r)ai">^(r)^(0),a5V*(0)5(l/r,  0))dT,  (4.34) 

In  the  zero*'*  order  of  interaction  S{l/T,  0)  =  1,  the  electron  spin  susceptibility  is  given 
(using  "Wick's"  theorem)  by 

Xo^°'\i-^n)  =  {^^f  f^'^  e-"^G(-r)G(T)TVace  {a^'^^a^'^)dT.  (4.35) 

Since  Trace  (cr^^'a^^^)  =  25ab,  the  relation  is  simply  expressed  in  terms  of  the  single  particle 
Green's  function  G. 

Xo^''\i^n)  =  ^-^Sa,  f^'^  e'--^G{-T)G{T)dr,  (4.36) 

Thus  Xo^''''\ii^n)  =  ^^^Sabn{iuJn),  where  n(za;„)  is  the  usual  Lindhard  function. 

Due  to  electron-electron  interaction,  the  susceptibility  has  the  Lindhard  form  only 
in  the  lowest  order  of  interaction.  In  what  follows  we  obtain  corrections  to  the  free 
susceptibility  by  including  interaction  effects. 

Higher-order  corrections.  We  incorporate  higher  order  terms  by  including  the  correc- 
tions due  to  the  vertex  function  as  in  Fig.  4-6  (Bl)  and  due  to  the  self-energy  insertion 
to  one  of  the  Green's  function  as  in  Fig.  4-6(B2).  The  effective  interaction  in  both  cases 
considered  is  the  RPA  potential.  Before  evaluating  the  additional  contributions  to  the 
susceptibility,  let  us  briefly  discuss  the  relevant  range  for  the  momentum  and  frequency 
transfer.  The  typical  momentum  transfer  k  in  the  dense  single  component  limit  is  of 
the  order  of  Vskp  which  is  much  less  than  kp.  The  case  of  an  A^- valley  system,  however, 
greatly  differs  from  the  single  component  limit.  Due  to  the  presence  of  N  valleys,  the 
electrons  get  distributed  over  the  valleys  and  hence  the  Fermi  energy  of  each  valley  gets 
suppressed  by  a  factor  of  l/y/N.  Typical  momentum  transfers  becomes  much  greater 
than  the  Fermi  momentum  for  the  interaction  regimes  such  that  iV"^/-  <     <  1.  Con- 
tributions from  the  regions  with  small  momentum  transfers  k  <^  kp,  me  sub-leading  to 
those  from  the  large-momentum-transfer  regions.  Details  of  the  calculation  are  provided  in 


99 


Appendix  C.  The  correction  to  the  Lindhaxd  form  of  the  susceptibiUty  is 


n*   =  Uo  +  5U 


f  qidqi  duji 


qf/m 


m* 


(4.37) 


Thus  from  this  form  of  the  susceptibiUty,  we  can  say  that  the  corrections  to  the  polariza- 
tion due  to  small  momentum  regimes  will  give  the  same  -1/NlnN  correction  as  the  mass 
term.  The  g*  factor  remains  unaffected.  This  result  is  in  qualitative  agreement  with  the 
experiment  on  Si-MOSFETs.  A  similar  behavior  is  exhibited  by  ^He  where  the  renormal- 
ization  of  the  spin  susceptibility  is  essentially  due  to  the  effective  mass.  Near  the  critical 
point  the  spin  susceptibility  of  ^He  diverges,  due  to  ^He  being  "almost  locaUzed"  (m* 
diverges)  rather  than  being  "almost  ferromagnetic"  (where  g*  diverges).  It  can  similarly 
be  argued  that  the  divergence  of  the  spin  susceptibility  of  Si-MOSFETs  near  the  critical 
point  is  due  to  the  divergence  of  effective  mass  rather  than  the  g*  factor. 

However,  a  recent  experiment  on  g*  (AlAs  sample)  has  shown  that  the  g*  factor  is 
affected  by  lifting  the  valley  degeneracy  [91].  We  beheve  more  work  is  required  to  qualify 
this  behavior  to  a  many-body  effect. 


Theoretical  arguments  based  on  general  grounds  predict  that  in  a  generic  spin- 
polarized  Fermi  liquid,  the  effective  mass  of  spin-up  and  spin-down  electrons  will  be  dif- 
ferent. Contrary  to  these  expectations,  the  effective  mass  of  electrons  in  2D  Si-MOSFETs 
remains  polarization-independent  even  under  full  spin-polarization.  This  indicates  that 
Fermi  statistics  do  not  play  a  role  in  the  mass  renormahzation  and  that  the  mass  renor- 
mahzation  is  due  to  classical  degrees  of  freedom.  For  a  system  with  spin-valley  degeneracy 
N  and  in  the  strong  screening  regime  such  a,  polarization-independent  effective  mass 
is  shown  to  result  from  the  interaction  of  electrons  with  the  high  energy  plasmons.  In 
the  strong  screening  regime  plasmons  are  the  main  contributors  to  the  mass  renormal- 
ization  and,  since  they  are  not  affected  by  the  change  in  A^,  the  effective  mass  remains 
polarization-independent  to  the  leadmg  order  in  1/A^.  We  used  a  large  A-expansion  and 
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showed  that  the  leading  correction  is  independent  of  N,  the  information  about  the  degen- 
eracy A''  comes  only  in  the  next  to  leading  term.  The  expansion  of  the  effective  mass  in 
1/A''  turns  out  to  be  rapidly  convergent  even  for  non- valley  degenerate  systems  (A^  =  2). 
In  particular,  we  find  that  within  this  model,  the  percentage  change  in  the  effective  mass 
of  a  2- valley  system  under  full  spin-polarization  is  not  more  than  1-3%,  a  result  consistent 
with  experiment. 

The  spin  susceptibility  within  this  model  is  almost  entirely  renormaUzed  by  the 
mass.  This  finding  is  also  in  qualitative  agreement  with  the  experimental  results  on 
Si-MOSFET's  . 


CHAPTER  5 
CONCLUSIONS 

In  this  work,  we  have  investigated  the  origin  of  processes  which  lead  to  the  deviation 
from  Fermi-liquid  universality  in  two  dimensional  interacting  systems. 

In  the  first  part  of  this  work  we  have  analyzed  the  low-energy  properties  of  Fermi 
liquids  in  two  dimensions.  We  have  found  that  even  at  second  order  in  perturbation 
theory  (for  a  short-ranged  potential),  two  interesting  results  are  obtained:  a  nonanalytic 
term  in  the  imaginary  part  of  the  self-energy  and  a  logarithmic  singularity  on  the  mass 
shell  for  fermions  with  a  linearized  dispersion.  The  nonanalytic  term  and  the  log-singular 
term  originate  from  special  scattering  events.  Scattering  events  that  involve  exchange  of 
small  momentum  between  two  parallel-moving  quasi-particles,  also  known  as  forward- 
scattering  events,  give  rise  to  the  log-singular  term  on  the  mass  shell.  We  find  that  at 
higher  orders  in  the  interaction  the  singularity  due  to  forward  scattering  processes  is 
enhanced  and  becomes  a  power-law  causing  the  perturbation  theory  to  formally  diverge. 
We  have  dealt  with  the  divergences  by  a  re-summation  process  which  involves  adding 
up  all  the  maximally  divergent  diagrams.  An  outcome  of  this  resummation  process  is  a 
zero-sound  mode.  We  find  that  the  interaction  of  the  fermions  with  the  zero-sound  mode 
results  in  a  non-Lorentzian  shape  of  the  spectral  function. 

We  have  investigated  the  contribution  of  the  zero-sound  mode  to  the  nonanalytic 
r2  term  in  the  specific  heat.  We  have  found  that  for  a  short-ranged  potential,  both  the 
imaginary  and  the  real  part  of  the  self-energy  emanating  from  the  zero-sound  mode  re- 
gion contribute  to  the  nonanalytic  7^  term  in  the  specific  heat.  In  fact  the  contributions 
from  the  real  and  imaginary  parts  exactly  cancel  each  other,  thus  we  conclude  that  the 
nonperturbative  zero-sound  mode  does  not  contribute  to  the      term  in  the  specific  heat. 
We  have  verified  the  result  by  calculating  the  specific  heat  directly  from  the  thermody- 
namic potential  using  both  the  real-  and  imaginary-frequency  formalisms.  This  work  was 

extended  to  the  case  of  finite-ranged  potential.  The  main  result  here  is  that  the  term 
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can  be  represented  in  terms  of  the  Fermi-liquid  parameters.  We  have  also  considered  the 
case  of  the  long-ranged  Coulomb  potential  and  show  that  nonperturbative  processes  do 
contribute  a      term  to  the  specific  heat  and  for,  sufficiently  low  T,  this  contribution  is 
interaction  independent. 

Direct  evidence  for  the  Fermi-liquid  behavior  (or  its  potential  breakdown)  of  a  two- 
dimensional  electron  system  comes  from  the  analysis  of  Shubhnikov-de  Haas  oscillations  of 
the  conductivity  in  the  presence  of  both  perpendicular  and  tilted  magnetic  fields.  A  recent 
experiment  [68]  on  Si-MOSFET  has  revealed  an  unusual  feature  of  the  electron  effective 
mass  (m*):  while  m*  exhibits  a  strong  dependence  on  the  electron  density  (r^),  it  does 
not  depend  on  the  degree  of  spin  polarization.  Also,  the  masses  of  electrons  with  up  and 
down  spins  are  the  same.  We  have  made  use  of  two  independent  approaches,  perturbation 
theory  in  the  weak-couphng  limit  and  Landau's  phenomenological  approach,  to  show  that 
the  experimental  results  are  in  apparent  contradiction  with  Fermi-hquid  theory,  which 
predicts  two  different  and  polarization-dependent  masses  in  a  partially  spin-polarized 
regime.  To  explain  the  experimental  observation,  we  have  considered  the  effective  mass 
of  the  Coulomb  gas  in  the  large-N  limit  (for  Si  MOSFET,     =  4).  We  show  that  the 
effective  mass  is  renormaUzed  primarily  due  to  a  polaronic  effect:  emission  and  absorption 
of  high  energy  plasmons.  As  plasmons  are  classical  objects,  the  quantum  degeneracy, 
and  hence  polarization,  does  not  affect  the  effective  mass  to  the  leading  order  in  1  /N . 
Polarization  dependence  shows  up  only  at  the  next-to-leading  order.  We  find  that  for 

=  2  —  6  the  change  in  effective  mass  between  the  unpolarized  and  fully  polarized  states 
is  within  1  -  3%,  a  result  which  is  consistent  with  the  experiment. 


APPENDIX  A 

SPECIFIC  HEAT  CONTRIBUTION:  COULOMB  POTENTIAL 
In  the  following  analysis  we  will  reconsider  the  plasmon-pole  and  particle-hole  con- 
tribution to  the  specific  heat  from  an  alternate  approach.  The  thermodynamic  potential 
given  in  Eq.  (3.48)  can  be  rewritten  in  the  following  form 

A  2       7  (27r)2  I         kl-  V{k)U{k,  iWn) '  ^  ^ 

The  frequency  summation  can  be  converted  to  frequency  integration  by  using  the  Ehash- 
berg  approach.  The  contour  encloses  all  the  discrete  frequencies  except  iWn  —  0,  the 
branch  cut  and  the  two  poles.  The  branch  cut  corresponds  to  the  particle-hole  region,  for 
a  given  value  of  k  it  extends  from  —w{k)  to  +w{k)  (Fig.  3-5).  The  origin  of  the  branch 
cut  is  due  to  the  11  term  which  satisfies  n(A;,a;  -I-  iS)  7^  ]l{k,uj  -  id)  for  a;'s  such  that 
u  <  u!{k),  where  u}{k)  =  kvp  +  k'^/2m  is  the  boundary  of  the  particle-hole  continuum.  The 
poles  correspond  to  the  solution  of  Soik,  w)  =  0,  hence  the  solutions  are  the  plasmons  for  a 
system  of  charge  €„,  which  is  beyond  the  particle-hole  region.  For  an  arbitrary  value  of  el 
(<  e^)  the  pole  hes  between  the  particle- hole  and  the  plasmon  curve.  We  can  express  the 
firequency  summation  as  follows 

/^trl-^n(fc,m„)         [   J_^,,,  27ri''''^'>i-2^U(k.e  + 


w(k)  27ri       'i_^n(/j,e  +  z^) 

r+w(k)  J 

-Pj  ^ns{e) 


de  U{k,£-i5) 

7: — nBis)  

^(k)  27ri  'i-?^u{k,£-iS) 

1  Il{k, —Wpi{k)) 

^e-0w„i(k)  _  i'       27re2  dn{k,E)  U^-Wpi{k) 
k  de 

I  (  1  Mk,Wj,i{k))  1 

k  de 

The  first  expression  in  the  summation  corresponds  to  the  particle-hole  region,  where  as 
the  second  region  would  correspond  to  the  plasmon-pole  region  (where  Upi  =  VFVme'^k). 
The  above  expression  can  be  written  in  a  further  simplified  form  by  making  use  of  the 
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Mowing  properties  1).  n«(fc,£)  =  [U^{k,e)]*,  2).  U{k,±Wpi{k))  =  1/V{k)  and 

dU{k,e). 


3). 


2  2 

m  9 


is  odd  in  Upi.  The  simplified  form  is 


U{k,iWn) 

-VU{k,iWn) 


=  P  — nB(£:)Im  

V-u;(fc)     TT  1  - 


U^{k,e) 


27re2 


^n«(/c,e) 


+ 


■       (3wpi{k)  1 


(A.3) 


thus, 


n(fc,?u;„) 

-yn(A;,ra„) 


d_ 

dT' 


"  /        — nB(e)Im— 

J-w{k)     ^  1 


n''(A;,£) 


+ 


27re2 


1 


n«(A;,e) 


^^COtn^      2      ^^2^e2^2gnfi(fc,e)  l^-""pKfc) 


/+«)(*;) 


1 

4T2 


-Im- 


n«(A;,e) 


+ 


TT  sinh^(£/2T)  i-2^n«(A:,£) 


As  before,  we  will  first  calculate  the  contribution  from  the  Plasmon-pole, 
dQpi/A 


dei 


dT  IGTrmvlT^Jo       Jo      "  sinh^Upi /2T)ei' 

where  the  fc-integration  is  limited  till  ko  =  me^.  At  ko  the  plasmon-pole  enters  the 


(A.4) 


particle-hole  spectrum  and  is  no  longer  free.  This  integral  is  easily  evaluated  by  making 
use  of  the  result 

y'dy  3C(3) 


^0 


we  obtain 


/o  sinh2(y)  2 
dQpi/A   _  3C(3)r2 


(A.5) 


dT 

Cv-pi 


2-Kvl 


d_dnpijA 

dT  dT 
3C(3)r^ 


(A.6) 


The  remaining  contribution  is  from  the  particle-hole  region, 
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-Y  — 


IWn 


VU{k,iWn) 


PH 


i 


xlm 


TT  2T^smh^{e/2T) 
n«(A;,e) 


l_2^n«(A;,£)' 


(A.7) 


where  the  Im  term  is  odd  in  s.  Plugging  this  in  Eq.  (A.l),  we  obtain 


dT 


PH 


de 


-Im- 


n  2T2  sinh2(e/2T)     1  -  ^m{k,  e) ' 
next  performing  the     integration  we  obtain  an  identical  expression  as  in  Eq.  (3.21), 
with  the  frequency  and  the  momentum  integration  limited  to  be  on  the  PH  region.  The 
evaluation  of  this  term  proceeds  exactly  as  before  and  we  will  not  pursue  it  any  further. 


APPENDIX  B 
SPIN-DEPENDENT  EFFECTIVE  MASS 

In  this  Section  we  will  calculate  the  effective  mass  of  spin-polarized  electrons.  The 

first  term  of  Eq.  (4.14)  is  frequency  independent,  thus  it  does  not  contribute  to  the 

dT,^^/diu)  term  in  the  effective  mass.  We  need  to  consider  only  the  2"''  term  of  Eq.  (4.14) 

in  the  self-energy  expansion 


^    '^''^^  ^  'L  WT 


kdk    V^{k)n{k,  iu)  ''^^ 
Vik)U{k,iu) 


f 

Jo 


dOG]L{p  +  k,iij  +  iuj) 


'(0) 


f+°°dw  r  kdk   V\k)l{{k)    f'\A^n^     ,  ■     •  N 


+  G\l^{p  +  k,  -iu)  +  iu) 


(B.l) 


The  angular  integration  yields 


(p  +  k,  -iu  +  iu)   =  - 


d^G(Q)(p  +  k,  -iu  +  iu) 


'2m  r0(-a;  -\- u)  -  Q{u  -  u) 


4 


1  + 


(B.2) 


We  will  differentiate  E^^  by  taking  the  limit  a)  ^  0+.  The  result  remains  the  same  if 
instead  we  had  taken  the  limit  a;  ->  0~.  The  6  functions  are  defined  as  follows 


G(x)  =  1;  for,  x  >  0; 
e(a;)   =  0;    for,  x  <  0. 


Hence  S^^^  can  be  written  as 


2m  r 


lU] 


+  2 


Jo   2tt  Jo  (27r 


^V{k,iu) 


2m 


^(1)  +  ^(2)  ' 


J 


(B.3) 
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where  V{k,  iu)  has  been  used  to  represent 


Vik,iuj) 


V^{k)U{k,iLO) 


l-V^{k)U{k,iuj) 


The  contribution  from  the  first  term  is: 


(1) 


diu 


<i=0 


2nm 


We  will  make  the  following  substitution  in  the  above  integral  uj  —>■  uj\jAuJkE^  ,  k  ->  kk]} 
and  =  (1  -  ^)/(l  +      where  ^  is  the  polarization  (we  will  assume  ^  «:  1).  This 

results  in  the  following  equation, 


"(1) 


diu 


u/=0 


-  r 

Trplr  Jo  Jo 


dk 


UJ 


k  (l+u;2)5 


+ 


+ 


2) 


(B.4) 


Expanding  Eq.  (B.4)  about  5,  where  1  +  ^Ti  ^     ^        ±  0  =^      ~  ^2^,  we  obtain 


^^(1) 


diuj 


U!=0 


=  /     duj  / 

7r\/2  io  io 


a; 


k  (l  +  a;2)| 


- 1 


rr  +  l 


dk- 


(l+a;2)3 


(: 


+  1  - 


(B.5) 


To  prevent  ultra-violet  singularities,  an  upper  cut  off  on  the  A;-integral  at  A;  =  A;^  is 
imposed.  Such  a  cut-off  is  needed  since  the  form  of  E  was  obtained  by  considering 
contributions  due  to  small  /c-values.  Thus  we  have 


diw 


Q=0 


=  /      duj  / 

jTi     /-oo  /.I 

/    duj  I  di 
Stt  70  Jo 


dk 


UJ 


k  {1  + 


■7^ 


W  2 


- 1 


'  %/2i 


(l+u;2)3 
+  1 


(B.6) 


The  second  term  of  S^^^  represented  by  E^^''^  gives  the  following  contribution 

d    r  duj  ['^  kdk    V'^{k)T[{k,uj)  2m 


diljj 


ii->0 


2—  r  —  /■°°_^_ 


V{k)Uik,u)p^  n 

m  Y 

dk 


Adding  up  S^j'^  >       result  assumes  the  following  form 


^(2) 


diu} 


diuj 

UD-.0 


as?!! 


'  dk 


+ 


i±-  r  /  —  — 

+—  /  dw— 
Sn  Jo  (1 


iu!— >0 

1  -  -7^ 


\/2: 


^  +  1- 


+ 


+  a;2)3 


1  - 


\/2i 


simplifying  it  further  and  keeping  terms  of  the  order  of  r]^,  Eq.  (B.8)  reduces  to, 


diu 
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(B.7) 


(B.8) 


(B.9) 


Next  we  will  calculate  dT,/d^p.  Unlike  the  earlier  case,  where  we  needed  to  perform 
frequency  differentiation  on  only  the  H^^lj,  both  terms  E^^^  and  E^^^,  need  to  be  differenti- 
ated with  respect  to  ^p.  The  contribution  due  to  the  first  term  is: 

/d^k  /"°° 
^F(A;)y    duGll^ip  +  k,iu  +  iQ),  (B.IO) 

completing  the  frequency  integral  we  have 

=  -/(02^W©(f^'-b  +  ^l) 

2Lj]}r'^^  f°°  dx 
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notice  we  have  expressed  the  above  integration  variables  in  dimensionless  units.  The  0 
function  imposes  a  restriction  on  the  range  of  values  of  x  (=  cos  6)  and  k.  The  restriction 
are  as  follows 

1  -        +  p2  _  2kpx  >  0 
=>x>      +  l)/2fcp, 
and  the  restriction  on  a;,  -1  <  x  <  1,  imphes  {k'^  +     —  2kpx)/2kp  <  1 

^p-l<k<p+l 

Finally  when  we  take  the  derivative,  dH^^^/d^^^,  we  will  take  the  limit  p  — »  1+  ;  thus  the 
range  of  values  k  can  take  is 

p-l<A:<p  +  l, 
and  the  range  of  values  x  can  take  is 

k'+p'-l 

Since  (A;  +  p)^  >  1  the  lower  limit  of  x  in  Eq.  (B.ll)  can  be  replaced  from  -1  to  (k"^  +  p^  - 
l)/(2pfc). 

Evaluating  the  x  integration,  we  obtain 

2ijl^rl^  ..  dx 


Jp-l      V2pfc  -  (p2  +     -  1) ^2pk  +  (p2  +  k^-l) 
substituting  k  =  t+pmEq.  (B.12),  we  obtain 

E"'  =  /'  d<-,=£±i£±i=.  (B.13) 

Next  we  will  differentiate  the  above  term  with  respect  to  the  momentum,  thus  we  obtain 


dp 


84M^  f\^_  P  +  t 


nV2  7-1 


(B.14) 
p=i 


no 


substituting  ^  =  2it  -  1  in  Eq.  (B.14),  Eq.  (B.14)  modifies  to 


dp 


7rV2   io  u(l +u)(l 


du 


p=i 


Ti 


24V]i  r 

7rv/2 


log2-  /   ^+  /  —  .  B.15 

Jo  u{l- 11^)2     Jo    u\  ^  ' 


4i  r 


log2-l  + 


(B.16) 


The  infra-red  divergence  of  Eq.  (B.16)  is  cancelled  by  a  similar  divergent  term  originating 
from  the  SE^^V^Cp^  term.  The  contribution  from  the  E^^^  term  is  calculated  below.  The 
contribution  to  dE^^/d^pn  ,  due  to  the  second  term  E^^^  is  given  by 

j^2n  ^  _         r  kd^_v\k)n{k)  '•^^ 

Jo      2itI  {27r)H-V{} 


V{k)U{k) 


f 

Jo 


d9 


+  G]L{p  +  k,  -iu)  +  iu) 


(B.17) 


We  will  a  priori  set  the  external  frequency  w  =  0,  the  angular  integration  gives 


/■27r 


de 

2mm 
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pk 
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2mm 


(B.18) 


we  differentiate  the  above  expression  with  respect  to       to  obtain 
dll^  _      __.r    Aujk~2{uJk  +  (l})  +  2iu 


n 


'P  - 


-m 


+ 


4iOk  -  2{u;k  +  ^l^)  -  2iu 


(B.19) 


By  linearizing  the  spectrum  (discarding  Uk  terms  compared  to  ~  kvF  terms),  Eq.  (B.19) 
acquires  a  simpler  expression 


i^pUUk+UJ^' 


(B.20) 


Ill 


Using  the  above  relation,  we  obtain  the  following  form  for  5E^^V^^p^ 

5S2Ti         rdw  r  kdk    V^{k){I{^^{k,iu)  +  Y[^{k,iuj)) 


Ti 


f'^  dw  kdk 


V{k){nn{k,iLu)  +  m{k,iuj)) 


(B.21) 


Similar  change  of  variables,  as  in  Eq.  (B.12)  is  made,  w  — >  iwkw]}  ,  k  — >  kk]}  and 
(wp/wp)'^  =  (1  —  0/(1  +  0  =  1  +  expanding  in  the  small  parameter  6,  and  keeping 
terms  of  the  order  of  r,  we  obtain 
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7r>/2L 


dk  TT 

Jo 


+  1  -  log(v^r]^) 


f 

Jo 
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TT  Jo 


UJ 


(l+u;2)3 


UJ 


x/TT 
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duj 


(l+u;2)3 


4i 


+  ... 


(B.22) 


The  infra-red  divergence  is  eliminated  once  we  add  contributions  from  dT}^^ /di^} 
[Eq.  (B.16)]  and  dT?^^/d(l}  [Eq.  (B.22)], 


5^ 


4i 


TT       ,         1  1 


\  r 

TT  Jo 


duJ 


(l+a;2)3 


vot) 


r 

1  ^ 
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V2  V  v^TT^ 


+ 


.(B.23) 


Finally,  results  from  Eq.  (B.9)  and  Eq.  (B.23)  are  plugged  back  into  Eq.  (4.13)  to  obtain 
the  effective  mass  in  terms  of  the  interaction  strength 


m 


m 


(B.24) 


APPENDIX  C 
CORRECTIONS  TO  THE  SPIN  SUSCEPTIBILITY 

In  this  Section,  we  will  evaluate  the  corrections  to  the  spin  susceptibility  due  to  the 

contributions  from  vertex  and  self-energy  insertions. 

Vertex  function.  We  will  go  beyond  the  Lindhard  form  of  the  susceptibility  and 

calculate  the  corrections  to  the  spin  susceptibility  from  the  vertex  function  (Fig.  4-6(Bl)), 

2       r  r  roo 

6Ui{q,iuj)  =  J  d^q\  J    dui  J  (fp  J    deGM{p,i£)GM{p  -  q^ie  -  iui) 

Gm{v-  q\  -q,i£  -  ii^i  -  iuj)GM{p-  q,i£  -  iuj)D{qi,ui),  (C.l) 

where  q,  uj  are  the  external  momentum  and  frequency  transfer  respectively.  The  effective 
potential  is  y  ■  * 


  l-n{q,uj)V{qy 

where  V{q)  =  l-K^jq  is  the  2D  Coulomb  potential  and  ^ 

2n(gV2m) 

is  the  polarization  operator  in  the  limit  of  large  momentum  and  frequency  transfer.  The 
external  frequency  a;  will  be  set  to  a;  =  0,  as  we  are  interested  in  the  static  and  long 
wave-length  form  of  the  spin  susceptibility,  though  we  refrain  from  setting  g  =  0,  and  take 
the  long  wave-length  hmit  at  the  very  end  of  calculation.  Thus  we  have, 

2  /•  /"OO  n  poo 

6Ui{q)   =  J  d^qi  J    duJi  J  d^p  J    deGMip,i£)GM{p-  qi,ie  -  iui) 

Gm{p-  qi  -q,ie-  ioJi)GM{p-  q,ie)D{qi,uJi).  (C.2) 

Integration  on  the  frequency  variable  e  is  performed  by  closing  the  contour  of  integration 
from  above. 
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We  obtain 

e(-ep) 


+ 


+ 


+ 


 e(-ep--,-j  

 Q(~^p-qi-q)  

(<^1  -  'i^p-qi-q  +  i^p){'iip-qi  "  i^p-q-i-q){^l  "  ^^p-q-qi  +  i^p-q) 

Qi-ip-q) 


(C.3) 


By  appropriate  choice  of  variables,  the  9  function  is  made  to  have  the  same  variable  in  all 
the  four  expressions  inside  the  big  bracket, 
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i^p  +  i^p^q) 

0(-e->) 

{ui  +  i^. 

-P+Oi  ~ 

-p){i^-p-q  -  i^-p)i^i  +  i^- 

-p+qi-q  ~  i^-p) 

e(-e-p) 

{ux  +  ii- 

-p+qi+q  ~ 

'^-p){i^-P+q  -  i^-p){^l  - 

^-p  +  i^-p^-ql ) 

e(-e-p) 

(C.4) 


{i^-p  -  i^-p+di^^i  +  i^-p  -  i^-p+q-q-i){'i'i-p-q-i  -     -  i^-p) 
The  above  expression  can  be  represented  in  a  compact  form,  for  this  purpose  we  utilize 
the  following  property  of  the  integral.  The  integral  remains  invariant  under  the  following 
change  of  variables  ui  — <•  —  Wi,^  — >  — >  —p  [where  we  have  utilized  the  property  of 
Matsubara  potential  D{qi,uJi)  =  D{\qi\,  \oJi\)].  Thus  we  obtain 

e(-ep)Wi,a;i) 


(UJI  +  i^p  -  i^p-q-,){u)i  +  iip  -  i^p-qi-g) 


(C.5) 


{-iip+iip-q) 

This  contribution  along  with  the  contribution  from  the  self-energy  insertion  will  be  added 
together  in  the  following  analysis,  to  obtain  the  final  result. 
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Self-energy  insertion.  In  this  Section,  we  will  calculate  the  corrections  to  the  spin 
susceptibility  from  the  self-energy  insertion  to  the  green's  function.  As  shown  in  Fig.  4- 
6B,  there  are  two  such  diagrams  which  contribute  equally  to  the  corrections.  The  effective 
potential  as  before  is  D{q,iu)  =  V{q)/ (1  -  Ii{q,uj)V{q)),  where  V{q)  =  2-Ke'^/q  is  the 
2D  Coulomb  potential  and  the  polarization  operator  in  the  limit  of  large  momentum  and 
frequency  transfer  is  given  by 

2n(gV2m) 

Tl{q,u>i)  =  p   \  ,  X,- 

u\  +  {q^/2my 

The  correction  term  is 

2  X  2  /■  .  f  f°° 

5U2{q,iuj)  =   --(2^  J        J    duji  j  d^p  j    deGM{p  -  q,i£)Glj{p,is) 

Gm(P- qi,i£  -  iu!i)D{qi,Ui) 

^  (C.6) 


As  a  reminder,  the  factor  of  2  is  due  to  the  contributions  from  two  identical  bubbles. 
Continuing  fiurther,  we  will  perform  the  e  integration  by  closing  the  contour  from  the  top 

g^^"         /*  poo  r 

6U2{q,iu)   =  j  d?qi  J  J  d'^pD{qi,uJi) 


X 


-f 


+ 


 e(-Cp-<rJ  

 e(-^p)  


(C.7) 
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As  before,  trivial  change  of  variables  is  made  to  keep  same  argument  in  the  9  function. 

6U2{q,iuj){2TTf    ^  _ 

J  J  — OO 

1 


+ 


{i^P  -  iip-qY{^l  +  i^P  -  iip-q-qi) 
1 

(wi  +  iip  -  %_,-,_g-)(a;i  +  iip  -  i^p-qi)'^ 

 1  

(z^p-_,-  -  i^p)2(wi  +  %  -  i^p^g-,) 
1 


i-i^p  +  i^p~q){uJi  +  iip  -  iip-qiY 
Corrections  from  the  vertex  and  self-energy  insertion:.  In  the  analysis  below  we  will 
make  use  of  the  following  notation: 

iy    =    i^p  —  i^p^qi 


(C.8) 


IX   =   i{eq  — -) 


m 


1 


1 


oJ\  +  iip  -  iip-q^         wi  +  iy 
1  1 


+  i^p  -  i^p-qi-q        uji  +  iy-ix' 
where  l/{uJi  +  i^p  —  i^p-^i-^)  will  be  taylor  expanded  in  the  small  parameter  ix  = 


iieq  -  ^^^)  {x  is  small  since  in  the  long  wavelength  limit  q  ^  0).  Adding  up  the 


m 

contributions  from  the  vertex  function  and  from  the  self  energy  insertion  we  obtain 


6U  =  SU1+SU2 
Sni 


1       /       iqx.qlm    \  If  {iqi.q/mf 


_\  _  1  /  [iqi.q/my  \ 
f,y)  {uj,+iy)A(i^._,-iF.v) 


L  (a;i  +  iyf  \    {i^^^^  -  i^^f  )     (a;i  +  iyf  V  {iip_^  -  i^^f 
1       /  ix{iq[.q/vfif 


/  ix{iqi.qlmy  \ 
'\iii^-n--iE^'^^) 


+ 


(C.9) 


(a;i  +  i?/)n(?ep-_q--%) 

Large  momentum  transfer  regions,  qi  '>  pp,  are  the  dominant  contributors  to  SU,  we  will 
Taylor  expand  l/(a'i  +  i^p  -  i^p-gi)  in  the  small  energy  term  p.qi/m  (where  p  f»  pp)  to 
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obtain 


5U 


8ni  r 


1 


d^pd^qi(LjxQ{-ip)D{quiOi) 
iq{.q/m 


+ 


iq[.q/m 


2iqi  .p/m 


1         /   {iqy.q/mf  \ 


+ 


(  {tqi.q/my 


+ 


The  first  term  of  the  above  equation  is  zero  since  <  qi.q  >6i=  0,  where  6i  is  the  angle 
between  vectors  qi  and  q.  The  second  term  is  also  zero  : 

iq[.q/m 


(C.IO) 


2"''term 


8m 


2igl  .p/m 
{uii  -ie,i ) 


^^(e,  —  pQ/mcos0)2J ' 


)3 /  ii^{£q  —  pq/m COS 9) 


where 


/ 


=  0. 


(C.ll) 


.i^{eq  -  pq/mcos6)     P{£q  -  pq/m  cos  9)'^. 

The  last  term  of  Eq.  (C.IO)  is  zero  as  well,  since  <  {qi-p}{qi-q)^  >ei=  0.  The  third  term  is 
the  leading  term  which  survives, 

1        /  {iqi-q/mf 


611  = 


(2 


J dFpd%d^,Q[-ip)D{quuj,) 


01 


m  f  qidqi  dui  r 


ql/m 


Hence 


n*  = 


m 


-/ 


qidqi  duj 
27r  27r 


-D{qi,ui) 


(C.12) 
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